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This report was prepared for the National Aeronautics and Space 
Administration, Langley Research Center, Hampton Virginia by Arvin/Calspan 
Advanced Technology Center, Buffalo, New York, it covers the preparation, 
conduct, and analysis of an in-flight simulation program investigating the 
flying qualities of aircraft in the final approach and flared landing task. 
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Department of Calspan. Dr. Philip Reynolds was the TIFS program manager and 
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Calspan, as well as Boeing, Lockheed, and DFVLR. 

The authors wish to acknowledge the contributions of the many indi- 
viduals who contributed to this program* 

John Ball, Michael Parrag, Scott Buethe - Safety Pilots 
Robert Gavin, Ralph Siracuse - Computer and Systems Preparation 
Edward Frantz, Thomas Franclemont - Aircraft and Systems Maintenance 
Ms. Christine Turpin - Report Preparation. 


'^actsma not slaw not FUiH) 


iii 



TABLE OF CONTENTS 


Section Pa 9 e 


1 INTRODUCTION 1-1 

2 EXPERIMENT DESIGN 2-1 

2.1 CONFIGURATION DEFINITIONS 2-2 

2.1.1 Configuration Selection Rationale 2-3 

2.1.2 Pole-Zero Patterns 2-4 

2.1.3 Configuration Parameter Selection 2-10 

2.2 DETAILED CONFIGURATION DESCRIPTIONS 2-17 

2.2.1 Introduction 2-17 

2.2.2 Longitudinal Configurations 2-21 

2.2.2. 1 Transfer Functions 2-26 

2. 2. 2. 2 Matrix - Vector Representation 2-32 

2.2.3 Lateral/Oirectional Description 2-40 

2.2.4 Thrust 2-44 

2.2.5 Feel System and Actuator Dynamics 2-44 

2.2.6 Turbulence Sensitivity 2-47 

2.2.7 Ground Effect. 2-48 

2.2.8 Bank Compensation 2-48 

2.3 EVALUATION TASKS AND PROCEDURES 2-53 

2.3.1 Evaluation Sequence 2-55 

2.4 PILOT COMMENT CARD AND RATING SCALES 2-56 

2.5 EVALUATION PILOTS AND FLIGHT SUMMARY 2-60 

3 EXPERIMENT MECHANIZATION 3-1 

3.1 EQUIPMENT 3-1 

3.2 SIMULATION GEOMETRY 3-5 

3.3 EVALUATION COCKPIT CONFIGURATION 3-5 

3.4 TIME DELAY COMPENSATION 3-5 

3.5 MODEL FOLLOWING VERIFICATION 3-7 

3.6 DATA RECORDING 3-9 

4 DATA 4-1 

4.1 INTRODUCTION 4-1 

4.2 EVALUATION CHRONOLOGY 4-1 

4.3 PILOT RATINGS 4-13 

4.4 PILOT COMMENTS 4-18 

4.5 APPROACH TIME HISTORIES 4-18 

‘ if- 


i 

4 


TABLE OF CONTENTS (Cont’d) 


Section 


5 INTERPRETATION OF THE RESULTS USING MIL-F-8785(C) AS A 

FLIGHT CONTROL DESIGN CRITERIA 5-1 

5.1 INTRODUCTION TO ANALYSIS 5-1 

5.2 ANALYSIS OF THE RESULTS 5-3 

5.2.1 Comparison with (Dp vs. n/a Requirements 5-3 

5.2.2 Comparison with 0 ^X 02 vs. c Requirements 5-6 

5.2.3 Time Domain Analysis 5-6 

5.2.4 Design of Rate Command Systems 5-14 

5.3 EFFECT OF PILOTING TECHNIQUE 5-16 

5.4 CONCLUSIONS AND RECOMMENDATIONS 5-26 

6 PREDICTIVE CRITERIA RESULTS AND ANALYSES 6-1 

6.1 TIME DOMAIN RESULTS AND ANALYSES 6-1 

6.1.1 introduction 6-1 

6. 1.1.1 Data Base Description 6-2 

6.1.2 Time Domain Results 6-3 

6. 1.2.1 Pilot Ratings 6-4 

6. 1.2.2 Discussion of Results 6-8 

6.1.3 Time Domain Analysis 6-9 

6. 1.3.1 Time Delay and Sensitivity Metric 

Development 6-9 

6. 1.3.2 Predictive Criteria Refinement 6-12 

6. 1.3.3 Application and Limitations of the Time 

Domain Predictive Criteria 6-14 

6.1.4 Time Domain Predictive Criteria Testing 6-16 

6. 1.4.1 Landing Program Results 6-16 

6. 1.4.2 Pitch Rate Program Results 6-20 

6. 1.4.3 LAHOS Program Results 6-23 

6. 1.4.4 Large Aircraft Program Results 6-27 

6. 1.4.5 NT-33A/Ames Study Results 6-31 

6. 1.4.6 SST Program and X-29 Data Results .... 6-34 

j 6. 1.4.7 Criteria Testing Summary 6-34 

6.2 FREQUENCY DOMAIN RESULTS AND ANALYSES 6-38 

6.2.1 Low-Order Equivalent Systems 6-38 

6.2.2 Pitch Attitude Bandwidth Criterion 6-39 

6.2.3 Altitude Rate Bandwidth Criterion 6-42 

6.2.4 Neal-Smith Analysis 6-46 


v 



TABLE OF CONTENTS (Cont’d) 


Section Page 


6.3 COMPARISON OF TIME DOMAIN CRITERIA AND FREQUENCY DOMAIN 

CRITERIA PERFORMANCE 6-49 

6.4 WASHOUT INVESTIGATION 6-52 

6.5 CONCLUSIONS 6-53 

6.5.1 Time Delay and Sensitivity Considerations 6-53 

6.5.2 Time Domain Criteria Considerations 6-54 

6.5.3 Frequency Domain Criteria Considerations 6-54 

6.5.4 Overall Predictive Criteria Considerations .... 6-55 

6.6 RECOMMENDATIONS 6-55 

6.6.1 Time Delay Sensitivity 6-55 

6.6.2 Time Domain Criteria 6-56 

6.6.3 Frequency Domain 6-56 

7 REFERENCES 7-1 

8 ABSTRACT 8-1 


vi 



LIST OF FIGURES 


Figure No . Page 

1 ANGLE OF ATTACK COMMAND SYSTEM 2-4 

2 PITCH RATE COMMAND SYSTEM 2-5 

3 SHORT TERM ANGLE OF ATTACK, LONG TERM PITCH RATE 

COMMAND SYSTEM 2-8 

4 SHORT TERM PITCH RATE, LONG TERM ANGLE OF ATTACK 

COMMAND SYSTEM 2-9 

5 ROLL CONTROL SYSTEM 2-42 

6 YAW CONTROL SYSTEM 2-43 

7 ROLL STEP RESPONSE, 20 DEG 6 AW STEP 2-45 

8 YAW STEP RESPONSE, \ INCH 5 RP STEP 2-46 

9 APPROACH AND LANDING TASK 2-54 

10 PILOT COMMENT CARD 2-57 

11 COOPER-HARPER HANDLING QUALITIES RATING SCALE 2-58 

12 PIO TENDENCY CLASSIFICATION 2-59 

13 USAF/CALSPAN TOTAL IN-FLIGHT SIMULATOR (TIFS) 3-2 

14 TIFS MODEL FOLLOWING SIMULATION 3-3 

15 EVALUATION PILOT CONTROLS AND INSTRUMENT PANEL 3-6 

16 MODEL FOLLOWING - FLIGHT 893, APPROACH 1, CONFIGURATION B . 3-8 

17 MODEL FOLLOWING - FLIGHT 893, APPROACH 10, CONFIGURATION 28 3-9 

18 COMPARISON WITH MIL-F-8785(C) SHORT PERIOD REQUIREMENTS 

7 OF 12 CONFIGURATIONS CORRECTLY PREDICTED 5-5 

19 CONFIGURATIONS 1-10, 13 AND 14 uv,t 92 vs. £ S p REQUIREMENTS 

4 OF 12 CONFIGURATIONS CORRECTLY PREDICTED 5-7 

20 TIME HISTORY RESPONSE ENVELOPE 5-7 

21 NORMALIZED ANGLE OF ATTACK RESPONSE, CONFIGURATION 1 . . . . 5-9 

22 NORMALIZED ANGLE OF ATTACK RESPONSE, CONFIGURATION 4 . . . . 5-9 

23 NORMALIZED ANGLE OF ATTACK RESPONSE, CONFIGURATION 5 . . . . 5-10 

24 NORMALIZED ANGLE OF ATTACK RESPONSE, CONFIGURATION 8 . . . . 5-10 

25 NORMALIZED ANGLE OF ATTACK RESPONSE, CONFIGURATION 9, 13 . . 5-11 

26 NORMALIZED ANGLE OF ATTACK RESPONSE, CONFIGURATION 16 ... 5-11 

27 SHORT TERM PREFERENCE AS SPECIFIED BY MIL-F-8785(C) .... 5-15 

28 TYPICAL APPROACH, CONFIGURATION 1, PILOT A 5-17 

29 TYPICAL APPROACH, CONFIGURATION 1, PILOT B 5-18 

30 TYPICAL APPROACH, CONFIGURATION 1, PILOT G 5-19 

31 PSD, CONFIGURATION 1, PILOT A 5-21 

32 PSD, CONFIGURATION 1, PILOT B 5-22 

33 PSD, CONFIGURATION 1, PILOT G 5-23 


vii 



LIST OF FIGURES (Cont’d) 


Figure No. 


46 

47 

48 


Page 


PSD, CONFIGURATION 14, PILOT A 5-24 

PSD, CONFIGURATION 14, PILOT G 5-25 

INDIVIDUAL PILOT RATINGS, FLARE AND TOUCHDOWN (FLIGHT 

CONTROL CONFIGURATIONS) 6-5 

INDIVIDUAL PILOT RATINGS, FLARE AND TOUCHDOWN (TIME DELAY 

CONFIGURATIONS) 6-6 

EFFECTS OF TIME DELAY AND PITCH SENSITIVITY FOR WHEEL 

CONTROLLERS 6-7 

REFINED TIME DOMAIN CRITERIA APPLIED TO THE LANDING PROGRAM 

(TIME DELAY CONFIGURATIONS) 6-17 

REFINED TIME DOMAIN CRITERIA APPLIED TO THE LANDING PROGRAM 

(ALL CONFIGURATIONS) . . . . 6-18 

REFINED TIME DOMAIN CRITERIA APPLIED TO THE PITCH RATE 

PROGRAM 6-21 

REFINED TIME DOMAIN CRITERIA APPLIED TO THE LAHOS PROGRAM . 6-25 

REFINED TIME DOMAIN CRITERIA APPLIED TO THE LARGE AIRCRAFT 

PROGRAM 6-28 

REFINED TIME DOMAIN CRITERIA APPLIED TO THE NT-33/AMES 

STUDY (NOT YET PUBLISHED). . 6-32 

REFINED TIME DOMAIN CRITERIA APPLIED TO THE SST PROGRAM 

AND X-29 UN-PUBLISHED DATA 6-35 

THETA BANDWIDTH RESULTS (REFERENCE 5 CONFIGURATIONS INCLUDED) 6-43 

ALTITUDE RATE BANDWIDTH RESULTS (CURRENT PROGRAM) 6-45 

NEAL SMITH RESULTS (PITCH ATTITUDE BANDWIDTH = 2.0 rad/sec) 6-48 


» 


i 


viii 



LIST OF TABLES 


Table No. 


1 TIME DELAY/SENSITIVITY MATRIX 2-12 

' 2 CONFIGURATION CHARACTERISTICS 2-22 

3 MEASURED EFFECTIVE TIME DELAY 2-39 

. 4 BANK COMPENSATION GAINS 2-52 

5 FLIGHT SUMMARY 2-62 

6 TIFS RECORDING LIST 3-10 

7 TIFS/ELARED LANDING APPROACH EVALUATION FLIGHT DATA SUMMARY 4-3 

8 CONFIGURATION - FLIGHT NUMBER/ElLOT CROSS REFERENCE .... 4-12 

9 PILOT RATINGS, APPROACH 4-14 

10 PILOT RATINGS, FLARE AND TOUCHDOWN 4-15 

11 PILOT RATINGS, OVERALL 4-16 

12 PIO RATINGS 4-17 

13 PILOT RATINGS SUMMARY 5-4 

14 LANDING PROGRAM - TIME DOMAIN CRITERIA APPLICATION 6-19 

15 PITCH RATE PROGRAM - TIME DOMAIN CRITERIA APPLICATION. . . . 6-22 

16 LAHOS PROGRAM - TIME DOMAIN CRITERIA APPLICATION 6-26 

17 LARGE AIRCRAFT PROGRAM - TIME DOMAIN CRITERIA APPLICATION. . 6-29 

18 T-33/AMES STUDY - TIME DOMAIN CRITERIA APPLICATION 6-33 

19 SST PROGRAM AND X-29 DATA - TIME DOMAIN CRITERIA APPLICATION 6-36 

20 REFINED TIME DOMAIN CRITERIA TESTING RESULTS 6-37 

21 LOW-ORDER EQUIVALENT SYSTEM RESULTS 6-40 

22 PITCH ATTITUDE BANDWIDTH CRITERION RESULTS 6-41 

23 ALTITUDE RATE BANDWIDTH RESULTS 6-44 

24 COMPARISON OF FREQUENCY DOMAIN AND TIME DOMAIN PREDICTIVE 

TECHNIQUES 6-50 



LIST OF SYMBOLS 


ADI 

b 

c 

CAP 

C D 

C D 0 

c Di 

c.g. or CG 

C L 

C C 0 

C Li 

c Cj 

Cn 


Cq 

c *j 

Cm 



Cni 

Cn j 


c yi 



attitude director indicator 
wing span, ft 

mean aerodynamic chord, ft 

control anticipation parameter = w§p/n z /a 

drag coefficient = D/qS 

drag coefficient at zero angle of attack 

3Co/3i, i = a, 6 e , deg _1 

center of gravity 

lift coefficient = L/qS 

lift coefficient at zero angle of attack 

3 Cl/ 3 i, i = a, 5 e , 6 DL c, deg - 1 

— 3 Cl/3 J f J = a, q, deg - 1 
c 

rolling moment coefficient = L/qSb 

3C/3i, i = B, 6 g* *5r» 9eg -1 
2V 

5 - 3C /3 j, j = p, r, deg - 1 
pitching moment coefficient = M/qSc 
pitching moment coefficient at zero angle of attack 
3Cm/3i> i = a, 6 e , deg - 1 

— 3Cm/3j» J = a, q, deg - 1 
c 

yawing moment coefficient = N/qSb 

3Cn/3i> i = B* deg 

2 v 

5“ 3C n /3j, J = P, r f deg - 1 
side force coefficient = Y/qS 

3Cy/3i, i = B, 6 a > deg - 1 

2V 

b" 3 c y / 3 J> j = P, r, deg - 1 


x 


LIST OF SYMBOLS (Cont'd) 


center of rotation 
drag, lb 

direct lift control 
digital reversion mode 

decibel units for Bode amplitude = 20 log 10 (amplitude) 

aileron wheel force, lb 

elevator wheel force, positive aft, lb 

rudder pedal force, lb, positive, right pedal forward 

gravitational constant = 32.17 ft/sec 2 

altitude of airplane c.g., ft 

altitude of airplane at pilot station, ft 

altitude of airplane at model wheels, ft 

Handling Qualities Rating 

Horizontal Situation indicator 

moments of inertia about X, Y, Z body axes, slug-ft 2 
product of inertia about X, Z body axes, slug-ft 2 
Instrument Flight Rules 

steady state pilot gain altitude loop closure, rad/ft 

steady state pilot gain in attitude loop closure, lb/rad 

loop gain in pitch augmentation system, deg/deg/sec 

command gain 

integrator gain 

loop gain 

lift, lb 

moments about X, Y, Z body axes, ft-lb 



LIST OF SYMBOLS (Cont’d) 


m 

ms 


N.D. 


n y > n z 

"j 

Nq, N a , Ny 
N ZCG 

N Z P or n z p 
P, q or Q, r 
PI OR 
PR 
PSD 

q 

9c 

s 

S 

T 

TD 

T Z 

Tf 

Two 

tl 

1/t 01 

1/ t 02 ~ La 
TIFS 


mass of airplane, slugs 

milliseconds 

normal digital mode 

lateral, normal acceleration, g’s (n y + right, n z + down) 

transfer function numerators 

transfer function numerators of q, a, v per 6 e 

normal acceleration at CG 

normal acceleration at pilot station 

roll, pitch, yaw rates, deg/sec 

pilot-induced oscillation rating 

pilot rating 

power spectral density 

dynamic pressure = ^ pv 2 , lb/ft 2 

pitch rate command 

Laplace operator, sec -1 

reference wing area, ft 2 

total thrust, lb 

equivalent time delay from equivalent systems analysis, sec 
lead + lag pre-filter zero 
lead + lag pre-filter pole 
washout pre-filter pole 

effective time delay from maximum slope intercept method, sec 

lower frequency zero in q/6 e transfer function 

higher frequency zero in q/6 e transfer function, £ ^*~ a 

V C L 

Total In-Flight Simulator 


xii 



LIST OF SYMBOLS (Cont’d) 


At 

V 

Vi 

V X > V y , V z 

w 


X, Y, Z 


XMP 

X PCR 


Yp h 

Yp e 

a 

a g 

a T 

°I 


e 

B g 

Pi 

Y 

An z 

5a 

5 AW 


rise from time history criteria analysis, sec 
true airspeed, ft/sec 
inertial speed, ft/sec 

velocity components along x, Y, Z body axes, ft/sec 
airplane weight, lb 

body axes, X-Z plane is in plane of symmetry with X directed 
forward parallel to the fuselage reference line, Z directed 
downward, and Y directed out the right wing 

distance along X-body axis between c.g. and pilot station, ft 

distance along X-body axis between instantaneous center of 
pitch rotation and pilot station, ft 

pilot describing function in altitude loop closure 

pilot describing function in attitude loop closure 

angle of attack, deg 

turbulence component of angle of attack, deg 

total angle of attack with respect to true airspeed, deg 

inertial angle of attack (with respect to inertial velocity), 
deg 

sideslip, deg 

turbulence component of sideslip, deg 

total sideslip with respect to true airspeed, deg 

inertial sideslip (with respect to inertial velocity), deg 

flight path angle, deg 

incremental normal acceleration, g’s 

aileron surface deflection, included angle positive 
left T.E. down, deg 

aileron wheel deflection, positive clockwise, deg 


xiii 



LIST OF SYMBOLS (Cont'd) 


5 AWC 

5 e 

6 EC 

«r 

fiRP 

fiRPC 

fith 

C 

^ph 


X 

P 


°i 


TL 

tr 

4 > 


“BW 


0)(j 


<*>ph 

^sp 

(i) rs 


aileron command positive left T.E. down, deg 

elevator surface deflection, positive T.E. down, deg 

elevator column deflection, positive aft, inch 

rudder surface deflection, positive T.E. left, deg 

rudder pedal deflection, positive right pedal forward, inch 

rudder, command positive T.E. left, deg 

throttle lever position, deg 

damping ratio 

damping ratio of Dutch roll mode 

damping ratio of phugoid mode 

damping ratio of short period mode 

damping ratio of coupled roll-spiral mode 

pitch attitude, deg or rad 

aperiodic real root magnitude, sec -1 

air density, slug/ft 3 

mean square gust intensity, i = a, B, deg 

time constant of pilot’s lead element, sec 

time constant of roll mode, sec 

time constant of spiral mode, sec 

bank angle, deg 

bandwidth frequency, rad/sec 

undamped natural frequency of Dutch roll mode, rad/sec 
undamped natural frequency of phugoid mode, rad/sec 
undamped natural frequency of short period mode, rad/sec 
undamped natural frequency of coupled roll-spiral mode, rad/sec 


xiv 



LIST OF SYMBOLS (Cont’d) 


c.g. 

DLC 

e 

9 

G.E. 

I 

m 

MP or PM 

MTCG 

P 

TIFS or 
unsubscripted 

sp 

WH 

m c.g. 

stab 


SUBSCRIPTS 

center of gravity 
direct lift control 

equivalent parameter from equivalent system analysis 

turbulence component 

ground effect 

inertial quantity 

model quantity 

model quantity at pilot station 
model quantity transformed to TIFS c.g. 
quantity at pilot station 
TIFS quantity at its c.g. 

short period mode 
model wheel height 
model center of gravity 
stability axis 


xv 



Section 1 
INTRODUCTION 


Much research using in-flight simulators has been conducted in recent 
years in the area of flying qualities for the flared landing approach phase of 
flight [References 1 through 5). Most of these studies dealt with obtaining 
data to verify or refine old criteria and to develop new criteria. Because of 
the introduction of highly augmented aircraft with the resulting non- 
conventional airplane response, effective time delay, and higher-order dynamic 
behavior, some flying quality criteria that concentrated on one airplane 
state, namely pitch rate, have been shown not to work well. The primary 
emphasis of this in-flight flying qualities experiment utilizing the USAFAlFS 
(Total In-Flight Simulator), was to generate a consistent set of data to find 
out what the pilot requires in order to be able to flare and land an airplane. 
Two separate areas of analysis were performed on the data. One was to 
investigate what kind of commanded response (e.g., angle of attack or pitch 
rate) and its characteristics that the pilot preferred. The other area was to 
refine a time history criterion that took into account all the necessary 
variables and their characteristics that would accurately predict flying 
qualities. The result of the first part of the program would be to provide 
guidelines to the flight control system designer in developing systems using 
MIL-F-8785(C) as a guide that would yield the responses that pilots prefer in 
flared landings. The second part of the program provides the flying qualities 
engineer with an accurate predictive tool which would tell him how good the 
resulting system would be. 

Flight Control Design Configurations 

The flying qualities specification published in 1969, MIL-F-8785(B) 
(Reference 6) represented the culmination of nearly 20 years of experimental 
flight research involving many variable stability aircraft such as the 
USAF/Calspan NT-33A, the Princeton Navion, and the NASA/Boeing 367-80. This 
specification for the first time quantitatively defined satisfactory and 
acceptable regions of specific modal parameters of an airplane, such as short 
period, phugoid and Dutch roll frequencies and damping ratios. A more recent 
revision, MIL-F-8785(C), (Reference 7) for the most part preserves the modal 
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requirements specified in MIL-F-8785(B) without significant revision both in 
terms of the modal requirements and definitions of the modes. MIL-F-8785(C), 
for instance, states that " short period response of angle of attack shall meet 
the requirements . The proposed MIL Standard and Handbook (Reference 8), 
however, is oriented almost totally in the direction of defining an 
"equivalent** short period mode from the pitch rate rather than the angle of 
attack response of the vehicle. If the vehicle is unaugmented or if the 
response order has not been increased by compensation or other dynamic ele- 
ments, then it does not matter whether the short period is defined from the 
angle of attack or from the pitch rate behavior of the vehicle. However, if 
dynamic elements have been introduced in the control system, such as a pre- 
filter or integral plus proportional compensation in the loop with pitch rate 
feedback, then it is unlikely that the short period frequency, equivalent or 
otherwise, may be accurately obtained from the pitch rate response of the 
vehicle. 


The flying qualities specification, MIL-F-8785(C) and the preceding 
(B) version, unfortunately have seldom served as direct criteria for the 
design of a flight control system for the longitudinal-vertical degrees of 
freedom of motion of an aircraft. There appears to be two primary reasons for 
thisi 

1. The effect of compensation networks and other dynamic elements 
that produce higher order response effects in the flying quali- 
ties frequency range of interest were not explicitly addressed 
in the specification. The implicit message of the specifica- 
tion is that higher order response behavior almost always 
degrades the flying qualities of the vehicle. 

2. The flight control system designer is oriented more towards the 
command input-controlled output response philosophy. Generally 
speaking, the flight control system designer has often con- 
sidered the feedback quantities to be the "controlled 
variable(s)" and seeks design criteria based upon the response 
of the "controlled variables" to a command input. The modal 
approach of MIL-F-8785(C) addresses dynamic requirements in 
terms of commanding the entire airplane rather than a specific 
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state variable(s) of the airplane. Although the command- 
response approach to control system design is not incompatible 
with the requirements, the tendency has been to generate new 
criteria based upon the feedback rather than the modal para- 
meter requirements. Several significant examples of this 
approach to control system criteria exist. The most prominent 
among them are C*(t) (Reference 9), the Shuttle Orbiter pitch 
rate envelope criteria (Reference 10), and the angle of attack 
time history response envelope (Reference 11). 

There has been a profusion of alternate flying qualities criteria 
since the publication of MIL-F-8785(B) . Criteria are open-loop, closed-loop, 
man-in-the-loop and both frequency domain and time domain oriented. Some 
attempt directly to transform from the modal specification format of MIL-F- 
8785(C) while others are empirically derived from the analysis of sets of 
flight data. Although most are successful to some degree, they are primarily 
compliance or evaluation methods rather than direct design criteria that can 
aid the designer in trying to decide what to feedback, how much and whether 
higher order networks such as an integral plus proportional compensation 
should be added to the system. All of these criteria can be useful, but most 
do not appear to have the objective of direct extensions to the solid foun- 
dation of MIL-F-8785 (C) or interpretation of MIL-F-8785(C) for the flight 
control system designer. It would seem to be a desirable ultimate objective 
to help the control system designer make direct use of the results of the 
twenty years of experimental flight test that are embedded in the flying 
qualities specification. 

It is an objective of the experimental flight program described in 
this report to interpret the flying qualities requirements of MIL-F-8785(C) in 
terms of command/response configurations. In this sense, then, the effort 
i described in this report is directed toward the first steps required to trans- 

form MIL-F-8785(C) from a pure flying qualities specification into a joint 
“ flying qualities/flight control design criteria. The two commanded longitu- 

dinal state variables chosen for investigation are pitch rate and angle of 
attack, the two "natural** states of constant speed aerodynamic flight. By 
evaluating pitch rate command and angle of attack command configurations a few 
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basic guidelines are established to help the flight control system designer 
decide whether the system should behave as a conventional aircraft (angle of 
attack command system) or as a rate command, attitude hold system. Rate 
command and angle of attack command configurations were specified totally 
within the context of MIL-F-8785(C) in order to determine pilot preference for 
one or the other as a function of l/x^ and phugoid mode characteristics. The 
results suggest that pilot preference for angle of attack or pitch rate 
command is a function of 1 /t 02 . 

Time Domain Criteria (Pitch Axis, Flared Landing Task) 

The primary objective of this phase of the experiment was to refine 
the time domain predictive criteria of Reference 5. (This criteria is appli- 
cable to the pitch axis in the flared landing task.) A major limitation of 
the original criteria was an inability to contend with time delays. An 
additional limitation was concern as to the ability to handle excursions of 
short period frequency that fell outside the borders of Level 1 boundaries of 
MIL-F-8785 N z /a vs. Wn-.. requirements. 

Experience with the Calspan Learjet had indicated that time delay 
effects were strongly dependent on control sensitivity (sensitivity defined, 
in this instance, as maximum slope of the pitch rate response per pound of 
pitch controller input). Consequently, it was determined that a valid time 
delay experiment must include sensitivity effects and these sensitivity 
effects might be applicable to high frequency short period cases from the 
MIL-F-8785 requirements. 

A time delay test matrix was chosen which utilized two baseline 
flight configurations that were Level 1 in flying qualities performance (one a 
"conventional" angle of attack command, and the other a pitch rate command). 
The matrix was completed by applying three levels of pitch sensitivity and 
three levels of time delay to the baseline configurations. 

The in-flight data provided by this matrix resulted in the desired 
time delay and sensitivity matrices for criteria refinement. The resultant 
refined criteria was applied to applicable configurations of this experiment 
as well as to the flight data of the pitch rate program (Reference 5), the 
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LAH05 program (Reference 1), the Large Aircraft program (Reference 4), the 
Ames NT-33 Study (unpublished), the SST program (Reference 12), and the X-29 
program (unpublished), for a total of 129 configurations. These 129 con- 
figurations consisted of variations in short period frequency and damping, 
angle of attack command, pitch rate command, 1A 02 variations, time delay 
variations, prefilter variations, sensitivity variations, wheel and stick 
controllers, and aircraft sized from the X-29 to 1,000,000 lb. gross weight. 

The results were that the time domain predictive criteria predicted 
60% of the 129 configurations from seven different programs to within one 
Cooper-Harper rating, 88% to within two HQR ratings, and 81% by flying quali- 
ties Level. The predictions for this subject program, where seven pilots and 
many repeated evaluations provided a more valid statistical base, were* 77% 
within 1 HQR, 96% within 2 HQR’s, and 100% by flying qualities Level. 

A number of frequency domain predictive techniques were applied to 
the data of this program and to selected configurations from other programs. 
The most promising of these was a technique using Neal-Smith pilot lead 
compensation angle as the prime parameter. For this program the results were* 
65% predicted within +1 HQR, 81% within 2 HQR and 73% by flying qualities 
level, however, 19% missed by more than 2 HQR’s (compared to 4% for the time 
domain criteria) and 4 Level 3 configurations were missed (compared to none 
missed by the time domain criteria). This technique did not successfully 
account for time delay and sensitivity effects and needs to be applied to a 
wider data base. 

The balance of the report is organized in two volumes as followsi 

Volume I 

Section 2 Experiment Design - includes the detailed con- 

figuration descriptions and evaluation procedures. 

Section 3 Experiment Mechanization - describing the TIFS 

aircraft and equipment. 

Section 4 Data - includes the flight chronology and raw data 
obtained in the experiment (with some references to 
the Appendices). 
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Section 5 Interpretation of the Results Using MIL-F-8785(C) as a 
Flight Control Design Criteria. 

Section 6 Predictive Criteria Results and Analysis - includes 
time domain and frequency domain analyses. 

Sections 5 and 6 each contain their own conclusions 
and recommendation subsections. 


Section 7 References 


Volume II 

Appendix A 
B 
C 
D 
E 
F 


Pitch Step Responses 

Throttle Step Responses 

Model Following Verification Steps 

Pilot Comments 

Approach Time Histories 

Neal-Smith Parameter Plane Plots 
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Section 2 
EXPERIMENT DESIGN 


The motivation for the flight experiments performed using the 
USAF/AFWAL Total In-Flight Simulator (TIFS) and discussed in this report are 
many. The flying qualities specification MIL-F-8785(C) presents flying quali- 
ties requirements in terms of satisfactory and acceptable regions of specific 
modal parameters of an airplane such as short period and phugoid frequency and 
damping ratio. Modal residues (or zeros of transfer functions) are not speci- 
fically addressed, therefore, the responses between aircraft having the same 
modal characteristics can vary over a wide range. 

This experiment explores, in a preliminary way, a range over which 
the response of the aircraft may vary yet satisfy flying qualities require- 
ments as defined by MIL-F-8785(C). The range of dynamic behavior is related 
to flight control system design considerations in the sense that purely angle 
of attack command and purely pitch rate command configuration are compared in 
as directly a form as possible. For instance, all the configurations were 
designed to have a short period natural frequency of 2 rad/sec. For the angle 
of attack configuration, the short period damping ratio was held constant at 
? S p = 0.7. The pitch rate command configuration required, because of the 
''command*' requirements, variations in £ sp as a function of 1 /t 02 . A more 
complete discussion of the rationalization for the selection of these con- 
figurations is given in Reference 13. 

The flight control system designer is often oriented toward the 
command/response philosophy of control system design, so the experiments are 
designed to provide, in a preliminary sense, an interpretation or guideline 
for the control system designer in the sense that he will be able to decide 
whether to design a control system that will result in an angle of attack 
command (or conventional aircraft) responding system or a pitch rate command, 
attitude-hold system. In this sense, the attempt is to demonstrate to the 
flight control system designer that not only are the specifications in 
MIL-F-8785(C) an appropriate way to judge compliance with flying qualities 
requirements, but it may also be directly used as a flight control system 
design criteria. In this way, the results of the twenty years of experimental 
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flight tests used to define the flying qualities specification can also be 
directly used by the flight control system designer. 

2.1 CONFIGURATION DEFINITIONS 

The experiment described below considers angle of attack and a pitch 
rate command systems in terms of the locations of the poles of the system with 
respect to the zeros in various transfer functions. In this respect, the 
intent is to try to determine whether the modal parameter approach as spe- 
cified in MIL-F-8785(B) has been properly interpreted in the MIL-F-8785(C) 
standard in terms of aircraft axes. The short period and the phugoid will be 
considered separately because it is possible to design a pitch rate command 
system for the short term, but an angle of attack command system in the long 
term, or vice versa. The idea is to try to determine pilot preference both in 
the short term and long term. 

The angle of attack or the pitch rate command systems can be defined 
solely in terms of the locations of the vehicle poles with respect to the 
zeros of the transfer functions. A pitch rate command, attitude hold system 
will produce a pole-zero cancellation such that three poles are placed at the 
zeros of the transfer function located at the origin of the s plane, at 
-lA 0 p and at -1 /tq 2. Therefore, the response in pitch rate is dominated by 
the one remaining pole. In the angle of attack command system two poles are 
located at the low frequency zeros <%, of the angle of attack transfer 
function. The response is dominated by the remaining two poles, which define 
the short period mode. These systems are briefly described below without 
regard to how they may be mechanized. The mechanization problem is not a 
difficult one and is discussed in Section 2.2.1 and Reference 13. 
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2.1.1 Configuration Selection Rationale 


For this experiment, command configurations such as pitch rate, angle 
of attack, or rate of change of flight path were designed with direct refer- 
ence to MIL-F-8785(C). The configurations were defined as if the particular 
quantity commanded or fed back was constrained by the requirements of the 
modal parameters of the MIL-F-8785(C) format. Consider the sketch below, 
which shows the root locus plot for a two degree of freedom vehicle represen- 
tation using pitch rate feedback. 



PITCH RATE FEEDBACK - 2 DOF VEHICLE REPRESENTATION 

As the feedback gain is increased, the damping ratio of the short 
period mode becomes high and eventually the short period mode separates into 
two real roots; one that will terminate at the zero located at s = -1 /t 02 and 
the other at infinity if not constrained . The short period mode is defined by 

(s + p,)(s + p-) and indicates a damping ratio greater than critical with short 

( Pi + P2 ) 

period frequency of u) sp = /pjp 2 and damping ratio of c s /p^ 

Therefore, in the limit one of the short period poles will be located 
at s = -l/r 02 and the pitch rate response will be a pure first order response 
characterized by the transfer function q/P s (s) = K/(s + p^), even though the 
short period mode itself is characterized by (s + p^Ks + 1/t 02 ). 
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2.1.2 Pole-Zero Patterns 


The sketches below show the pole-zero pattern representation of the 
different “controlled variable” configurations. 

Angle of Attack Command System (Configurations 1, 5) 

POLE-ZERO CONFIGURATION RESPONSE TO STEP COMMAND 




ANGLE OF ATTACK POLE-ZERO ANGLE OF ATTACK RESPONSE 

CONFIGURATION 




PITCH RATE POLE-ZERO PITCH RATE RESPONSE 

CONFIGURATION 


Figure 1. Angle of Attack Command System 

As shown in Figure 1, the response of the angle of attack command 
system is dominated by the short period poles u) S p, s S p. The phugoid poles are 
located at the low frequency zeros q, of the a/6p(s) transfer function. 
The result can be a quick, smooth and well behaved angle of attack response as 
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defined by the short period mode. Theoretically there is no residue in the 
angle of attack response in the phugoid mode; i.e., a(t) = 0 after the short 
period response. 

The pitch rate response of the angle of attack command system is 
typical of a conventional aircraft. The transfer function zero at -I/t^ 
produces an overshoot in the pitch rate response to a step command input, and 
a significant phugoid mode residue with zero ultimate steady state value is 
evident. 

Pitch Rate Command System (Configuration 2, 6) 

As the angle of attack command system showed pole-zero cancellation 
in the angle of attack transfer function, the pitch rate command system 
indicates pole-zero cancellation in the pitch rate transfer function. The 
pole-zero pattern showing these cancellations are indicated in Figure 2 below. 


* 




ANGLE OF ATTACK POLE-ZERO ANGLE OF ATTACK RESPONSE 

CONFIGURATION 



PITCH RATE POLE ZERO 
CONFIGURATION 



Figure 2. Pitch Rate Command System 
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The response of pitch rate in the pitch rate command system is 
dominated by the single order pole p^ shown in the above figure. The response 
shows no residue in the phugoid mode, and the zero at the origin is cancelled 
by a pole, which indicates that the system will be an attitude-hold system. 
The angle of attack response is generally sluggish, dominated by the poles at 
- 1 /t 0 2 and -l/t ei that are not cancelled by numerator zeros of the angle of 
attack transfer function. The pole at the origin also contributes to the 
response and leads to a steady state ramp response in angle of attack. After 
the step input is returned to zero, the pitch rate returns to zero but the 
change in angle of attack does not. The pitch rate command, attitude hold 
system is also an angle of attack M hold M system, although the response in 
angle of attack is normally so sluggish that steady state angle of attack 
would likely be rarely seen in actual flight. Speed change will also exhibit 
neutral stability. 

The two types of system described in Figure 1 and Figure 2, mainly a 
pitch rate command and angle of attack command system involve both short 
period and phugoid dynamic behavior of both of the response variables. Each 
part contributes significantly to the dynamic behavior of the system. The low 
frequency behavior of the angle of attack command system is such that after 
the angle of attack reaches steady state, then changes in flight path are 
equal to changes in pitch angle* i.e., by = Ae since di(t) = 0. The pilot can 
judge changes in flight path by observing changes in pitch angle. Because 
pitch rate eventually goes to zero following a step command, the pitch atti- 
tude and the flight path reach new steady state values. In the pitch rate 
command, attitude-hold system, the angle of attack responds sluggishly and 
never reaches a steady state value to a step command input. The change in 
flight path angle is not equal to changes in pitch angle* i.e., Ae / by and 
the pilot has more difficulty in judging changes in flight path by observation 
of changes in pitch angle. The result of the sluggish angle of attack 
response is frequently an overcontrol tendency by the pilot during flare and 
landing. A correction of the overcontrol leads to pilot complaints of 
"non-monotonic" stick forces. 
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The differences in the short period response are more obvious. In 
the angle of attack command system, the numerator zero in the pitch rate 
transfer function may be considered a lead term in the pitch rate response. 
In the pitch rate command system, the singularity that previously was a pitch 
rate lead becomes a pole or lag in the angle of attack response. 

Hybrid Systems 

Simple variations in the types of pure controlled variable systems 
designed for this experiment should allow both the flying qualities engineer 
and the flight control system designer to determine whether or not the 
“controlled variable** philosophy of control system design applies to both the 
short term and the long term or phugoid mode. For instance, by the indepen- 
dent placement of the short period and phugoid poles, it is a relatively 
simple matter to obtain a short term angle of attack command, long term pitch 
rate command system. This can be done as shown in Figure 3 below, in which 
the short period poles are placed as if the system were angle of attack 
command, while the phugoid poles are placed as if the system were pitch rate 
command. The converse, as shown in Figure 4 below, can also be accurately 
evaluated using the TIFS variable stability airplane. 

In the past it has been often stated that the pilot is little 
affected by the long term or phugoid motion of the vehicle. It has been 
assumed that the pilot either ignores these long term effects or corrects for 
them more-or-less subconsiously. If this hypothesis is true, then it should 
make no difference if the phugoid poles were located at either the zeros of 
the numerator of the angle of attack transfer function (u> a , or at the 
origin and at -l/t ei , two of the numerator zeros of the pitch rate transfer 
function. It is expected that the hybrid variations depicted by Figures 3 
and 4 whould help significantly to settle the question of the importance of 
phugoid dynamics with respect to flying qualities. 
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Short Term Angle of Attack - Long Term Pitch Rate (Configuration 3, 7) 


The short-term angle of attack, long-term pitch rate command system 
pole-zero configuration is shown in Figure 3. 



ANGLE OF ATTACK POLE ZERO 
PATTERN 



ANGLE OF ATTACK RESPONSE 
TO STEP COMMAND 




PITCH RATE POLE-ZERO PITCH RATE RESPONSE 

PATTERN TO STEP COMMAND 


Figure 3. Short Term Angle of Attack, Long Term Pitch Rate 
Command System 


The behavior of this system is characterized by the smooth and rapid 
short period angle of attack response and the modal residues of the poles 

located at 1 /t 6 ^ and at the origin. The pitch rate response is characterized 
by an initial pitch rate overshoot followed by a steady state pitching rate; 
no phugoid mode residue is evident. 
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Short Term Pitch Rate - Long Term Angle of Attack Command System 
(Configuration A, 8) 


The short term pitch rate, long term angle of attack system is shown 
in Figure A below* 



ANGLE OF ATTACK POLE-ZERO 
PATTERN 



ANGLE OF ATTACK RESPONSE 
TO STEP COMMAND 



PITCH RATE POLE-ZERO 
PATTERN 



Figure 4. Short Term Pitch Rate, Long Term Angle of Attack 
Command System 


The behavior of this system is characterized by an angle of attack 
response dominated by the pole at -1/t 02 and can be sluggish. The angle of 
attack response remains steady in the long term. The pitch rate response is 
initially rapid and dominated by the single pole at but then exhibits the 
effects of a significant residue at the phugoid mode frequency. 
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2.1.3 Configuration Parameter Selection 


All configurations with one exception, Configuration 12, were 
deliberately designed with the pilot located at the center of rotation (i.e., 
percussion) in order to eliminate the initial acceleration experienced by the 
pilot during maneuvering when located somewhere other than the center of 
rotation. For all configurations the short period natural frequency was held 
constant at = 2 rad/sec, while 1 A 01 was held constant at 1/t 01 = 0.1 
rad/sec and the high frequency zero of the a/F s (s) transfer function was 
located for ail configurations at 1 /t = -3.7. This value was relatively low 

but yielded a vehicle configuration having the pilot and center of rotation 
co-located . 

The first four configurations 1-4 were chosen to have a value of 
1 /t 6 2 = 0.5 and the low frequency numerator zeros of the angle of attack 
transfer function were defined to have a frequency Wq of 0.3 rad/sec with 
S = 0.10. For the pitch rate command configuration, the result was a very 
rapid response in pitch rate dominated by a pole at s = - 8.0 while the angle 
of attack response was sluggish, dominated by the pole at s = -0.5 (or 

-IA 92 ). The zeros of the angle of attack transfer function meant that large 
phugoid mode residues would appear in the attitude and flight path response 
for the angle of attack command configurations. The purpose of these designs 
were to test the following hypothesis* 

1. Although the pitch rate response could be very rapid, the 
vehicle may not be Level 1 if the angle of attack response in 
the short term is sluggish. The purpose is to demonstrate that 
the MIL-F-8785(C) short period requirements pertain to the 
angle of attack response of the vehicle. 

2. The phugoid mode may be relatively high frequency and low 
damping (u> p h = 0.3, c p h = 0.1) and still be acceptable if the 
vehicle is configured as an angle of attack command system. 

The large phugoid mode residues appearing in the attitude and 
flight path responses are acceptable to the pilot because the 
angle of attack response is well behaved in the long term, a(t) 
exhibits no phugoid mode residue and, therefore, Ae = Ay in the 
long term. 
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Configurations 5-8 were designed to have a value of 1/t 02 = 0.9 and 
low frequency angle of attack transfer function zeros defined by <% = 0.1, 

5 = 0.1. The short term pitch rate response of the pitch rate command con- 
figurations would be slower than those of Configuration 2 or 4, but the angle 
of attack response would be more rapid in the short term. The lower phugoid 
mode frequency located at a frequency of 0.1 rad/sec is identical to the value 
of 1/t 01 . Therefore, the phugoid mode residue in the pitch attitude and 
flight path angle responses would be small for all the Configuration 5-8. The 
purpose of these selections were to test the following hypothesis! 

1. A pitch rate command configuration would be rated Level 1 if the 
angle of attack response satisfied the vs n/a requirement of 
MIL-F-8785(C). If true, the flying qualities of a pitch rate 
command system would be a direct function of the value of 1/t 02 . 

2. The flying qualities of a pitch rate command system would 
improve if changes in attitude more closely corresponded with 
changes in flight path angle in the long term. This is mani- 
fested by a smaller long term residue in the angle of attack 
response. 

Configurations 9 and 10 were designed to have no phugoid residue in 
either the pitch rate or angle of attack response in the long term. The low 
frequency zeros of the a/F s (s) were given values equal to those of the qA" s (s) 
transfer functions, i.e., s = 0, -l/x 0 ^. Therefore, the system would be both 
angle of attack and pitch attitude “hold” in response to pilot commands. This 
kind of configuration is feasible with an additional control surface such as a 
canard surface with relatively low power and low bandwidth characteristics. 
The neutrally stable phugoid response would appear only in the speed change 
degree of freedom of motion. The purpose of these configurations is to test 
the following! 

1. The pilot would find quite acceptable, and even prefer a 
vehicle that exhibited no phugoid mode behavior in both angle 
of attack and pitch rate. 
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2. By eliminating the phugoid mode to stick command inputs, the impor- 
tance of the phugoid mode to flying qualities can be demonstrated. 

Configurations 11 and 12 were defined as y C g command configurations 
with respect to the center of gravity of the aircraft but with the center of 
rotation located either at the pilot station (34 ft forward of the eg) 
(Configuration 11) or at the center of gravity (Configuration 12). Both 
configurations have y transfer functions with two zeros at the origin of the 
y/F s (s) transfer functions, so the denominator contained a s 2 = 0 term. The 
response in y C g for the configuration in which the center of rotation was 
located at the pilot station (Configuration 11) was non-minimum phase, with a 
transfer function zero at s = +3.86. A pole was, therefore, placed at this 
value reflected about the jio axis, i.e., at s = -3.86. This is the pole value 
that would minimize the non-minimum phase response effects and is the value 
that would approximately be achieved if the design were accomplished using 
linear optimal control methods. Because the y/F"s(s) transfer functions were 
not rational, i.e., the numerator and denominator polynomials had the same 
number of singularities, an addition pole was added at s = -2 to yield a 
rational response to a pilot command input. The purpose of these config- 
urations was to demonstrate that 

1. The s 2 = 0 terms, which dominate the long term response, are 
very detrimental to flying qualities even if the n z (t) (and in 
this case q(t) also) are smooth and well behaved in short term. 

2. The response closely resembles the behavior of a system with a 
transfer function of a/F s (s) = K/s, which is very detrimental 
to flying qualities. 

3. A "command*' configuration should be designed independently in 
the short period and long period dynamics of the vehicle. 
Although it is believed that very good Level 1 y(t) command 
configurations can be designed, the design should likely be for 
the short term only, with the long term most effectively 
designed as an angle of attack command system, i.e., phugoid 
poles at or near the values of 1 %, zeros. 
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For Configurations 13 and 14, a third value of 1 /t 02 =2.0 was chosen 
as a large value of 1/t 02 equal to the value of the short period frequency and 
resulting in a value of CAP that would place the configurations in the Level 1 
area for the Category C precision requirement but in the vicinity of the Level 
1 - Level 2 boundary with respect to Category A precision requirements. A 
large value of 1/t 02 was chosen for several reasons* 

1. To determine where the lower constant CAP boundary was accur- 
ately defined for both approach and for flare and landing. 

2. To try to help settle the controversy of whether the constant 

CAP, top vs n/a requirement of the MIL-F-8785(C) specification 
or the vs hypothesis of the MIL-F-8785 Handbook and 

MIL Standard more accurately defines approach, and flare and 
landing requirements. 

3. To test the hypothesis that a pilot might favorably accept a 
dynamical configuration in which the flight path and attitude 
responses were in close coordination during the short period 
response of the vehicle. Not only did the higher 1/t 02 reduce 
the short period steady state angle of attack change required 
to maneuver, but also made the pitch rate and angle of attack 
responses dynamically more similar in the short term. 

It was felt that the configurations designed and flown in this program 
would go a long way toward helping the control system designer properly 
interpret the MIL-F-8785(C) requirements. Most of the configurations are 
designed to test the principle that* 

1. The short period requirements of MIL-F-8785(C) apply to the 
angle of attack response of the vehicle. 

2. In the long term the system should be designed as an angle of 
attack command system. The phugoid poles should be placed at 
or near the low frequency zeros of the a/F s (s) transfer 
function . 
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3. If these objectives are achieved, then the pilot can use the 
pitch response of the vehicle as an effective surrogate for the 
flight path response because the vehicle then flies propor- 
tional to the direction it is being pointed by the pilot. 

u> Verification ( Configurations 15,16) 

5 P 

Two angle of attack command configurations with different short 
period frequencies were added as part of the experiment to provide extra data 
to further refine the time domain criteria. These configurations were similar 
to Configuration 1 except <D sp = l rad/sec for Configuration 15 and (D gp = 3 
rad/sec for Configuration 16. 

Washout Investigation (Configurations 17-20) 

The previous TIFS/^itch Rate program (Reference 5) briefly investi- 
gated the effects of a washout prefilter on a specific pitch rate command 
configuration (Shuttle Orbiter). It was decided to make a systematic investi- 
gation of washouts in this experiment. Configuration 1-2-2 of the Reference 5 
experiment was chosen as a baseline Level 2 pitch rate command configuration. 
This was called Configuration 17 in this program. A washout prefilter was 

added to the pilot command path with various washout time constants 
(Configuration 18-20). These configurations differed from all others in that 
they were mechanized with a q feedback path and a proportional-plus integral 
in the command path. 

Time Delay Control Sensitivity Matrix Configuration Description 

The purpose of this part of the test matrix was to determine the 
effect, if any, of pitch sensitivity on time delay flying qualities effects. 
The basis of the matrix was Configuration B, a Level 1 conventional angle of 
attack command system (see Appendix A for step responses). 

Three variations in pitch sensitivity were chosen during the calibra- 
tion test flights. The mid sensitivity was 0.42 deg/sec 2 /lb. This value had 
been found to be a near optimum sensitivity of wheel controllers in previous 
programs and was verified as such during the calibration flights. Minimum and 
maximum sensitivity values were chosen, by flight tests, that would still 
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yield Level 1 or borderline Level 1 flying qualities. The minimum value 
chosen was 0.25 deg/sec 2 /lb and the maximum value was 0.63 deg/sec 2 /lb. These 
values were selected by controlling the command gain of the flight control 
system. 


There are a number of ways to change pitch sensitivity, i.e., 
changing command gain, changing short period frequency, changing short period 
damping ratio, the addition of prefilters, etc. It was felt then no matter 
what the method of changing sensitivity, the result is much the same to the 
pilot. The command gain method of changing sensitivity was chosen as it 
tended to better isolate the effects of sensitivity while keeping other criti- 
cal factors constant. 

The minimum time delay of Configuration B as implemented in TIFS was 
150 ms (measured from the time of wheel force application to maximum slope 
intercept of the resultant TIFS pitch rate response). Additional transport 
time delays of 100 ms and 200 ms, respectively, were added to provide three 
levels of time delay; 150 ms, 250 ms, and 350 ms. 

The three levels of sensitivity were used with the three levels of 
time delay to obtain the matrix of Table 1. 

There was also interest in observing the difference, if any, in time 
delay effects between conventional angle of attack command and pitch rate 
command flight control systems. A Level 1 pitch rate system, Configuration 
17+LL, was chosen as a baseline for an additional 3x3 matrix. During the 
early part of the program the 350 ms time delay portion of the matrix was 
flown and the results were so similar to those of the angle of attack command 
matrix that the balance of this pitch rate command matrix was not flown. 
Table 1 also shows those three high time delay configurations. 
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Table 1 

TIME DELAY/SENSITIVITY MTRIX 


CONFIGURATION 

NUMBER 

FLIGHT CONTROL 
SYSTEM 

SENSITIVITY 

(q/lb) 

TIME DELAY 
(ms) 

B 

B 

0.42 

150 

21 

B 

0.42 

250 

22 

B 

0.42 

350 

23 

B 

0.25 

150 

24 

B 

0.25 

250 

25 

B 

0.25 

350 

26 

B 

0.63 

150 

27 

B 

0.63 

250 

28 

B 

0.63 

350 

17+LL 

17+LL 

0.42 

150 

22A 

17+LL 

0.42 

350 

25A 

17+LL 

0.25 

350 

28A 

17+LL 

0.63 

350 


Throttle Response 

The matrix of control effectiveness terms for this experiment was 
chosen such that a throttle term appeared only in the V equation, in other 
words, it was assumed that the thrust acted through the vehicle center of 
gravity aligned with the longitudinal axis of the airplane. The responses of 
the vehicle variables to a throttle command then depended only upon the aero- 
dynamic coupling between the speed change degree of freedom of motion and the 
other degrees of freedom. The effect of the different flight control system 
designs is to produce near-decoupling of the speed response from the remaining 
states for some of the configurations. Configurations 9 and 10, for instance, 
are completely decoupled dynamically j a throttle input produces only a speed 
change with no change in attitude or flight path angle. 
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2.2 DETAILED CONFIGURATION DESCRIPTIONS 

2.2.1 Introduction 

As previously mentioned, most of the longitudinal characteristics of 
the configurations were defined by exact pole and zero placement and not by 
defining an aerodynamic model with various control systems. Once specific 
poles, zeros, and gains were chosen for each transfer function of a con- 
figuration, a phase variable canonical transformation was performed on these 
transfer functions to convert them to linear, state space matrix-vector repre- 
sentations which defined the longitudinal equations of motion of the aircraft. 
The configurations which were based on specific aerodynamic configurations 
were first linearized and then linear F and G matrices were obtained for them. 

The Configurations 1-14 were first defined in the frequency domain as 
transfer functions. It was a relatively easy matter to define the pole-zero 
cancellation transfer function and the short period frequency constraints of 
the configurations using this method. The set of transfer functions for a 
particular configuration were then assembled as a transformation associated 
with a phase variable canonical form, with the transformation defined, the 
system of transfer functions were transformed into the time domain for mech- 
anization in the TIFS computer in the familiar state space or matrix-vector 
format 

x(t) = F x(t) + G u(t) 

Equation (1) can be transformed into the phase variable canonical form 

y(t) = A y(t) + B u(t) 


using a linear transformation 

x(t) = T y(t) 


(3) 
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in which the matrices A and B are of the form 


010 0 


0 

0010 0 

1 

1 

B = 

1 

1 

1 

o 01 


0 

; a o - a i Vi 


1 


( 4 ) 


The coefficients of the matrix A are determined from the character- 
istic polynomial of the equations of motion 

1 1 s-F l =s n + a_ - n ' 1 


n-1' 


ou s + 

1 o 


(5) 


The transformation matrix T is obtained from the relationship 

TS = [IS-F] adJ G cs) 

where [ls-F] ad j represents the ad jugate of the matrix [is-F] and S is the 
column matrix S T = [1 s s 2 .... s n_1 ] . 

A row of the matrix T, therefore, is composed of the coefficients of 
the numerator polynomial of a transfer function of the system. If the system 
of Equation (1) is completely controllable and observable, the square trans- 
formation matrix T is non-singular and the inverse of T exists. The transfor- 
mation from the transfer function or phase variable form of Equation (2) and 
(3) to the state space form of Equation (1) is simply given by 

x(t) = T A T _1 x(t) + T B u(t) 

= F x(t) + G u(t) 

Systems with multicontroller inputs, such as a stick and throttle can 
be easily accomodated by defining a separate transformation Ti for each 
controller u^. The original system of Equation (1) is written as 



For each input u* a separate transformation 7* is defined repre- 
senting the transfer functions of the states for a particular input Ui. A 
total of p transformations are defined from the transfer functions of the p 
inputs. The transformation from these transfer functions to the state space 
equations is then defined by 

P P 

x(t) = ( J T) A( J T) xCt) + J l B u x + T 2 B u 2 ... + T B u 

1=1 U1 ( 9 ) 

Therefore, transfer functions with respect to a throttle command 
input can be defined independently of those with respect to a stick command. 
The only requirements (other than controllability and observability) is that 
the poles of the stick and throttle transfer functions be in common. 

Example 


A simple example of the use of the phase variable transformation is 
given below. Assume it is desired to find the augmented equations of motion 
that would yield the following transfer functions to an elevator and a direct 
lift surface 


M (s + 1 /t ) 

q (s) 5 e 9 2 

s ! + 2t»n s ♦ 



M 


6 


e 


s 2 + 2i;u) n s + ojn 


From these transfer functions, the matrices A and T are defined 


then 



“ 

■ 


" 

A = 

0 

1 

T = 

Mr A. Mr 

6 e e 2 6 e 



- w n ' 


Mr 0 

°e 


• 

q 

-1 

q 


0 


= T A T 


+ T 


a 


a 


1 
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or 


• 

q 


UA e2 -200 n ) + 2 


q 

+ 

'V 

• 

a 


-l - 1 /t s 2 


a 


0 



Not only can the transformation described above be used to find the 
equation of motion equivalent to a set of transfer functions, but as described 
in Reference 14, a transformation similar to the phase variable transformation 
be used to define filters or observers. Those observers are then used to 
define a control law for any aircraft that can place all the poles of a system 
in accordance with the requirements of MIL-F-8785(C) using a minimum number of 
sensors and without increasing the order of the closed-loop response to a 
pilot command or other input to the system. 

Though most of the configurations were defined originally by poles 
and zeros, the general characteristics of all of the configurations were those 
of a medium transport aircraft. The lateral/directional characteristics were 
those used in a previous TIPS program (Reference 5). A transport type wheel/ 
column and rudder pedal feel system were also used. 

Command gains were chosen during the checkout phase of the program to 
yield initial pitch sensitivities that were typical of a transport aircraft. 
It turned out that a nominally good value for the pitch sensitivity yielded a 
maximum pitch acceleration of approximately 5 deg/sec 2 for a 10 pound step 
input. Command gains were adjusted on all of the configurations (except those 
which specifically had sensitivity variations as an experiment variable) to 
yield this nominal initial sensitivity. It should be noted that for con- 
figurations with large pitch rate overshoot the constant initial pitch accel- 
eration resulted in these configurations having higher steady state forces in 
a pitch maneuver. Steady state pitch forces were eliminated in turns with a 
system that automatically inserted the proper amount of pitch command to yield 
the normal acceleration as function of bank angle to hold altitude. 
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The gain and real root of the angle of attack transfer function were 
chosen for most of the configurations to put the instantaneous center of 
rotation at the pilot position, which was 33.8 ft. forward of the center of 
gravity in the TIFS. This was done to eliminate from the evaluations the lead 
that one gets from q acting through the lever arm distance of the center of 
rotation to pilot position. The configurations which had different centers of 
rotation are indicated in the following summary. 

2.2.2 Longitudinal Configurations 

A brief summary of the specific configurations are now presented, 
followed by the complete transfer functions and F and G matrices. Table 2 
lists a summary of the configuration characteristics. 


Configuration B 

This was a baseline conventional airplane configuration that was 
chosen to yield Level I flying qualities, about which time delay and pitch 
sensitivity variations could be made to investigate their effects on flying 
qualities. It was based on a TIFS aerodynamic model with increased ^ and Mq 
derivatives in order to achieve an u> sp = 2 rad/sec and a s sp = .7, lA 62 = .75 
N z /a = 5.3. The instantaneous center of rotation was 22.2 ft aft of the 
pilot. 


Configurations 1-8 

Form basic set to evaluate different command response types. 


1. a cmd (short/long term) 

2. q cmd (short/long term) 

3. a (short term)/q (long term) 

4. q (short term)/a (long term) 

5. a cmd 


***sp = 2, £ = .7, uiph = .3, 1 Aq 2 = .5, 
N z /a = 3.5 

w sp = 2, £ = 2.1, <% = .3, 1 /tq 2 = .5, 
N z /a = 3.5 

^sp = 2, Z = .7, (4 a = .3, 1 /tq 2 = .5, 
N z /a = 3.5 

^sp = 2, t = 2.1, a>ph = *3, 1 /t Q 2 = .5, 
N z /a = 3.5 

u sp = 2, Z - .7, u)ph = .1, 1 /tq 2 = .9, 
N z /a = 6.3 
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Table 2 

CONFIGURATION CHARACTERISTICS 



Note* Configurations with phugoid characteristics 

listed as (-) have real roots where indicated. 

1/tqj =0.1 for Configurations 1-14 
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6 . 


q cmd 


7. 


ct/q cmd 


^Sp = C — l»-^» ^ot = 1/t©2 = 

N z /a = 6.3 

^sp = 2 , £ = .7, Wqj = • 1 j 1 / Tq 2 - « 9 , 
N z /a = 6.3 


8. q/a cmd w S p = 2, 

l/x 0i =0.1 for all configurations. 
Configurations 9, 10 


C = 1*3, Wph - .l f 1/ x 0 2 = 
N z /o = 6.3 


•9, 


The low frequency (phugoid) response is eliminated in both q and a by 
proper zero locations based on Configurations 1, 2. 


9. a cmd 
10. q cmd 


low frequency zeros of q and a at t-lA 01 , 0) = (.1, 0), 

1/t 02 = .5, N z /a = 3.5 

II 


Configurations 11, 12 

Investigate y Cor N z ) command systems referenced to different 
locations, all with o) sp = 2, z, = .7, 1/t 02 = .5, N z /a = 3.5. 

11. y cmd w/r/t CG - CR at pilot 

12. y cmd w/r/t CG - CR at CG 33.8 ft aft of pilot 


Configurations 13, 14 

investigate the use of direct lift control to increase effective l/x 02 . 

13. a cmd hi 1 /t 02 = 2, u) sp = 2, z = .7, N z /a = 14 

14. q cmd hi 1 /t 02 = 2, <u S p = 2, z, = 1, N z /a = 14 

Configurations 15, 16 

Short period frequency variations to gather data for flying qualities 
time history criteria refinement. 

15. ct cmd ***sp = 1» C = •7, Wph = I/X92 = *55, N z /s = 2.5 

(Was not flown and will not be discussed further) 

16. ot cmd ^sp ~ £ ~ ^ph = l/"r 0 2 = 1»0» N z /bi = 7.0 
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Configurations 17-20 


Investigate the. effect of various washout prefilters on a Level 2 
configuration from the previous pitch rate program (Reference 5, Configuration 
1-2-2). Instantaneous center of rotation was 14 ft aft of the pilot. 

These configurations differ from all the others in that they are 
mechanized with a q-feedback path and a proportional plus integral in the 
command path: 


Command 



K L = 4 
Kj = 2 

Washout = (s)/(s + 1 /t w0 ) 


17. 

(1-2-2) 

w sp = 2 

18. 

W.O. 

II 

o 

* 

< 

H 

19. 

W.O. 

1/t W o = 

N> 

O 

• 

W.O. 

> 

s 

II 


8 f 1 /tq 2 = .7, Nz/a = 
.05 
.10 
.20 


5 No washout (W.O.) 


Configurations 21-28 

Time delay and command gain (sensitivity) variations on a baseline 
Level 1 conventional airplane configuration to gather data for criteria 
refinement. 


21. 

Configuration B - 

nominal sensitivity with 

100 ms extra 

delay 

22. 

U 

M 

200 ms 

•• 

23. 

M 

.6(Nominal) sensitivity 

0 ms 

N 

24. 

H 

M 

100 ms 

N 

25. 

M 

M 

200 ms 

N 

26. 

II 

1.5(Nominal) sensitivity 

0 ms 

•1 

27. 

U 

M 

100 ms 

M 





The time delay/sensitivity configuration matrix is shown below* 



Extra time delay (ms) 

Sensitivity 

0 

+100 

+200 

.6(Nominal) 

23 

KB 

25 

Nominal 

B 

■9 

22 

1.5(Nominal) 

26 

27 

28 


The baseline configuration had an effective time delay of 150 ms so 
the increased time delay configurations had effective time delays of 250 and 
350 ms. The three values of sensitivity were defined during checkout phase of 
the program and were chosen to span a range of sensitivities which would still 
yield Level I flying qualities with the minimum time delay value. 

Four additional configurations were added during the evaluation phase 
of the program to gather time delay/sensitivity data on a pitch rate command 
type airplane. The baseline configuration chosen for this set was 
Configuration 17 with a lead/lag filter. . This was Configuration 4-2-2 of the 
previous TIFS/fcitch Rate program. The block diagram of this configuration is 
the same as that shown above for Configuration 17 with a lead/lag filter* 
(1.22 s + l)/( .5 s + 1) replacing the washout block. The nominal command gain 
had to be reduced by a factor of .5/1.22 to keep the initial sensitivity the 
same. 


The time delay/sensitivity configuration matrix based on 
Configuration 17 + LeadAag is* 



Extra time delay (ms) 

Sensitivity 

0 

+100 

+200 

•6( Nominal) 

23A 

24A 

25A 

Nominal 

17 + LA 

21A 

22A 

1.5( Nominal) 

26A 

27A 

28A 


Only Configuration 17 + Lead/Lag, 22A, 25A, and 28A were flown. 
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2.2.2. 1 Transfer Functions 


This section presents a tabulation of the transfer functions of each 
configuration. It is written in the shorthand notation where* 

K(a)[s, to] is equivalent to K(s + a)[s 2 + 2?o)S + uj 2 ] 

The transfer functions are all with respect to force input (F^s). 
The following factors are found in each configuration and are not listed in 
the individual transfer functions: 

Feel system: 20 2 /[.7, 20] 

Gradient: 1/12 in/lb 

Actuator: 20 2 /[.7, 20] 

The command gain, « c (deg/inch), is listed separately for each 
configuration and should multiply the gain of each numerator. The command 
gain was chosen to yield the same maximum q for each configuration for a step 
input (see time histories in Appendix A). A pure transport time delay of 60 
ms should also be added to each transfer function to take into account model 
following lags in the TIFS. This 60 ms delay has been added in the time 
histories shown in Appendix A. 

The velocity numerator to stick command input is the same for all 
configurations (except for B, 17-20, 17 + Lead/Lag, 21-28 and 21A-28A). 

N v = -.05(1)(-15) 

Configuration B Baseline conventional airplane = 2, = .7, 

1/tq 2 = .75, o)ph = .16 

K c = -3.3 

N q = 19.5(0)( .069)( .75) 

N a = 1.(21. 1) [.09, .19] 

N v = -.33 (1 .68) (5.03) (-3. 10) 

D = [.7, 2] [.095, .16] 
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Configuration 1 a-cmd, w sp = 2, ? sp = .7, 1 /t 02 = .5, iOp h = .3 

Kq = -7.8 

N q = -K0X.1X.5) 

N a = -.15(3.7) [.1, .3] 

D = [.7, 2][.l, .3] 

Configuration 2 q-cmd, w sp = 2, ? sp = 2.1, 1 /t 02 = .5 

K c = -11.7 

N q = -1(0)(.1)(.5) 

N a = -. 15(3.7) [.1, .3] 

D = (0)(.1)(.5)(8) 

Configuration 3 a/q, w sp = 2, Csp = .7, 1 /t 02 = .5 

K c = -7.8 

N q = -1(0)(.1)(.5) 

N a = -. 15(3.7) [.1, .3] 

D = (0)(.1)[.7, 2] 

Configuration 4 q/a, a) sp = 2, c sp = 2.1, 1 /t 02 = .5, u> p h = .3 

K c = -11.7 

N q = — 1 (0 ) ( .1)( .5) 

N a = -. 15(3.7) [.1, .3] 

D = ( . 5) (8) [.1, .3] 

Configuration 5 a-cmd, t») sp = 2, t s p = * 7 , 1A© 2 = .9, Wph = .1 

K c = -7.8 

N q = -1(0)(.1)(.9) 

N a = -. 15(3.7) [.1, .1] 

D = [.7, 2] [. 1 , .1] 
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Configuration 6 q-cmd, w sp = 2, e sp = 1.3, 1/t 02 = .9 

K c = -9.1 

N q = -1(0)(.1)(.9) 

N a = -.15(3.7) [.1, .1] 

D = (0)(.1)(.9)(4.4) 

Configuration 7 a/q, w sp = 2, £ sp = .7, 1/t 02 = .9 

K c = -7.8 

N q = -1(0)(.1)( .9) 

N a = -.15(3.7) [. 1, .1] 

D = (0 ) ( . 1 ) [.7, 2] 


Configuration 8 q/a, w sp = 2, c = 1.3, 1/t 02 = .9, w ph = .1 

K c = -9.1 

N q = — 1(0)(.1)( .9) 

N a = -.15(3.7) [.1, .1] 

D = (.9)(4.4)[.l, .1] 

Configuration 9 o, q decoupled from phugoid a-cmd, a) sp = 2, l/x 02 = 

K c = -7.8 

Nq = — 1(0) ( . 1)( .5) 

N a = -. 15(0 ) ( .1)(3.7) 

D = (0)( .1) [.7, 2] 


Configuration 10 a, q decoupled from phugoid, q-cmd, u> sp = 2, l/x 02 

Kq = -11.7 

Nq = — 1(0)(.1)( .5) 

N a = — . 15(0 ) ( . 1)(3.7) 


D = (0)(.1)(.5)(8) 



Configuration 11 


.5 


Configuration 12 


Configuration 13 


Configuration 14 


Y-cmd pilot at CR, 34* fwd of CG, u> sp = 2, 1 /t 02 = 
K c = -40 

N q = -1(0)(.1)(.5) 

N a = -.15(3.7) [.1, .3] 

N nzcG = -• 018(0 )(0)( .955) (-3.86) 

D = (0)(0)(.955)(3.86)(2) 


Y-cmd pilot 34 • fwd of CR, CR at CG, u> sp = 2, 1/t Q2 = .5 

K e = -25.48 

Nq = — 1.(0)(. 1)( .5) 

N a = -.555 [.1, .3] 

N nZcG = -• 055(0 )(0)( 1.274) 

D = (0)(0)(1.274)[.7, 2] 

High 1 /tq 2 = 2, a-cmd, a) sp = 2, ^ sp = .7 

K c = -7.8 

Nq = —1(G) ( . 1 ) (2 ) 

N a = -.15(3.7) [.1, .3] 

D = [.7, 2] [. 1 , .3] 

High 1 /tq 2 = 2, q-cmd, a)s P = 2, £sp = 1 

K c = -7.8 

N q = —1(0) ( . 1) (2) 

N a = -.15(3.7) [.1, .3] 

D = (0)(.1)(2)(2) 
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Configuration 15 


Not flown 


Configuration 16 a-cmd, a ) sp = 3, £sp = .7, 1 /t 02 = 1, u> p h = .3 

K c = -7.8 

Ng = -1(0)(.1)(1) 

N a = 15(3.7) [.1, .3] 

D = [.7, 3] [. 1, .3] 

Configuration 17 Configuration 1-2-2 of previous TIFS/^itch Rate 

program with no washout 

K c = -1.7 

N q = -1417(0 )( .056)( .72) (2) 

N a = -129. 5(11. 4)(2)[.106, .2] 

N v = -12. 3( .878) (2) (-54) 

D = (0 ) ( . 052) ( .82) [.78, 2.9] [.7, 17.1] 

NOTE: [.7, 17.1] is the actuator pole which has migrated from [.7, 20] 

as it is now in the loop. 


Configuration 17 + |_ead/Lag Configuration 4-2-2 of previous TIFS/^itch Rate 

program where lead/lag is: 

(1.22 s + l)/(.5 s + 1) = 2.44 (s + .82)/(s + 2) 


Kc 

Nq, 


D 


= -.74 


Na. 


N v - same as Configuration 17 
with 2.44 (.82) factors 

> 

same as Configuration 17 
with (2.) factor j 


2.44 (s + .82 ) 
(s + 2) 


l 
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Configuration 18 Configuration 17 with 1 /t wo = .05 
K c = -1.7 . 

N q> N a , N v - same as Configuration 17 with (0) factor 
D = same as Configuration 17 with (.05) factor 

Configuration 19 Configuration 17 with 1 /t w0 = .1 

K c = -1.7 

N q> N a , N v - same as Configuration 17 with (0) factor 
D = same as Configuration 17 with (.1) factor 

Configuration 20 Configuration 17 with 1 /t w0 = .2 
K c = -1.7 

N q» N a , N v - same as Configuration 17 with (0) factor 
D = same as Configuration 17 with (.2) factor 

Configurations 21-28 Time Delay/Sensitivity variations based on a 

conventional airplane 




) — 

J s + .05 


These configurations are identical to the previously described 
Configurations B except for the command gain, K c and extra time delay above 
the baseline 150 ms. 


Config 21 
Con fig 22 
Config 23 
Config 24 
Config 25 
Config 26 
Config 27 
Config 28 


Config B with 
Config B with 
Config B with 
Config B with 
Config B with 
Config B with 
Config B with 
Config B with 


K c = -3.3, 
K c = -3.3, 
K c = -2.0, 
K c = -2.0, 
K c = -2.0, 
K c = -5.0, 
K c = -5.0, 
K c = -5.0, 


+100 ms delay 
+200 ms delay 
+0 ms delay 
+100 ms delay 
+200 ms delay 
+ 0 ms 

+100 ms delay 
+200 ms delay 
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Configurations 21A-28A Time Delay, Sensitivity variations based on the pitch 

rate command system. 


These configurations are identical to the previously described 
Configuration 17 + Lead/Lag except for the command gain, Kc, and extra time 
delay above the baseline 150 ms. 


Config 21A 
Con fig 22A 
Config 23A 
Config 24A 
Config 25A 
Config 26A 
Config 27A 
Config 28A 


Config 17 + 
Config 17 + 
Config 17 + 
Config 17 + 
Config 17 + 
Config 17 + 
Config 17 + 
Config 17 + 


L/L with K c 
l_/L with K c 
L/l_ with K c 
L/_ with K c 
l_/L with K c 
L/L with K c 
l/L with K c 
L/t with K c 


-.74, +100 ms 
-.74, +200 ms 
-.44, +0 ms 
-.44, +100 ms 
-.44, +200 ms 
-1.11, +0 ms 
-1.11, +100 ms 
-1.11, +200 ms 


Only Configurations 22A, 25A, and 28A were flown. 


2. 2. 2. 2 Matrix-Vector Representation 

This section presents the linearized model equations of motion in the 

form of: 


X = Fx + Gu 


where the state vector X = 


q, deg/sec 
e, deg 
a, deg 
V, ft/sec 


the control vector u = [5 , deg] 

e 


In addition, each configuration has a second order feel system 
U = .7, u) = 20 rad/sec) and a second order actuator (c = .7, w = 20 rad/sec) 
and a command gain K c (deg/inch). The feel system gradient is 12 lb/in for 
all configurations. There is also an additional 60 ms of pure transport time 
delay to take into account model following lags in the TIFS. 
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The F and G matrices for each of the configurations are presented; 


• 


■ 

q 


q 

• 



0 

= F 

e 

a 


a 

• 



V 


V 


Configuration B Baseline conventional airplane, w sp = 2, 1 /t 02 = .75 



-1.9000 

.0000 

-2.2900 

.0246 



1.9500 

F = 

1.0000 

.0000 

.0000 

.0000 


G = 

.0000 


1.0000 

.0000 

-.9010 

-.0711 



-.1010 


.0000 

-.5614 

.2618 

-.0382 



-.0326 

Configuration 1 

a-cmd 

CM 

II 

Ql 

if) 

3 


1/t 02 = .! 


-1.4047 

-.1536 

-5.7022 

.0896 " 


-1.0000 

F = 

1.0000 

.0000 

.0000 

.0000 

G = 

.0000 


.3528 

.0165 

-1.4520 

-.0237 


-.1500 


-.1440 

-1.3427 

1.2930 

-.0033 


.0000 

Configuration 2 

g-cmd 

w sp = 2 


1/t 62 = .5 


’-8.0000 

.0000 

.0000 

.0000 ' 


’-1.0000 

F = 

1.0000 

.0000 

.0000 

.0000 

G = 

.0000 


-.6365 

.0396 

-.5967 

-.0371 


-.1500 


-.1440 

-1.3427 

1.2930 

-.0033 


.0000 

Configuration 3 

a/g 

CM 

II 

Ql 

CO 

3 


1/t 02 = .5 


’-1.5868 

-.5784 

-4.7550 

-.3556 ’ 


’-l.oooo' 

p _ 

1.0000 

.0000 

.0000 

.0000 

G = 

.0000 


.3255 

-.0472 

-1.3099 

-.0905 


-.1500 


-.1440 

-1.3427 

1.2930 

-.0033 


.0000 


K c = -3.3 


K c = -7.8 


K c = -11.7 


K c = -7.8 
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Configuration 4 


q/a 


w sp = 2 


1/tq 2 = .5 



’-8.1720 

-.4709 

1.4131 

.5428 ' 


-1.0000 ' 

F = 

1.0000 

.0000 

.0000 

.0000 

G = 

.0000 


-.6623 

-.0311 

-.3847 

.0443 


-.1500 


-.1440 

-1,3427 

1.2930 

-.0033 


.0000 


K r = -11.7 


Configuration 5 


a-cmd 


u sp 


= 2 


F = 


1.0000 

-.1166 

.0216 


.0000 

.0400 

.8224 


.0000 

-.9027 

.1894 


.0000 

-.0115 

-.0973 


1/Tq 2 = .9 


F = 

-1.0030 

1.0000 

.3930 

.0216 

.0652 

.0000 

-.0302 

-.8224 

-5.4467 

.0000 

-1.7197 

.1894 

.0205 " 
.0000 
-.0085 
-.0973 

G = 

-1.0000 " 
.0000 
-.1500 
.0000 

Configuration 6 

q-cmd 


w sp = 2 

1/t 02 = .9 


r -4 . 4000 

.0000 

.0000 

.0000 ] 


r - 1.0000 1 


G = 


.0000 

-.1500 

.0000 


K r = -7.8 


K r = -9.1 


Configuration 7 


a/q 


^sp = 2 


1Aq 2 = .9 



’-1.1634 

-.2828 

-4.9108 

-.0706 " 


’-1.0000 ' 

F = 

1.0000 

.0000 

.0000 

.0000 

G = 

.0000 


.3689 

-.0824 

-1.6393 

-.0221 


-.1500 


.0216 

-.8224 

.1894 

-.0973 


.0000 


K r = -7.8 


Configuration 8 


q/a 


^sp = 2 


1/tq 2 = .9 



’-4.3232 

.3206 

.0216 

.0914 ’ 


’-1.0000 ' 

F = 

1.0000 

.0000 

.0000 

.0000 

G = 

.0000 


-.1051 

.0081 

-.8995 

.0022 


-.1500 


L .0216 

-.8224 

.1894 

-.0973 


[ .0000 


K r = -11.7 
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Configuration 9 


a-cmd uncoupled from phugoid 


-1.4094 

1.0000 

.3436 

-.1922 


.0000 

.0000 

.0000 

-1.5000 


-5.9376 

.0000 

-1.3906 

1.6146 


- 1.0000 

.0000 

-.1500 

.0000 


Configuration 10 q-cmd uncoupled from phugoid 


-8.0000 

.0000 

.0000 

.0000 


-1.0000 

1.0000 

.0000 

.0000 

.0000 

G = 

.0000 

-.6450 

.0000 

-.5000 

.0000 


-.1500 

-.1922 

-1.5000 

1.6146 

-.1000 


.0000 


Configuration 11 


Y-cmd CR at Pilot 34* fwd of CG 


-3.7598 .3405 -3.0348 -.1796 

1.0000 .0000 .0000 .0000 

-.0005 .0907 -1.0519 -.0641 

[ -.1440 -1.3427 1.2930 -.0033 

also requires 2/(s + 2) prefilter 


- 1.0000 

.0000 

-.1500 

.0000 


Configuration 12 


Y-cmd CR at CG 34* aft of pilot 


-2.2000 

.0815 

-4.8569 

-.3184 " 


- 1.0000 ' 

1.0000 

.0000 

.0000 

.0000 

G = 

.0000 

.5550 

.0261 

-.5837 

-.0372 


.0000 

.0500 - 

■1.3336 

1.0891 

-.0163 

1 

.0000 

also 

requires 

1.274/(s 

+ 1.274) 

prefilter 
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Configuration 13 

a-cmd 

high 1 /tq 2 

= 2 





.6294 

-.1467 

-9.2630 

-.1766 i 



’-1.0000 ' 

F = 

1.0000 

.0000 

.0000 

.0000 


G = 

.0000 


.6266 

-.0514 

-3.2772 

-. 1602 



-.1500 


.2751 

-.4206 

-1.5010 

-.2122 



.0000 

Configuration 14 

q-cmd 

high 1 /tq 2 

= 2 





’-2.0000 

.0000 

.0000 

.0000 " 



’-1.0000 ' 

F = 

1.0000 

.0000 

.0000 

.0000 


G = 

.0000 


.2322 

-.0294 

-1.8878 

-.1337 



-.1500 


.2761 

-.4206 

-1.5010 

-.2122 



.0000 

Configuration 15 

a-cmd 

w sp = 1 

II 

CM 

(D 

i-H 

.35 




’ -.8546 

-.1724 

-1.0357 

.0429 




"- 1.0000 " 

F = 

1.0000 

.0000 

.0000 

.0000 


G = 

.0000 


.4445 

.0461 

-.6635 

-.0241 



-.1500 


-.2672 

-1.7754 

2.1148 

.0581 




.0000 

Configuration 16 

a-cmd 

w sp = 3 

1A 02 = 

1.0 




’-1.0121 

-.0585 

-14.3193 

.1176 




’-1.0000 ' 

F = 

1.0000 

.0000 

.0000 

.0000 


G = 

.0000 


.3972 

-.0135 

-3.1479 

-.0479 



-.1500 


.0500 

-.7500 

.0000 

- . 1000 



.0000 


K r = -7.8 


K r = -7.8 


K r = -7.8 


-7.8 


Configurations 17-20 Based on previous TIFS/^itch Rate Configuration 

1-2-2* It has loop gain K|_ = 4, integral gain 
Kx = 2 (see block diagram in previous section) 
Unaugmented F and G* 



’ -.4880 

.0000 

.5200 

.0053 ' 


’ -.8860 ' 

F = 

1.0000 

.0000 

.0000 

.0000 

G = 

.0000 


1.0000 

.0000 

-.6870 

-.0740 


-.0810 


.0000 

-.5610 

.0950 

-.0460 


L .0000 
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Configurations 20-28 Time delay/sensitivity variations based on conventional 

Configuration B. F and G matrices are the same as those 
for Configuration B. 

Configuration 21 = Configuration B with K c = -3.3, +100 ms 

Configuration 22 = Configuration B with K c = -3.3, +200 ms 

Configuration 23 = Configuration B with K c = -2.0, +0 ms 
Configuration 24 = Configuration B iwth K c = -2.0, +100 ms 

Configuration 25 = Configuration B with K c = -2.0, +200 ms 

Configuration 26 = Configuration B with K c = -5.0 +0 ms 

Configuration 27 = Configuration B with K c = -5.0 +100 ms 
Configuration 28 = Configuration B with K c = -5.0 +200 ms 
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Configurations 21A-28A Time delay/Sensitivity variations based on pitch rate 

command Configuration 17 + Lead/Lag. F and G matri- 
ces are the same as those for Configuration 17 + 
Lead/Lag. 

Configuration 21A = Configuration 17 + L/L with K c = -.74, +100 ms 

Configuration 22A = Configuration 17 + L/L with K c = -.74, +200 ms 

Configuration 23A = Configuration 17 + L/L with K c = -.44, +0 ms 

Configuration 24A = Configuration 17 + L/L with K c = -.44, +100 ms 

Configuration 25A = Configuration 17 + L/L with K c = -.44, +200 ms 

Configuration 26A = Configuration 17 + L/l with K c = -1.11, +0 ms 

Configuration 27A = Configuration 17 + L/L with K c = -1.11, +100 ms 

Configuration 28A = Configuration 17 + L/L with K c = -1.11, +200 ms 

Time histories for a 10 pound F^s step (5 sec in, 5 sec out) are 
shown in Appendix A. Also presented are the first one second of the time 
histories on an expanded scale, so the effective time delay in the responses 
can be seen. The effective time delay in the pitch rate response (maximum 
slope intercept method) is shown in Table 3. 
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Table 3 

MEASURED EFFECTIVE TIME DELAY 
(PITCH RATE MAXIMUM SLOPE INTERCEPT) 

(Includes All Systems Delay and Model Following Lags, 60 ms) 


Configuration 


Time Delay (ms) 


1 

160 

2 

150 

3 

160 

4 

140 

5 

160 

6 

150 

7 

160 

8 

150 

9 

160 

10 

150 

11 

250 

12 

290 

13 

160 

14 

160 

16 

160 

17 

150 

17 + Leao/Lag, 23A, 26A 

150 

18 

150 

19 

150 

20 

150 

B, 23, 26 

150 

21, 24, 27 

250 

22, 25, 28 

350 

21A , 24A, 27A 

250 

22A, 25A, 28A 

350 
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2.2.3 Lateral/tJlrectlonal Description 


The lateral/directional aerodynamics and control system were the same 
as those used for the previous TIFS/fcitch Rate program. The aircraft model 
had the following characteristics. 


Constant Characteristics 


Weight (W) = 

Mass (M) = 

Wing Area (S) = 

Wing Span (b) = 

Wing Chord (c) = 
Ixx = 

% 

Izz 

Ixy = 

Trim Conditions 
v Trim - 

q 


193,000 lbs. 

5999.4 slugs 
2147 ft 2 
157 ft 
15.074 ft 

4,003,900 slug-ft 3 
5,408,550 slug-ft 3 
9,184,470 slug-ft 3 
223,410 slug-ft 3 


132 K1AS (223 fps) 
59.14 psf 


Lateral/t)irectional Non-Dimensional Derivatives (per degree) 


c y& 

-0.03136 

C VP 

0.00563 

C Yr 

0.01345 

c y6 r 

0.00536 


-0.00256 

C *P 

-0.01022 

Ct r 

0.00749 

CL 6 a 

0.00148 


0.00023 
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Lateral/ftirectional Non-Dimensional Derivatives [per degree) (Cont'd) 


Ci 6 r 

0.00050 

c n e 

0.00394 

Cn P 

-0.00074 

Cn r 

-0.00552 

Cn «a 

0.00023 

Cn,s SP 

0.00024 

Cn 6 r 

-0.00169 


The lateral/directional flight control system is shown in Figures 5 
and 6. The gains were adjusted to achieve the following characteristics, 
which were found to be Level 1 and consequently "transparent'* to the longitu- 
dinal investigation. 


Dutch Roll Mode 


ion (rad/sec) 

0.768 

K 

0.297 

tojp (rad/sec) 

0.228 

P (sec) 

8.57 

4>/6 

0.188 


Roll Spiral Mode 


(sec) 

- 

x s (sec) 

- 

w rs (rad/sec) 

4.741 

Crs 

0.369 

C rs wrs ( rad/sec) 

1.748 

P rs (sec) 

1.47 


Effective roll mode time constant 
(time to 63 % max roll rate) 1.0 sec 
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Figure 5 ROLL CONTROL SYSTEM 




2 ™ 







Roll Control Parameters 


aty/wd 1 . 000 

ty/5d 1.001 

Time histories for a 20 degree 6/w step and a 1/2 inch 6pp step are 
shown in Figures 7 and 8. 

2.2.4 Thrust 

The throttle control system was the same for all of the configura- 
tions, though the response of the airplanes differed because of the different 
dynamics in the rest of the longitudinal axes. The pilot controlled the 
thrust through a throttle lever on the center console. It had a ±7056 travel 
about its mid-trim position, it acted through a first order lag with a time 
constant of .2 sec to move a power lever angle (PLA) command which produced a 
V of ±8.4 ft/sec 2 at ±7056. 


th 



Time histories for 1056 6th steps are shown in Appendix B for each config- 
uration. 


2.2.5 Feel System and Actuator Dynamics 

The feel system parameters were chosen during the calibration 
flights. Known good feel system parameters from the previous TIFS/Pitch Rate 
program were used as the starting point and were slightly modified to provide 
better characteristics. Idealy the feel system should have been “transparent** 
to the flight test. The lack of evaluation pilot comments concerning the feel 
system indicates that this was achieved. 

Below are the model feel system parameters that resulted from the 
calibration flights and were used throughout the investigation t 
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Column 

Pitch 

■H 

Pedals 

Yaw 

Gradient 

12 (lbs/in) 

0.37 (lb/deg) 

36 (lb/in) 

Wn (rad/sec) 

20* 

25 

12 

C 

0.7 

0.7 

0.6 

Breakout (lbs) 

4.0 

1.0 

12 


* (25 rad/sec was used in the TIFS to help reduce 

the effective time delay to that of the model). 


The actuator dynamics used in the model for the elevator, aileron, 
and rudders were second order with uip = 20 rad/sec, 5 = 0.7. As previously 
mentioned the throttle had a first order lag with t = .2 sec. 


2.2.6 Turbulence Sensitivity 


External disturbances were imposed upon the evaluation in the form of 
a discrete vertical gust (1-cosine with a maximum amplitude of 7.6 ft/sec or 
ogust = 1.9 deg with a duration of 4 sec). To make all of the configurations 
respond to this discrete gust in a consistant manner, a constant turbulence 
sensitivity vector was chosen for the oguST* The values chosen were those for 
the NASA TCV/737 with neutral static margin (Reference 15). The following 
turbulence sensitivity terms were added to the state equations? 



Though the turbulence sensitivities were the same for all configurations, the 
responses differed because of different free responses of each configuration. 
Any natural turbulence, if present, was not fed to the model, but did disturb 
the model-following system primarily in the heave axis where it was felt as a 
light chop and was seen in the airspeed indicator which used the TIFS 
airspeed. Most of the flights had light or no turbulence, so it was not a 
factor. 
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2.2.7 Ground Effect 


The ground effect chosen for the model was similar to that used on 
the previous TIFSAitch Rate program, though no X-force due to ground effect 
was used. The following ground effect sensitivity terms were added to the 
state equations; 


• 


■ 

q 


-.08 

• 



e 

• 


0 

a 

• 


-.25 

_ V _ 


0 


• f (h) 


where f(h) was the following ground effect function j 


h 

f(h) 

>100 

0 

50 

.1 

40 

.14 

30 

.2 

20 

.32 

10 

.5 

SO 

1 



2.2.8 Bank Compensation 

The steady state forces to hold a given normal acceleration varied 
greatly from one configuration to another. To relieve pilot pitch forces in 
turns and remove this factor from the evaluation which was primarily a longi- 
tudinal task, an automatic altitude-hold-in-turns system was used. For the 
configurations which had actual pitch rate feedback and an integrator 
(Configurations 17, 17 + Lead/Lag, 18-20, 21A-28A), this system was essen- 
tially a (-gA) $> 2 feedback to the pitch command summer. This had been used 
in the previous TIFSAitch Rate program. For the balance of the configu- 
rations, which did not have an integrator in the loop, the elevator command to 
hold the proper normal acceleration for a level turn as a function of bank 
angle had to be determined. The following derivation shows how this elevator 
required in turns or bank compensation was calculated. 
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BANK COMPENSATION 


(1) a = Z q + Z 0 e + Z a a + Z u u + Z fi 6 + (cos $-1) (57.3) 

e 

(2) q = M q + M Q o + M a a + M u u + Mg 6 g 

e 

in steady turn a = 0, q = 0 

q = (57.3) ^ (n - jlj-) n = (load factor) ^ 


From Equation (1) 


(3) 


a 


-Z g q -Z e e - z u u - Z 6 J e - (cos 4>-i) $ (57.3) 

S 


From Equation (2) 


(A) 


S e = 




V - M a“ - M u u 


M, 


Substituting q 


9 rn ^ 
= V (n " n 5 


57.3 


-Mq «7 [ n ' *1 > 57 ’ 3 - V ' V - M u u 

\ 

Inserting Equation (3) into Equation (5) yields 
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Ve = - M q ( v 57 - 3 - V - V - M o 


(- z q q - Z e e - Z u u - Z 6 6 e -(cos *-1) £ 57.3) 


Gather 6 terms on left side of equation 




M_ 


v« - 1 v« ■ l \ - i v *• 


s ,‘ 


1 ( V [ n " ol )57 - 3 - M e® - M u u " T c " z q ( V [ n - al )57 * 3 ' z e ® - z u u ” (cos *- 1) V 57 * 3 ^ 


M 6 e ' F Z 6 

e a e 


{ e = 


M 

-M + t 2 Z_ 
q Z a q 


7 <" " b 


M H M 

37,3 " M e e ” V + z^ Z e® + z^ Z u u + ^ cos V 37,3 
a a a 


M 

Mfi ■ £ \ 
a e 


Make following assumptions: Ae in turns is very small 

Au in turns is very small 

so drop out Mge, M u u, z e e, Z u u 


M 

+ T a Z * 




M a 9 

57.3 + 2 = (COS 4-1) ^ 57.3 
a 


M 

M 6 * T a Z 6 

e a e 


{ e = 


-57.3 7 M q (n - i ^ [z q (n - £) + 1 (cos 4* - 1)] 57.3 


M 

M 6 - T 2 6 
®e Z o ®e 
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-57.3 i 


* e = Z a\ “ M o Z 6 e IW" " ^ ’ M o [ Z q (n - £) + 1 (cos * - 1)]] 


fi e = 


-57.3 S 


= Z a M 6 - M o Z 6 " M °V (n " " 3 " M a (cos ♦ " «] 

e e J 


n-±. 1 


n cos 


x - cos 4 = 1 ~ c °s 2 <b _ (1 - cos 4X1 + cos 6 } 
4 T cos 4 cos 4 


g 

-57.3 i 


* V ,1 - COS 4, r , 

e = Z a M 6 - Z 6 M a ( cos 4 3 [ CZ a M q “ M a Z q 3(1 + cos + M a cos *] 


e e 


5 e ■ - 57 - 3 V ti -isif i) 


- h\ (1 + C ° S W * W°- z« M a (C0S * J 

e e e e 


S e = -V 1 cos°| ♦) ♦ K 2 (cos ♦)] 


K o * ^r- 2 * e - 14 


Z M - M Z 

X 2L_9 a Q 

1 Z M - - Z. M 

a fi e « e a 


M 


Ko = 


Z a\ - \* a 

e e 
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Using the values for Ma, M q , M$ e , Za, Z q , Z$ e from the F and G matri- 
ces, the bank compensation gains and K2 were calculated and are listed in 
Table 4 . 

Table 4 

BANK COMPENSATION GAINS 

f Configuration Kl «2 

B 2.62 -1.5 

1 6.79 -9.56 

2 8.0 0 

3 6.08 -7.97 

4 6.84 2.37 

5 4.28 -6.03 

6 4.4 0 

7 4.12 -5.44 

8 4.31 .02 

9 8.0 -11.88 

10 8.0 0 

11 6.63 -5.09 

12 6.82 -8.32 

13 1.98 -4.91 

14 2.0 0 

16 8.87 -14.32 

21-28 2.62 -1.5 



2.3 


EVALUATION TASKS AND PROCEDURES 


The evaluation task and procedures were the same as those used in the 
previous TIFS/Pitch Rate program, and are described below. 

The evaluation pilot was given control of the aircraft on the down 
wind leg and performed a visual turning approach to a 1.5 to 2 mile final 
approach. The ILS glide slope was intercepted in the turn and held to a point 
3500 ft. from the runway/glide slope intercept point. A constant speed of 132 
KIAS was held thoughout the approach until landing flare. 

Figure 9 details the final approach and flare geometry. A final 
approach "barrier" was defined as projecting up from the ground at a point 
3500 ft. short of the runway and glide slope intercept point and up to the ILS 
glide path. The evaluation pilot was not allowed to descend below the ILS 
glide slope until passing the "barrier" (the position 3500 ft. short of the 
runway/glide slope intercept is well marked by a railroad track). Peer 
pressure from the safety pilots and the flight test engineer was found to be 
quite sufficient to prevent barrier duck under. 

In addition to the altitude constraint of the barrier, lateral off- 
sets of 200 ft. (either left or right, and obtained visually by lining up on 
runway markers) were used to provide secondary tasking thus preventing pre- 
occupation with the pitch task, in order to further assure pitch task activ- 
ity a (1 - cosine) angle of attack gust was fed to the model in the zone 
depicted in Figure 9 (between 100 and 50 feet of altitude). 

The "Desired" touchdown area was defined as being 500 ft. long and 20 
ft. wide (+10 ft. of centerline) starting 250 ft. past the runway/glide slope 
intercept. The "Adequate" touchdown area was defined as 1000 ft. long, 40 ft. 
wide and starting at the same point on the runway. Airspeed requirements 
were: "Desired** 132 +3 KIAS, "Adequate" 132 ±5 KIAS, both at barrier passage. 

"Desired" sink rate at touchdown was defined as 0 to 3 fps and "Adequate" as 3 
to 6 fps (these values were obtained from the data records, however, experi- 
ence has shown that 0-3 fps touchdowns result in "smooth" landings, 3-6 fps 
touchdowns result in "solid" landings, ana touchdowns in excess of 6 fps can 
be recognized by any crew member with a 95% confidence level). 
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Figure 9 APPROACH AND LANDING TASK 
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Touchdown parameters were obtained from on-board recordings using DME 
and Localizer Deviation for x and y distance on the runway, and as described 
above for sink rate. Speeds were taken by the test engineer in the right 
evaluation seat and also recorded. These combined with the evaluation pilots 
comments and ratings provided the evaluation data. 

The design goal of the above task was to achieve sufficient pilot 
gain in the pitch axis to provide an adequate spread in the Handling Qualities 
Ratings (HQR’s) but not to be so difficult or easy so as to bias the HQR’s one 
way or another. 

2.3.1 Evaluation Sequence 

The evaluation pilot was given control of the aircraft on the down- 
wind leg. At this time automatic step inputs were made in order to verify the 
configuration. The evaluation pilot would then conduct a visual turning 
approach to intercept the glide slope approximately two miles from touchdown. 
He would visually line up for the offset and continue down the glide slope. 
At the "barrier” position (3500 ft. from touchdown and at approximate altitude 
of 200 ft.) the correction for the offset would commence and thrust and pitch 
would be adjusted for landing. At approximately 2000 ft. from touchdown a 
discrete (1 - cosine) angle of attack gust would be fed to the aircraft model 
to cause a standard flight path disturbance. The evaluation pilot would fly 
through the gust to flare and touchdown. At touchdown the safety pilots would 
take control of the aircraft (or at any time prior to touchdown if dictated by 
the situation). At this point the flight test engineer would record speeds 
and estimated touchdown dispersion and the evaluation pilot would begin his 
comments and give the HQR. This data was manually recorded by the flight test 
engineer as well as on voice tape. The safety pilots then executed the climb- 
out while the TIPS technical crew set up the next configuration (if required) 
to repeat the process. 

A normal evaluation consisted of two approaches, however, the eval- 
uation pilot had the option of repeats if desired. 
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2.4 


PILOT COMMENT CARD AND RATING SCALES 


The evaluation pilots were briefed on the general experiment purpose 
and flight task details. They had a general knowledge of what the test con- 
figurations were (had seen descriptions and time histories) but no knowledge 
of which of those configurations would be given on each flight. 

Pilot technique was necessarily different for pitch rate command 
flight control systems than for " conventional” systems (a-command) that 
require increasing average pull forces in the landing flare (i.e., monotonic). 
Consideration was given to informing the evaluation pilots, before hand, which 
type of system they had. It was decided to proceed with a “blind” experiment 
as it had been found successful in the previous TIFS/Pitch Rate program. The 
pilots adapted to technique changes rapidly and in some cases were unaware of 
using different techniques. 

An evaluation normally consisted of two approaches and landings. The 
pilot could make comments at any time, however, formal use of the comment card 
(Figure 10), Cooper-Harper scale (Figure 11), and the PIO scale (Figure 12) 
was made after the second landing for the configuration. The pilot had the 
option of a third landing on a configuration and in that case the comments and 
gradings were made after the third landing. 

The pilot comments and ratings were considered the primary data of 
the investigation, and were recorded on voice tape. In addition, the flight 
test engineer (in the right evaluation seat) manually recorded comment sum- 
maries, touchdown dispersion, and pilot ratings for use in the post-flight 
debriefing where pilot comments were elaborated in more detail. 
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tifs/Flared landing program 

PILOT COMMENT CARD 


Initial Overall Impression 
Approach 

1. Initial/Final response to control inputs 

2. Flight path control 

3. Pitch attitude control 

4. Airspeed control 

5. Offset correction 

6. Atmospheric disturbances 

7. Special pilot techniques 

Flare and Touchdown 

1. Pitch attitude and flight path control 

2. Control of touchdown parameters 

3. Atmospheric disturbances 

4. Special pilot techniques 

Pilot Ratings 

1. Approaches 

2. Flare and touchdown 

3. Overall 

4. PIO rating 


* E 


/igure 10 PILOT CO*CNT CARD 



HANDLING QUALITIES RATING SCALE 


ADEQUACY FOR SELECTED TAW OR 


^ It H > 

satisfactory without 
v improvement? > 


| Deficiencies 
► warrant 
improvement 


^isadequate\ 

/ performance \ No | Deficiencies 
Attainable with a toierabi ^ ^ require 
pilot workload 7 I improvement 


AMCRAFT 

CHARACTERISTICS 

DEMANDS ON THE PILOT 
m SELECTED TASK OR REQUIRED OPERATION* 

Excellent 

Pilot compensation not a factor lor 

highly desirable 

desired performance 

Good 

Pilot compensation not a factor for 

Negligible deficiencies 

desired performance 

Fair — Some mildly 

Minimal pilot compensation required for 

unpleasant deficiencies 

desired performance 


Minor but annoying 

Oetired performance requires moderate 

deficiencies 

pilot compensation 

Moderately objectionable 

Adequate performance requires 

deficiencies 

considerable pilot compensation 

Very objectionable but 


tolerable deficiencies 

pilot compensation 


Major deficiencies 

Adequate performance not attainable with 
maximum tolerable pilot compensation. 
Controllability not in question 

Major deficiencies 

Considerable pilot compensation is required 
for control 

Major deficiencies 

Intense pilot compensation is required to 
retain control 


it 

it controllable 7 


^ Improvement 
mandatory 


Major deficiencies 


Control will be lost during some portion of 
required operation 


► DetmiMy’ o* 'tqui'M OOerat*** d*»<fln*t*on Of tlfgM prvji* ®no c 

lutXKoMt aritn accompanying conOil’On* 


Figure li COOPER-HARPER HANDLING QUALILTIES RATING SCALE 


ORIGINAL PAGE IS 
Of POOR QUALITY 
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4 


5 


Figure 12 PIO TENDENCY CLASSIFICATION 
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2.5 


EVALUATION PILOTS AND FLIGHT SUMWWY 


This flight program was extremely fortunate to have the services of 
eight evaluation pilots with a wide variety of backgrounds. In addition to 
the NASA and Calspan pilots who were scheduled to participate in this NASA 
sponsored program, NASA invited other flight test organizations to provide 
pilots if they would pick up the costs for their flights. Boeing, Lockheed, 
and the German Aerospace Research Establishment (DFVLR) accepted this invita- 
tion and although their flights were conducted under a separate contract the 
results of their flights are presented in this report. 

The evaluation pilots werei 

Charles 3. Berthe, Jr. - 

Calspan checkout pilot (also co-investigator and safety pilot during 
evaluation phase) flew the calibration flights which determined the 
accuracy of the simulation, the efficacy of the task, fine tuned the 
command gains, time delays, and lateral/directional axes for the 
evaluation phase. 

Lee H. Person, Jr. - 

NASA/Langley evaluation pilot (A) - flying experience with aircraft 
handling qualities and IFR display programs including NASA’s Terminal 
Configured Vehicle - 737. 

John F. Ball - 

Calspan evaluation pilot (B) - flying experience as an evaluation 
pilot in many handling qualities programs, also a safety pilot in the 
USAF/TIFS, NT-33, and Calspan Learjet in-flight simulators. 

Roger E. Smith - 

NASA/tDryden evaluation pilot (C) - flying experience as an evaluation 
pilot in many handling qualities programs at Calspan and NASA, also 
test pilot on X-29 and AFTI/F-111. 
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Dale M. Ranz - 

Boeing evaluation pilot (D). 

J. Kenneth Higgins - 

Boeing evaluation pilot (E) 

- Production and research test pilots for all Boeing commercial 
transports including evaluation of advanced flight control concepts. 

Frank Hadden - 

Lockheed evaluation pilot (F) _ Chief Engineering Test Pilot at 
Lockheed - Georgia, with extensive flight test expeirence with 
transport category aircraft including the Jet Star C-130, C-141, C-5, 
and High Technology Test Bed (C-130). 

Hans Meyer - 

DFVLR (Germany) evaluation pilot (G) - Chief Test pilot at DFVLR with 
extensive flight test experience, including participation in in- 
flight simulation programs with the DFVLR-Hansa Jet. 

Most of the evaluation approaches were conducted at Niagara Falls 
AFB, with approximately one-fourth of the evaluations conducted at Buffalo due 
to traffic and weather constraints at Niagara. There was no discernible dif- 
ferences in the evaluations at the two airports, but x-runway distance was not 
as well marked at Buffalo. 

A summary of the flights versus pilot is shown in Table 5. 
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Table 5 

FUGHT SUMMARY 















r , Section 3 

***** * EXPERIMENT MECHANIZATION 

3.1 EQUIPMENT 

The USAF Total In-Flight Simulator (TIFS was used as the test vehicle 
in this experiment. TIFS is a highly modified C-131 (Convair 580) configured 
as a six-degree-of-freedom simulator (Figure 13). It has a separate evalu- 
ation cockpit forward and below the normal C-131 cockpit. When flown from the 
evaluation cockpit in the simulation or fly-by-wire mode, the pilot control 
commands are fed as inputs to the model computer which calculates the aircraft 
response to be reproduced. These responses, along with TIFS motion sensor 
signals, are used to generate feedforward and response error signals, which 
drive the six controllers on the TIFS (Figure 14). The result is a high fi- 
delity reproduction of the motion and visual cues at the pilot position of the 
model aircraft. More descriptions of the TIFS can be found in Reference 16. 

This experiment made use of the following major features inherent in 
the TIFS aircraft i 

1. Independent control of all six forces and moments by use 
of elevator, aileron, rudder, throttle, direct lift flaps 
and side force surfaces. 

2. Longitudinal and lateral/directional model-following 

systems to provide the evaluation pilot with motion and 
visual cues representative of the simulated aircraft. 

3. Separate evaluation cockpit capable of accepting appro- 

priate pilot controls, displays, and co-pilot assistance. 

(An observer but never co-pilot, was present there). 

4. Evaluation cockpit instruments included standard IFR 

instrument displays featuring an ADI and an HSI as the 

primary instruments, with angle of attack and slideslip 
displayed on indicators to the right hand side of the HSI. 
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Figure 14. TIFS MODEL FOLLOWING SIMULATION 




The vertical and horizontal bars on the ADI displayed 
command information for tracking localizer and glide 
slope, respectively. 

5. Digital magnetic tape recording system to record control 
inputs and appropriate aircraft responses. 

6. Two cassette tape voice recorders for recording evaluation 
pilot comments, and TIFS crew comments. 

7. A video camera pointing out the front window and video 
cassette recorder to record the forward view. 

8. The capability to simulate artificial or cancel actual 
crosswinds up to 15 kts incorporated in the model- 
following system. 

9. Turbulence simulated by playing pre-recorded random 
signals into the model through filters mechanized to 
produce the proper power spectrum of turbulence. (Used 
only on two evaluations by Boeing pilot.) 

10. A signal light located above the ADI and audio signal to 
indicate simulated or actual touchdown of main landing 
gear. 

11. Adjustable transport time delay circuits, available to 
simulate time delay in the pilot’s commands to the eleva- 
tor and aileron controls. 

12. Digital computing equipment to calculate model aero- 
dynamics and evaluate kinematic equations. 
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3.2 SIMULATION GEOMETRY 

,* .j,. 

The TIPS motion was configured to reproduce model motion at the 
evaluation pilot’s eye point. As this was a generic rather than specific 
simulation it was decided to superimpose the TIPS c.g. and Generic Transport 
c.g. Additionally, the cockpits of the two aircraft were superimposed. 
Consequently, no transformations were required from TIPS c.g. to model c.g. 

Approaches were made to touchdown and TIPS wheels and Generic 
Transport wheels were superimposed. This simplified geometry negated the 
requirement for eye position and wheel height transformations. 

3.3 EVALUATION COCKPIT CONFIGURATION 

The evaluation cockpit was configured as illustrated in Figure 13. 
The controls were standard wheel and rudders. No window masking was used. 
Thrust was controlled by four throttle levers tied together and total thrust 
was indicated on a single gage. Asymmetric thrust control was not provided. 

The evaluation pilot’s instrument panel is also shown in Figure 15. 
It was a standard configuration with raw glide slope and localizer data 
driving the flight director needles on the ADI. 

TIPS evaluation cockpit is a dual pilot side by side arrangement. 
For this investigation the right seat was occupied by a NASA flight test 
engineer. The engineer observed all approaches and landings, assisted in con- 
duct of the flight test card, recorded touchdown dispersion, and recorded sum- 
maries of evaluation pilot comments and handling qualities ratings (HQR) to 
provide timely post flight analysis. 

3.4 TIME DELAY COMPENSATION 

The TIFS model-following simulation mechanization has some inherent 
delays or lags due to computation delays, actuator dynamics, and structural 
dynamics. If no time delay compensation is done, the TIFS response will lag 
the model’s response Cdepending on frequency content) by approximately 60 to 
100 msec in pitch, 100 to 140 msec in roll and 100-140 msec in yaw. Two 
methods can be used to eliminate these lags. One is to use feedforward 
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Figure 15 EVALUATION PILOT CONTROLS AND INSTRUMENT PANEL 


signals from the evaluation pilot’s inputs or model actuator inputs directly 
to the TIFS command actuators. This will start to command TIPS surface 
motions before model repsonses have been generated. If these feedforward 
“jabber" signals are properly tuned, one can eliminate the lag without signi- 
ficantly changing the response shape. The second method that can be used is 
to eliminate delays or lags in the model itself so that it responds faster 
than the actual aircraft being modeled. When the TIFS’ lags are added to the 
speeded up model the final response occurs at the proper time. This latter 
method was used in the present simulation because of its ease of implemen- 
tation . 


In the pitch axis a 60 msec pure time delay was used in the analy- 
tical model but not in the model mechanized on the TIFS computer. In addition 
no pitch actuators were used in the TIFS model mechanization. The result was 
a TIFS pitch rate response which matched the time delay characteristics of the 
previously defined analytical model. 

In the heave axis the TIFS direct lift flaps have a higher bandwidth 
than the elevators (approximately 44 rad/sec versus 20 rad/sec), in order to 
slow down the direct lift flaps so that they would be "in-phase" with the 
elevators, a first order lag was inserted in the feedforward commands to the 
direct lift flaps. This lag had a time constant of .038 sec and slowed down 
the flaps to effectively be "in-phase" with the elevator. 

No time delay compensation was done in the lateral/directional axes. 
The lateral/directional responses were considered to have Level 1 flying 
qualities in the last TIFS/ Pitch Rate program without compensation. This was 
again verified in the checkout phase where slight modifications were made in 
the roll axis feel system to make the forces consistent with a transport 
aircraft feel. 

3.5 MODEL FOLLOWING VERIFICATION 

Step input responses were taken in flight for all of the con- 
figurations for comparison to the analytical model responses to validate the 
simulation. Pitch as well as roll and yaw step responses are shown in 
Appendix C. It can be seen that fairly good overlays were obtained. The only 
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characteristic of note is the 3. A Hz oscillation in N Z p that occurs with the 
very abrupt configurations (especially 2 and A). This is from a lightly 
damped structural mode. They appear large, but are being superimposed on a 
fairly small steady state incremental acceleration (~.lg). The magnitude of 
the structural oscillation remained the same with larger inputs. This type of 
structural oscillation also appears when flying in turbulence. It was not 
considered to be a factor in the evaluations as the pilots never noted the 
oscillatory characteristics. 

Typical model following responses during approaches are shown in 
Figures 16 and 17. In these figures the TIFS response is the solid line and 
the model response is the dashed line. 



Figure 16 MODEL FOLLOWING - FLIGHT 893, APPROACH 1, CONFIGURATION B 
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Figure 17 MODEL FOLLOWING - FLIGHT 893, APPROACH 10, CONFIGURATION 28 
3.6 DATA RECORDING 

Quantitative data was recorded on board the TIPS on a 58 channel 
digital recorder. Table 6 is the recording list of this data. Signals 
include TIFS responses (unsubscripted), model responses (M-subscript), pilot 
inputs, control surface positions and performance parameters. Signals which 
are increments from the engage value are indicated by a A. All angular units 
are in degrees, velocity in ft/sec (except as noted), accelerations in g’s. 

Two voice recorders were utilized to obtain evaluation pilot com- 
ments. One was generally used just for the evaluation pilot (plus anyone on 
hot microphone) and the other recorded all crew comments and radio 
transmissions. 

A video camera was installed in the evaluation cockpit and recorded 
the forward field of view. 
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Table 6 

TIFS RECORDING LIST 


CHANNEL 


VARIABLE 

1 

Si ' 

pitch rate, MODEL 

2 

q 

pitch rate, TIFS 

3 

P M 

roll rate, MODEL 

4 

P 

roll rate, TIFS 

5 

r M 

yaw rate, MODEL 

6 

r 

yaw rate, TIFS 

7 

Si - 

incremental angle of attack, MODEL 

8 

Aa 

incremental angle of attack, TIFS 

9 

6 M - 

sideslip, MODEL 

10 

8 

sideslip, TIFS 

11 

Ae M 

incremental pitch attitude, MODEL 

12 

Ae 

incremental pitch attitude, TIFS 

13 

Si 

roll angle, MODEL 

14 

* 

roll angle, TIFS 

15 

AN, 

incremental normal acceleration, pilot, MODEL 

16 

Pm 

AN - 

incremental normal acceleration, pilot, TIFS 

17 

p 

AN, 

Z M 

incremental normal acceleration, CG, MODEL 

18 

AN Z 

incremental normal acceleration, CG, TIFS 

19 

\ 

lateral acceleration, MODEL 

20 

N y - 

lateral acceleration, TIFS 

21 

4v m - 

incremental true airspeed, MODEL 

22 

AV 

incremental true airspeed, TIFS 

23 

• 

h 

altitude rate, TIFS 

24 

ERR-6 Z - 

error in 6 Z (internal model following signal) 
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Table 6 (Cont*d) 
TIPS RECORDING LIST 


CHANNEL VARIABLE 


25 

Y 

- 

flight path angle, TIFS 

26 

V IAS 

- 

indicated airspeed, knots, TIFS 

27 

V M 

- 

total true airspeed, MODEL 

28 

°M 

- 

total angle of attack, MODEL 

29 

“g 

- 

angle of attack, gust component added to Model 

30 

“m 

- 

angle of attack rate, TIFS 

31 

h WH + pulse TD 

- 

altitude of wheels plus pulse at touchdown 

32 

h (press) 

- 

pressure altitude 

33 

Loc. Dev. 

- 

localizer deviation 

34 

GS Dev. 

- 

glide slope deviation 

35 

f ec 

- 

pitch column force 

36 

f aw 


roll wheel force 

37 

f rp 

- 

rudder pedal force 

38 

6 ec 

- 

pitch column deflection, inch 

39 

6 AW 


roll wheel deflection, deg 

40 

6 rp 

- 

rudder pedal deflection, inch 

41 

PLA 


power lever angle, * 

42 

6 e M 

S 

h M 


elevator position, MODEL 

43 

- 

aileron position, MODEL 

44 

- 

rudder position, MODEL 

45 


altitude rate, MODEL 

46 


- 

incremental angle of attack, complimentary 
filtered, TIFS 

47 

%f 


feedforward command to $ 2 , TIFS 

48 

“m 


angle of attack rate, MODEL 
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Table 6 (Cont'd) 
TIFS RECORDING LIST 


CHANNEL 


VARIABLE 

49 


- flight path angle, PILOT 

50 

CONF 

- digital configuration number 

51 

y 

- lateral distance on runway (+ left of center ) 

52 

DME 

- distance from tacan, nautical miles 

53 

6 e 

- elevator, TIFS 

54 

6 a 

- aileron, TIFS 

55 

6 

r 

- rudder, TIFS 

56 

6 

z 

- direct lift flap, TIFS 

57 

5 

y 

- side force surface, TIFS 

58 

6 

- throttle, TIFS 
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Section 4 
DATA 


4.1 INTRODUCTION 

This section presents the data obtained from the evaluation flights. 
Included is a chronology of the evaluations, pilot ratings and performance. 
Pilot comments and approach time histories are presented in the appendices. 
Analysis dealing with the commanded response and time delay/ sensitivity 
aspects of the program are presented in Sections 5 and 6. 

4.2 EVALUATION CHRONOLOGY 

After a set of three checkout flights which took place during late 
November and early December of 1985, the evaluation flights commenced. A 
typical evaluation flight day started with a takeoff from Buffalo and 
proceeding to Niagara where evaluation approaches took place. After a config- 
uration was set up by the test engineer, the system was engaged, responses to 
an automatic step input were recorded for verification, and then control was 
given to the evaluation pilot. The final full stop landing was made at 
Niagara if another flight was going to be flown. A post flight debriefing was 
held during the refueling operation. After the final evaluation at Niagara, 
the TIPS proceeded to Buffalo for the final landing of the day. This 
generally was not an evaluation. On two days a separate non-evaluation ferry 
flight had to be flown from Niagara to Buffalo. 

During December of 1985 all of the evaluation flights for Pilot B and 
most of the evaluation flights for Pilots A and C were flown. In addition the 
pilots from Boeing (D and E) and the pilot from DFVLR (G) had their evaluation 
flights. In January of 1986 the evaluation flights for pilot A and C were 
completed as well as those for the Lockheed pilot (F). Most of these January 
flights were flown at Buffalo because of poor weather and visibility at 
Niagara. 
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A summary of the evaluation flight data in chronological order is 
presented in Table 7. Included in this listing arei 

Flight number 
date 

level of turbulence 
a note if approaches were at Buffalo 


Pilot: A 

Person 

NASA/Langley 

B 

Ball 

Cal span 

C 

Smith 

NASA/bryden 

D 

Ranz 

Boeing 

E 

Higgins 

Boeing 

F 

Hadden 

Lockheed 

G 

Meyer 

DFVLR 

Configuration: CONF 



Pilot Ratings: separately listed for approach portion (APP), flare and 

touchdown (TD), overall (OA). 

Pilot Induced Oscillation Rating: PIOR 

Approach number giving sequence order on flight. Approaches not listed were 

aborted due to traffic. 

Touchdown Performance: 

-li - sink rate, ft/sec 

V - airspeed, knots - IAS 

X - longitudinal runway position (obtained from TACAN-DME) 

(+) long, (-) short, ft (Resolution was 60 ft) (On a few 
flights, indicated by a " — ", the DME was not turned on, 
and approaches to Buffalo did not have this measurement. 

On these approaches Long or Short indicates the touch- 
down was out of the desired zone and was called out by 
the flight test engineer) . 

Y - lateral runway position, (+) left, (-) right, ft. 
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Table 7 

tifsAlared LANDING APPROAC 
EVALUATION FLIGHT DATA SUMMA 


FLT - DATE 

PILOT 

CONF 

PILOT RATING 

PI OR 

TURBULENCE 

APP 

TD 

OA 

891 - 12/4/85 
None 

A 

B 

3 

3 

3 

1 



1 

2 

2 i 

2 i 

1 



2 

2 

3 

3 

■ 

892 - 12/5/85 
None 

A 

24 

3 

■ 

■ 

■ 



21 

3 

■ 

■ 

■ 



27 

3 

5 

5 

■ 



28 

2 

3 

3 

2 



26 

2 

2 i 


■ 



25 

3 

3 

3 

1 


C * ^ 
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Table 7 (Cont’d) 
tifs/Flared LANDING APPROACH 
EVALUATION FLIGHT DATA SUMMARY 













































































































Table 7 (Cont’d) 
TIFS/FLARED LANDING APPROACH 
EVALUATION FLIGHT DATA SUMMARY 


FLT - DATE 

PILOT 

CONF 

PILOT RATING 

PIOR 

APP ♦ 

7D PERFORMANCE 

TURBULENCE 

APP 

TD 

OA 

Bl 

B 

a 

B 

896 - 12 / 1/85 

B 

3 

6 

6 

6 

■ 

1 

n 

in 

0 

8 







■ 

2 

n 

114 

60 

4 

Light 


■ 

■ 



B 


mm 






■ 


5 

5 


3 

1 

110 

1020 

-3 



1 




H 1 

4 

mm 

110 

0 

4 



11111m 

MU' 



II 

5 

D 

111 

0 

1 o 



17+ 

3 

3 

3 

H 

6 

B 

120 

0 

3 



Lead/ 





7 

H 

122 

60 

-2 



Lag 




Hi 

! 

■I 

1 





24 

3 

3 

3 

B 

t 8 

B 

120 

0 

3 







■1 

9 

D 

118 

60 

6 

891 - 12/7/85 

B 

21 

3 

6 

6 

Hi 

1 

B 

116 


0 







■1 

2 

n 

119 


-3 

Light 














27 

5 

8 

8 

5 

3 

i 

125 

Long 

-2 








4 

No 

T.D. 

- In PI 

0 



25A 

1 

5 

5 

n 

5 

B 

1 

121 

Long 

4 

! 



M 



■■ 

6 

a 

126 

— 

-6 



22A 

5 

9 

9 

5 ! 

7 

3 

H 


-8 








8 

5 

H 

— 

-8 

898 -12/10/85 

B 

23 

2 

3 

3 

m 

1 

fl 

H 

1020* 

-5 



1 





2 

mm 

H 

480 

-5 

None 






1 

3 

B 


240 

-2 



28A 

1 

10 

10 

6 

4 


«f!i 


B 




■1 




5 

hee 

BEB 

8022 

m 



26 

3 

El 

m 

B 

6 

B 

120 

0 

B 





■1 

■1 

■B 

7 

a 

120 

480 

n 



8 

3 

2 

3 

B 

8 

B 

122 

0 

5 







■l 

9 

H 

120 

120 

4 



17 

3 

3 

3 

B 

10 

B 

119 

0 

6 







■1 

11 

B 

116 

120 

2 


• Do not count in performance - pilot not set up 
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Table 7 (Cont’d) 
TIFS/FLARED LANDING APPROACH 
EVALUATION FLIGHT DATA SUMMARY 


n t nn tf 



PILOT RATING 



TD PERFORMANCE 


r L 1 • UH 1 1 

pti nT 

rruF 




PI OR 

APP * 






TURBULENCE 

riLU i 


APP 

TD 

OA 

D 

H 

n 

B 

900 -12/12/85 

B 

20 

5 

HI 

5 

H 

1 

2 

120 

0 

-4 








2 

1 

123 

0 

8 

None 




H 

a—Pf 









1 

3 

3 

3 

B 

3 

n 

118 

0 

5 







■ 

4 

■ 

117 

120 

-4 



2 

3 

3 

3 

B 

5 

n 

113 

120 

3 







■ 

6 

■ 

116 

0 

-3 



4 

2 

2 

2 

B 

7 

B 

116 

0 

-2 







■ 

8 

H 

118 

120 

2 



5 

5 

6 

6 

3 

n 

fl 

117 

60 

0 








Kfl 

■ 

116 

480 

-2 

901 -12/12/85 

D 

B 

2 

2 

2 

m 

i 

l 

120 

0 

4 








2 

2 

121 

0 

2 

None 






II 






(Boeing Fit.) 


25 

3 

H 

H 

B 

4 

3 

124 

240 

2 





■ 

■ 

■ 

5 

3 

123 

0 

0 



28 

5 

8 

8 

3 

6 

1 

125 

0 

3 








7 

2 

122 

0 [ 

4 

902 -12/12/85 

E 

B 

3 

3 

3 

m 

1 

B 

126 


-1 








2 

B 

123 

— 

-4 

Light 


1 




■ 

3 

B 

123 


2 

(Boeing Fit.) 






HH 








24 

2 

3 

3 

1 

4 

i 

120 

Long 

-1 







■ 

5 

2 

123 

— “ 

-5 



27 

!B 

8 

8 

5 

6 

Hi 

HH^Ti 

B 

BBS 

■ 




■ 




7 

B 


H 

fl 



B 

3 

2 

3 

B 

8 

i 

118 


-8 







■i 

9 

i 

120 

0 

3 










































































































Table 7 (Cont’d) 

TIF S/FLARED LANDING APPROACH 
EVALUATION FLIGHT DATA SUMMARY 


FLT - DATE 
TURBULENCE 


904 -12/13/85 
None 

(Boeing Fit.) 


with 

moderate 

taped 

turbulence 


905 -12/13/85 
Light 

906 -12/17/85 
Light 



B 

C 


PILOT RATING 


TD PERFORMANCE 


PIOR APP * 






































































































Table 7 (Cont’d) 
TIFS/lFLARED LANDING APPROACH 
EVALUATION FLIGHT DATA SUMMARY 


FLT - DATE 
TURBULENCE 

PILOT 

prior 

PILOT RATING 

PIOR 

APP 4k 

TD PERFORMANCE 

ujnr 

APP 

TO 

OA 

MrT “ 

D 

B 

D 

B 

907 -12/17/85 

G 

i 

3 

3 

3 

B 

1 

l 



4 








2 

2 

wSm 

Long 

6 

Light 







3 


Em 

— 

-2 

(DFVLR Fit.) 






m 


BBB| 






2 

3 

6 

6 

B 

4 

H 

120 

— 

5 

All app’s at 







5 

H 

124 

Long 

0 

Buffalo 








■ 






13 

3 

9 

9 

5 

6 

n 











7 

K 


Bui 




14 

3 

5 

5 

B 

8 

i 

125 

1 ^ 








■ 

9 

i 

120 



908 -12/18/85 

C 

17 

2 

2 

2 


1 

3 

126 

v __ 


Light to Mod. 












WMaisd. 






; 

2 

3 

120 

— 

-11 

All app’s at 












Buffalo 






IBB 






909 -12/20/85 

G 

11 

2 

2 

2 

2 

1 

In 

o T.D. 

- Heav 

y 








2 

No T.D. 

- Heavy 

None 












(DFVLR Fit.) 


12 

B 

6 

6 

3 

3 

B 

115 

0 

2 


j 


■ 




4 

H 

123 

120 

9 



3 

2 

2 

2 

B 

5 

n 


120 

8 







■■ 

6 

D 

Ba 

60 

6 




5 

5 

5 

B 

7 

2 


0 

8 




HH 



■■ 

8 

2 

W3M 

0 

6 



5 

3 

3 

3 

g 

m 

B 

120 

1020 

9 








BEB 


116 

480 

4 



6 

3 

3 

3 

B 

m 

i 


480 

-i 







■i 

o 

mm 


240 

7 



1 

2 

2 

2 

B 

13 

i 


480 

n 








■fl 

14 

i 

nvim 


n 
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EVALUATION 


FLT - DATE 
TURBULENCE 


PILOT 


CONF 


PIU 


APP 


910 -12/20/85 
None 

(DFVLR Fit.) 


10 


911 - 1/7/86 
Light 


26 


23 


28 


912 - 1/7/86 
Light to Mod. 
(Lockheed Flt)| 


All app’s at 
Buffalo 


Incomplete Eval. 


24 


21 


27 


913 - 1/8/86 
None 
(Lockheed Flt)| 


All app’s at 
Buffalo 


27 


25 


22 


Mh 


7 (Cont’d) 
i LANDING APPROACH 
LIGHT DATA SUMMARY 


3T RATING 



TD PERFORMANCE 



PIOR 

APP * 

1 




TD 

OA 

D 

D 

D 

ni 

2 

2 

1 

1 

H 


Long 

4 




2 

l 

H 

9 

5 

5 

B 

3 

i 

120 

Long 

19 



■i 

4 

l 

116 

— 

9 

3 

3 

B 

5 

l 

120 

n 

9 



■ 

6 

l 

119 

■ 

8 

3 

3 

i 

1 

2 

B 

0 

-5 




2 

1 

la 

0 

-5 

5 

5 

■ 

3 

1 

112 

0 

0 

B 

1 

B 

4 

1 

118 

0 

-6 

■fl 

■fl 

■I 

5 

1 

116 

0 

-7 

2 

2 

1 

1 

1 

B 

B 

n 




2 

3 

wEn 


II 




3 

3 

BMW 

■9 

n 


H 

1 

4 

4 

ffifff 

B 

-3 



5 

6 

tn 


-4 

5 

5 

3 

6 

2 

126 


El 




7 

2 

121 

Long 

El 

3 

3 

■ 

8 

3 

120 

— 

-2 

B 

H 

B 

1 

2 

120 


H 

■1 

■I 

2 

3 

125 

Long 

-l 

■ 

n 

H 

3 

4 

125 


-7 

wm 



4 

2 

127 

Short 

-3 

1 

| I 


5 

2 

116 

Short 

-5 

6 

6 

■ 

6 

2 

124 

Long 

fl 


-9 


















































































































Table 7 (Cont’d) 
TIFS/h-ARED LANDING APPROACH 
EVALUATION FLIGHT DATA SUMMARY 


FLT - DATE 
TURBULENCE 

914 - 1/8/86 
Light 

All app’s at 
Buffalo 



915 - 1/8/86 

Light to Mod. 

All app’s at 
Buffalo 

C 

22 

■ 

5 

5 

2 

1 

2 

25 

6 

8 

8 

3 

3 

4 


■ 

■ 

■ 

2 

5 

6 

8 

2 

2 

2 

1 

7 

916 - 1/9/86 

A 

25 

2 

m 


1 

1 








2 

Light 




1 

1 



All app’s at 


28 

3 

8 

8 

6 

3 

Buffalo 







4 



mm 

3 

9 

9 

6 

5 








6 



7 

n 

5 

5 

1 

7 




■a 




8 















































































Table 7 (Cont'd) 
TIFS/FLARED LANDING APPROACH 
DEVALUATION FLIGHT DATA SUMMARY 


FLT - DATE 
TURBULENCE 


PILOT 


CONF 


PILOT RATING 


APP 


TD 


OA 


PI OR 


APP * 


TD PERFORMANCE 


-h 


917 - 1/9/86 

Moderate 

(Lockheed Flt)| 

All app’s at 
Buffalo 


28 


1 

2 


5 

2 


124 

120 


13 


11 


3 

4 


2 

1 


118 

125 


H 


H 


5 

6 


4 

3 


130 

123 


Long 


7 

8 


No T.D. 
120 


In PIO 
Long 


-1 


918 - 1/10/86 
Light 


1 

2 


No 


T.D. 
113 


10 


2 \ 


2 i 


2 2 


3 

4 


1 

3 


110 

116 


14 


5 

6 


125 

128 


11 


7 

8 


117 

121 


12 


9 

10 


Heavy 
0 |-2 


-180 

0 


0 

0 


0 

0 


-1 

1 


-2 

-2 


-6 

0 


No T.D. - In PIO 
No T.D. - In PIO 


At Buffalo ( 

I 1 


11 

12 


1 

3 


121 

121 


Long 


-7 

-6 
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A cross reference between configuration, flight number and pilot is 
provided in Table 8. 


Table 8 

CONFIGURATION - FLIGHT NUMBER/RlLOT CROSS REFERENCE 


CONFIGURATION 


FLIGHT/PILOT 


B 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 
16 

17 

17 + Lead/Lag 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 


891/A, 893/8 , 906/C, 901/D, 902/E, 9 12 A 

891/A, 900/B , 917/F, 907/G 

891/A, 900/6, 907/G 

895/A, 896/6, 909/G 

895/A, 900/6, 915/C, 909/G 

918/A, 900/6, 909/G, (916/A - with +200 ms delay) 

895/A, 905/d , 909/G 

916/A, 896/6, 909/G 

898/6, 915/C, 910/G 

918/A, 910/G 

918/A, 910/G 

918/A, 917/F, 909/G 

918/A, 909/G 

917/F, 907/G 

918/A, 907/G 

904/D 

894/A, 898/6, 908/C, 904/D, 904/E 
914/A, 896/6, 904/D, 904/E 
894/A 
894/A 

894/A, 900/6, 904/1, 904/E 

892/A, 897/6, 906/C, 912/F 

895/A, 893/6, 915/C, 913/F 

895/A, 898/6, 911/C 

892/A, 896/6, 906/C, 902 A, 912/F 

892/A, 916/A, 893/6, 915/C, 901/D, 913/F 

892/A, 898/6, 911/C 

892/A, 897/6, 906/C, 902 A, 912/F, 913/F 
892/A, 916/A, 893/6, 911/C, 901/D, 917A 
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Table 8 (Cont’d) 

CONFIGURATION - FLIGHT NUMBER/iPILOT CROSS REFERENCE 


CONFIGURATION 


FLIGHT/PILOT 


22A 

25A 

28A 


914/A, 897/8 
914/A, 897/8 
914/A, 898/8 


4.3 PILOT RATINGS 

A summary tabulation of the pilot ratings versus configuration and 
evaluation pilot are shown in Tables 9, 10, and 11 for Approach, Flare and 
Touchdown, and Overall ratings, respectively. PIO ratings are tabulated on 
Table 12. 
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Table 9 

PILOT RATINGS, APPROACH 


CONFIGURATION 

A 

B 

C 

B 

3 

2 

3 

1 

2 

3 


2 

2 

3 


3 

3 

6 


4 

4 

2 

4 

3 

4,3* 

5 


6 

2 

6 


7 

4 

4 


8 


3 

2 

9 

3 



10 

2 i 



11 

5 



12 

6 



13 




14 

3 



16 




17 

2 

3 

2 

17 + L/L 

4 

3 


18 

3 



19 

2 i 



20 

2 

5 


21 

3 

3 

3 

22 

5 

7 

4 

23 

6 

2 

5 

24 

3 

3 

5 

25 

3,2 

3 

6 

26 

2 

3 

3 

27 

3 

5 

3 

28 

2,3 

7 

3 

22A 

3 * 

5 


25A 

4 

4 


28A 

2 

7 



PILOT 

D 

2 


2 

4 

2 


4 


3 


5 


* had +200 ms delay inserted inadvertantly 
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Table 10 

PILOT RATINGS, FLARE AND TOUCHDOWN 



* had +200 ms delay inserted inadvertantly 
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Table 11 

PILOT RATINGS, OVERALL 


PILOT 


CONFIGURATION 

A 

B 

C 

D 

E 

F 

G 

B 

3 

3 

3 

2 

3,3 

2 


1 

2 i 

3 




3 

3 

2 

3 

3 





6 

3 

3 

6 





2 

4 

4 

2 

4 




5 

5 

6,9* 

6 





3 

6 

2 

6 





3 

7 

5 

5 





2 

8 


2 

2 




2 

9 

3 






5 

10 

2 i 






3 

11 

7 





8 

2 

12 

9 






6 

13 







9 

14 

4 






5 

16 




2 




17 

2 

3 

2 

4 

2 



17 + L/L 

4 

3 


3 

4 



18 

3 







19 

2 i 







20 

2 

5 


4 

3 



21 

4 

6 

3 



5 


22 

7 

9 

5 



6 


23 

4 

3 

5 





24 

4 

3 

5 


3 

H 


25 

3,4 

5 

8 

4 


i 


26 

2 I 


4 





27 

5 


6 


8 



28 

3,8 


7 

8 


5 


22A 

H 

9 






25A 

4 

5 






28A 

5 


f 






* had +200 ms delay inserted inadvertantly 
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Table 12 
PIO RATINGS 


CONFIGURATION 

A 

B 

1 

1 

1 

2 

1 

3 

1 

4 

1 

5 

3,6* 

6 

1 

7 

. 1 

8 

9 

1 

10 

1 

11 

4 

12 

5 

13 

14 

3 

16 

17 

1 

17 + L/L 

1 

18 

1 

19 

1 

20 

1 

21 

4 

22 

4 

23 

1 

24 

1 

25 

1,1 

26 

1 

27 


28 


22A 


25A 

1 

28A 

4 



* had +200 ms delay inserted inadvertantly 
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4.4 


PILOT COMMENTS 


The reader is referred to Appendix D which presents the complete 
transcribed evaluation pilot comments for each evaluation. Included in these 
comments are some pertinent post flight comments as well as the in flight 
comments. 


4.5 APPROACH TIME HISTORIES 

Time histories of the approaches are presented in Appendix E. These 
time histories cover the final twenty seconds of an approach and show control 
activity and important responses of the aircraft. 
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Section 5 

INTERPRETATION OF THE RESULTS USING MIL-F-8785(C) 

AS A FLIGHT CONTROL DESIGN CRITERIA 

5.1 INTRODUCTION TO ANALYSIS 

The results of the experimental flight test program strongly indicate 
that the response of primary interest to a flight control system designer is 
angle of attack. Angle of attack is the fundamental and therefore, the most 
critical physical quantity associated with the maneuvering flight of an 
airplane. The purpose of an elevator on an airplane is to rotate the airplane 
with respect to the relative wind, therefore changing the angle of attack. 
The resulting change in lift enables the airplane to maneuver, thereby 
altering its flight path. 

A pilot, however, does not normally sense angle of attack and senses 
changes in flight path angle often only through vertical acceleration as a 
surrogate variable or perhaps directly when close to the ground. The variable 
the pilot senses most directly and consciously during maneuvers appears to be 
attitude changes. 

Considerable flight experimentation results suggest that the primary 
response of interest to a pilot is flight path; i.e., the pilot almost always 
manipulates his controls in order to change the direction of the flight of the 
airplane. The variable he senses most often, however, appears to be changes 
in attitude. It can be concluded, and considerable evidence exists to support 
the contention that the pilot wants the airplane to fly in a direction propor- 
tional to the direction being pointed. The key to whether or not changes in 
attitude can be used as a surrogate or substitute for changes in flight path 
angle is the behavior of the angle of attack response of the vehicle (or 
changes in lift if other than an elevator alone is used for maneuvering 
purposes). 

The results of the flight experiments reported upon in this document 
suggest that a pilot prefers a flight control system designed such that he can 
use changes in pitch angle as a surrogate for changes in flight path angle. 
If the correspondence between pitch angle and flight path angle is inadequate, 
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however, there are indications that he will judge the vehicle flying qualities 
directly on the basis of the y (and/or y) dynamic behavior of the vehicle. 

Flying qualities requirements, as defined by MIL-F-8785(C) concen- 
trate on the quality of the angle of attack response of the vehicle. The data 
used to formulate the specification was obtained primarily from conventional 
or angle of attack command airplanes. Because (to a fairly good approximation 
for conventional aircraft), y « ZqO. and because for a majority of the flight 
experiments the pilot was located near the center of rotation (percussion), 
the specification on angle of attack served equally well for y or n 2 , as well. 
A conventional or angle of attack command airplane has an a/6(s) transfer 
function characterized usually by phugoid poles placed relatively close to the 
lower frequency zeros, therefore, there is little residue of the phugoid mode 
in the angle of attack response. In the long term response of the vehicle 
there was little question of the use of pitch angle as a surrogate for flight 
path; since a(t) = 0 after the short term response, Ay = Ae and changes in 
pitch attitude indicated corresponding changes in flight path angle. 

The cop vs n/a (constant CAP) requirement in MIL-F-8785(C) defines the 
region of acceptable short period frequency as a function of n/a and there- 
fore, with the addition of an acceptable c sp range, define the ranges of 
acceptable a(t) dynamic behavior of the vehicle in the short term. The asso- 
ciated acceptance correlation parameter between angle of attack (or in conven- 

®U=0 + 

tional aircraft, n z ) and pitching motions is defined by CAP = fT] . Alter- 

2 z ss sp 

nately, the CAP parameter is defined as CAP = — — = — . For an 

n/a V/g(l/T 02 ) 

acceptable angle of attack (or n z ) response, the dynamics of the pitch atti- 

bp 

tude response is defined as a function of velocity and 1/t 02 . Therefore, the 
ton vs n/a item in MIL-F-8785 (C) describes not only the acceptability of the 
angle of attack behavior of the airplane, but also defines the harmony between 
rotation and translation (pitch and heave) in the short term. 

In the long term, the experimental results reported upon herein 
strongly suggest a requirement that in the long term the parameter of impor- 
tance is correspondence between changes in y and e. The smaller the phugoid 
or low frequency residue in a(t) the more likely will be the acceptance by the 
pilot of the dynamics of the vehicle. 
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It can be shown that many of the recently proposed flying qualities 
indicators correlate with the requirements on a(t). For instance, the pilot 
dislikes pitch “bobble** unless y also “bobbles** in adequate correspondence. 
The pilot generally does not prefer a high phugoid frequency except if the 
vehicle is o command, for then a(t) = 0 in the long term and Ae = Ay. If the 
control system is configured as a rate command system in the short term (i.e., 
no pitch “overshoot”) then 1 /t Q 2 should have a value sufficiently large such 
that the angle of attack response, (which is dominated by 1 /t 02 in this type 
of system), is sufficiently rapid to satisfy the aip vs n/a requirement of 
MIL-F-8785(C). Finally, if the phugoid frequency is very low (or approxima- 
tely equal to lA e ^) then the long term residue in both a and q is small and 
the flying qualities are enhanced. In fact, if the long term modes of motion 
were entirely decoupled from pitch and heave, it appears to make little dif- 
ference whether the flight control system is configured as an angle of attack 
or pitch rate command system, and in fact the pitch rate command system 
appears to be preferred because the pilot has confidence that flight path will 
be in correspondence with the pitch commands. 

5.2 ANALYSIS OF THE RESULTS 

5.2.1 Comparison with Up vs n/a Requirements 

Table 13 summarizes the pilot ratings for the first 14 dynamic 
configurations of the experimental flight test program. Separate evaluations 
were solicited for the approach and the flare/touchdown segments of the 
landing task. 

Figure 18 locates the configurations with respect to the <Dp vs n/a 
constant CAP requirement of MIL-F-8785(C) of the first 14 configurations. 
Configurations 2, 4 and 10 are predicted to be Level 2, because the damping 
ratio was higher than the allowable £ = 1.3, while Configuration 11 and 12 are 
not included because they are higher order. MIL-F-8785(C) accurately predicts 
the average ratings of 7 of the 12 configurations for the approach, and only 4 
of the 12 configurations during flare and landing. The rather poor predictive 
ability of MIL-F-8785(C) strongly suggests the need for criteria improvements 
relating to piloting technique and/or vehicle dynamics command configurations. 
Configurations 13 and 14 indicate strongly that the flare and landing task may 
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be quite demanding, requiring considerably more stringent constraints than 
those outlined in the Category C requirements of MIL-F-8785(C). As shown in 
Figure 18, Configurations 13 and 14 would be rated Level 1 according to the 
Category C requirements and for the approach phase of flight, the C require- 
ment accurately predicts the rating. For flare and landing, the ratings were 
Level 2 and 3, and as shown in the figure, more accurately corresponding to 
the flight precision requirements of Category A flight. 

Table 13 

PILOT RATINGS SUMMARY 


Con figs. 

Approach 

Flare and Land 

Actual 

Mean 

Actual 

Mean 

1 

2 , 3, 3, 3 


2^*3, 3, 3 

2.8 

2 

2,3,3 


3,3,6 

4.0 

3 

3,6,2 


3,6,2 

3.7 

4 

4,2,4, 5 

3.8 

2,2,4, 5 

3.2 

5 

4,5,3 

4.0 

6,6,3 

5.0 

6 

2,6,3 

3.7 

2,3,3 

2.7 

7 

4,4,2 

3.3 

5,5,2 

4.0 

8 

3,2,2 

2.3 

2 , 2,2 

2.0 

9 

3,4 

3.5 

3,5 

4.0 

10 

2 j »2 

2.3 

2-i 3 

2.8 

11 * 

5,5,2 

4.0 

7,8,2 

5.7 

12 

6,4 

5.0 

9,6 

7.5 

13 

3,3 

3.0 

5.1 o 

7.3 

14 

3,3 

3.0 

4,5 

4.5 


Predicted FQ Level 
(MIL-F-8785-C) 


1 

2 

1 

2 

1 

1 

1 

1 

1 

2 








NOTE: THE BOUNDARIES FOR VALUES OF n/Ot GREATER THAN 
100 ARE DEFINED BY STRAIGHT-LINE EXTENSIONS. 

THE LEVEL 3 BOUNDARY FOR n / a LESS THAN 1.0 IS 
ALSO DEFINED BY A STRAIGHT-LINE EXTENSION. 


5.6,7,8 


o 


13,14^,. 

.>o 


2 1.3.9.10 

CO 

CO 

< X 




7 OF 10 CONFIGURATIONS ACCURATELY PREDICTED FOR APPROACH 
4 OF 10 CONFIGURATIONS ACCURATELY PREDICTED FOR FLARE AND LAND 


n/a (g/rad) 


Figure 18. COMPARISON WITH MIL-F-8785(C) SHORT PERIOD REQUIREMENTS 
7 OF 12 CONFIGURATIONS CORRECTLY PREDICTED 




5.2.2 Comparison with aipT 02 vs c Requirements 


Figure 19 locates the configuration with respect to the u)pT e2 vs 

r requirements of the proposed MIL-F-8785(C) Standard and Handbook. This 
sp 

figure accurately predicts the mean flying qualities level for only 4 of the 
12 configurations for the flare and landing task whether it is considered to 
be Category C or if the task is considered to be Category A. 

5.2.3 Time Domain Analysis 

An alternate way to interpret the uip vs n/a requirement of MIL-F- 
8785(C) is to develop a normalized time history response envelope of a(t) 
based upon the range of allowable short period frequency and damping ratio for 
a particular value of n/a. Because, however, there is some question as to 
whether the most appropriate requirement is defined by Category A or Category 
C precision as applied to the flare and landing task, an angle of attack time 
history envelope was developed directly from the flight data found in 
Reference 17, the data base used to formulate the MIL-F-8785(B) requirements. 
The development of this time history envelope, shown in Figure 20 is detailed 
in Reference 11. The data used to develop, the time history envelope covered a 
range of n/a of 3.5 to 11.8, so the envelope would be expected to be less 
broad (i.e., more stringent) than an envelope obtained by transformation of 
the (Dp vs n/a requirements of MIL-F-8785(C) for each value of n/a. 

The flight data from which the a(t) envelope of Figure 20 was devel- 
oped has two major limitations or restrictions. First, the data were 
collected from variable stability aircraft that were mainly angle of attack 
command or conventional aircraft configurations in that the phugoid mode was 
normally lightly damped but dynamically stable. The second limitation is that 
the phugoid dynamics were not specifically documented with the data, so the 
effects of the phugoid mode residues could not be included - the a(t) time 
history response envelope considered only the part of the angle of attack 
response attributable specifically to the short period mode. Because of these 
restrictions it is expected that only those configurations that were angle of 
attack command in the long term, i.e., a(t) approached a constant value as 
t -*■ ® could be evaluated with respect to the envelope. Therefore, 
Configurations 2, 3, 6, 7, 10, 11, 12, and 14 were not evaluated with respect 


5-6 


•> 



a) CATEGORY A 


b) CATEGORY B 


c) CATEGORY C 


Figure 19. CONFIGURATIONS 1-10, 13 AND 14 <unT e2 vs. C sp REQUIREMENTS 

4 OF 12 CONFIGURATIONS CORRECTLY PREDICTED 



TIME (sec) 


Figure 20. TIME HISTORY RESPONSE ENVELOPE 
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to the time history response envelope. The responses of Configurations 1, 4, 
5, 8, 9, 13, and 16 are shown in Figures 21 through 26 with respect to the 
time history envelope criteria. The results indicate that the angle of attack 
time history response envelope criteria is as accurate an indicator of Level 1 
flying qualities as the uv, vs. n/a requirement of MIL-F-8785(C). In fact, the 
a(t) envelope, with the restriction that the responses represent only the 
short period mode, is a direct time domain representation of the oip vs. n/a 
requirement. 

It should be stressed that no single flying qualities indicator is 
fool-proof in the sense that exceptions do not exist. The purpose of an indi- 
cator such as the angle of attack time history envelope is to serve as a guide 
to the flight control system designer. The satisfaction of an indicator 
merely serves to increase the probability that his design will have Level 1 
flying qualities. Many other factors, including individual piloting technique 
have been shown to affect the pilot rating of a configuration. Agreement with 
indicators, however, tend to account for a wide range of piloting technique 
because the flight experiments flown to form a data base were performed with a 
variety of pilots and their associated varying techniques. In these experi- 
ments, Pilot G exhibited a smooth technique with minimum stick activity while 
Pilot B was constantly pumping or "dithering'' the stick. As is to be 
expected, Pilot G tended to prefer the angle of attack command configurations 
while Pilot B tended to prefer the pitch rate command configurations in which 
a pulsing stick input command technique produced approximately the same type 
of response as the smoother step command inputs of the angle of attack command 
configurations. It appears, however, that if the angle of attack response is 
well behaved as defined by the angle of attack response envelope, both 
piloting techniques can be accomodated. As seen from the figures, 
Configuration 1, an angle of attack command configuration in both the short 
and long term and Configuration 8, a pitch rate command configuration in which 
the angle of attack response almost satisfies the criterion, were both rated 
Level 1 by all three pilots who evaluated the configuration. 

One significant exception to the angle of attack response envelope 
rule is noted. Configuration 10, with no phugoid residue in either angle of 
attack or pitch rate was rated Level 1 although the angle of attack response 
did not satisfy the angle of attack response envelope criteria. 
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1 


2 


3 


4 


5 


Figure 22. 


NORMALIZED ANGLE OF ATTACK RESPONSE 
CONFIGURATION 4 
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4 


5 


Figure 23. NORMALIZED ANGLE OF ATTACK RESPONSE 
CONFIGURATION 5 



1 1 * 1 1 1 1 1 1 l * * * * * ■ * 1 * i * ■ 1 * * ■ ■ ■ ■ i ■ ■ ■ ■ ■ ■ ■ ■ ■ i ■ ■ ■ ■ ■ ■ ■ ■ ■ i 
1 2 3 4 5 


Figure 24. NORMALIZED ANGLE OF ATTACK RESPONSE 
CONFIGURATION 8 
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Figure 25. NORMALIZED ANGLE OF ATTACK RESPONSE 
CONFIGURATION 9, 13 



i—i. 1. 1 i i i i i i i i i i i i i i ■ i i ■ ■ i i . i i , , , , , i 

12 3 4 5 

Figure 26. NORMALIZED ANGLE OF ATTACK RESPONSE 
CONFIGURATION 16 
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The lessons learned by the analysis shown in this section is that the 
flying qualities are most likely to be rated acceptable by a range of pilots 
with differing piloting techniques if: 

1. The angle of attack response complies well with the 
angle of attack response envelope and 

2. The system is configured as angle of attack command 
in the long term or phugoid mode. 

It is difficult to express the angle of attack response requirements 
in terms of the other responses of the vehicle. For instance, Configurations 
11 and 12 showed smooth, well behaved short term responses in both pitch rate 
and normal accelerations, and suggest compliance with the short term require- 
ments. The lack of harmony between attitude and flight path, however, is 
demonstrated by the fact that the angle of attack response is not well 
behaved. An equivalent system analysis of angle of attack yielded results 
showing angle of attack to be approximated accurately by an integrator, K/s, 
and has poor flying qualities. 

It appears possible to formulate an angle of attack response time 
history envelope to include the effects of phugoid mode residue. Indications 
are that some residue is allowable, but not very much. It is suggested that a 
Level 1 criteria could allow a variation from the short-period-only steady 
state value of approximately +10X, -5% during the first 5 or 6 seconds of the 
angle of attack response to a pilot step stick command input. 

Phugoid or Long Term Dynamics 

Perhaps more clearly than any other result, the preference is that 
the airplane fly in the direction it is pointing, i.e., in the long term the 
preference is for Ae Ay. This observation seems to be much more important 
than specific values of phugoid frequency and damping ratio. Configuration 1, 
an angle of attack command aircraft, shows large phugoid residue in both pitch 
angle and flight path angle, but the two angles are in harmony, i.e., Ae = Ay 
because a(t) =0 in the long term. The configuration is Level 1. Configura- 
tions 11 and 12, the y command configurations, show large differences between 
changes in attitude and changes in flight path angle and are rated Level 2 and 
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Level 3. The only difference between Configurations 6 and 8 is the better 
harmony between e and y for Configuration 8, which is angle of attack command 
in the long term. Configuration 8 was rated Level 1 by all the pilots. One 
can easily conclude that the most important parameter in the long term for 
manual flight in the flare and landing is the harmony between attitude and 
flight path angle. Therefore, the most important criteria for long term 
behavior is the difference between s and y, which is angle of attack. 

The results of this experimental flight test program suggest that the 
flight control system should be designed to be angle of attack command in the 
long term regardless of the short term command configuration. If attitude 
"hold" is desirable, then the zeros of the angle of attack transfer function 
should be approximately the same as the zeros of the pitch rate transfer func- 
tion, thereby minimizing phugoid mode residue in both the angle of attack and 
pitch rate. To do this would require a second controller, such as a canard, 
used independently of the elevator. Configuration 10, designed as a pitch 
rate command system with no phugoid residue in either a or e was clearly rated 
Level 1 by both pilots for flare and landing even though 1 /t 02 , which domin- 
ates the a(t) response, would be considered too low for a Level l rating. 
This might be interpreted to mean that a small increase in y after e has been 
positioned by the pilot is acceptable. Configuration 2, which was not rated 
as highly but had the same short period dynamics as Configuration 10, showed a 
decrease or "droop" in y after the attitude was set by the pilot. In rate 
command systems, attitude might be used by the pilot as an anticipatory cue 
that y will follow. When this does not happen, a larger value of 
1 Aq 2 appears to be required to show the e-y relationship more quickly to the 
pilot. 


It appears then that if the designer were to try to obtain optimum 
flying qualities for a rate command/attitude hold system, he might require 
either a canard with low control power to minimize the angle of attack residue 
in the long term, or a direct lift flap with high control power to optimize 
the effective value of 1 /t Q2 . 
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5.2.4 Design of Rate Command Systems 


One of the major objectives of the approach and landing flying quali- 
ties experimental program was to try to transform the flying qualities 
requirements into a format more compatible with the techniques of the flight 
control system designer. Although the experimental program did not include 
higher order systems such as would be obtained if a proportional plus integral 
compensation network were used in the design process, the experiment did yield 
valuable information for lower order systems, such as would be attained if 
pole placement or model following methods were used in the design process. 
Because pole placement methods, using either state feedback or deterministic 
observer techniques, do not increase the order of the system response, the 
experiment was designed to define flying qualities/f light control system 
requirements under conditions that would be expected to yield optimum flying 
qualities. 


For a rate command system, one of the short period poles is placed at 
the zero of the constant speed, two degree of freedom q/F s (s) transfer func- 
tions, i.e., at s = -1 /tq 2. The other pole is defined from the relationship 
= P* Because the maximum short period damping ratio is specified to be 
5 S pmax = 1*30, a relationship can be obtained to directly relate the minimum 
value of 1 /t 6 2 and short period frequency to obtain Level 1 flying qualities 
as defined by MIL-F-8785(C). From the transfer function 
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the two relationships are obtained 
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With a maximum short period damping ratio of c S p = 1.30, the relationship 
between short period frequency and lA e2 for Level 1 flying qualities is 
therefore defined by MIL-F-8785(C) to be 


2 ‘ 131 “sp - 1/T e 0 - * 469 “sp 


(4) 
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The above equation specifies a required range of 1 /t 02 as a function 
of short period frequency. For the example of this experiment, w sp = 2.0 with 
a value of 1 /t 0 2 = 0.5, such as defined for the rate command Configuration 2, 
would not satisfy the above relationship and could not be expected to be rated 
Level 1 on a consistent basis. For this case, the angle of attack command 
system of Configuration 1 was indicated, and was shown to be Level 1. The 
rate command Configuration 6 was chosen with w S p = 2.0, 1 /t 0 2 = 0.9, a value 
close to the minimum 1 /t 0 2 requirement of, 

1A 92 = *469 Wsp (5) 

and was given improved pilot ratings. The relationship indicated by Equation 
(4) is depicted graphically in Figure 27. The two lines drawn in the figure 
show the two equality relationships of Equation (4). The angle of attack/ 
pitch rate preference boundary therefore represents an alternate way to define 
a maximum damping ratio of c = 1.3 specified in MIL-F-8785(C) . 



Figure 27. SHORT TERM PREFERENCE AS SPECIFIED BY MIL-F-8785(C) 


5-15 



5.3 


EFFECT OF PILOTING TECHNIQUE 


The evaluation pilots who participated in the program had varied 
backgrounds. Some were past participants in similar experimental flight test 
programs, had complete familiarization with the Cooper-Harper Pilot Rating 
Scale and past experience with pitch rate command flight control systems. 
Others were test pilots familiar with conventional large and medium sized 
aircraft. 


The evaluation pilots displayed what is believed to be a wide range 
of piloting techniques during the performance of this flight experiment. A 
wide range of stick activity is evident even by causual perusal of the raw 
data. An examination of typical flare and landing data for Configuration 1, 
(Figure 28, 29, and 30) show widely varying technique. Pilot G used smooth, 
step-like inputs (Figure 30) while Pilot B superimposed a nearly constant 1 
hertz " dither*’ pumping action (Figure 29). Pilot A flew the airplane in a 
fashion between Pilots G and B in terms of stick activity while relieving the 
required stick force bias by triming the elevator (Figure 28). 

The data tends to support the hypothesis that piloting technique is 
related to pilot preference for a rate command or angle of attack command 
configuration. The most appropriate stick command input for an angle of 
attack command system is a smooth, step or ramp type input as exhibited by 
Pilot G, who tended to prefer the short term angle of attack configurations. 
A comparison of ratings by Pilot G between short term a-command Configuration 
1 and 3, and short term q-command Configuration 2 and 4 clearly tends to 
support the association between technique and command configuration. Pilot B 
behaved oppositely. His natural tendency toward a dithering or pulsing input 
led him to tend to favor the rate command configurations because his inputs 
were more naturally in harmony with the way a rate command system is usually 
flown. Pilot A, on the other hand, displayed both lower frequency step and 
ramp-type inputs as well as a dithering input (at lower frequency than Pilot 
B). He flew configurations with precision and liked almost all of the 
Configurations 1-7. 
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Figure 28. TYPICAL APPROACH, CONFIGURATION 1, PILOT A 







Figure 29. TYPICAL APPROACH, CONFIGURATION 1, PILOT B 
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Figure 30. TYPICAL APPROACH, CONFIGURATION 1, PILOT G 


Power spectral densities of pilot stick activity for the three pilots 
are shown in Figures 31, 32, and 33 (Reference 18). The solid line shows 
stick activity for the approach, while the dashed lines indicate activity 
during the flare and touchdown phase for Configuration 1. These plots clearly 
show the low frequency spectrum of Pilot G, the narrow bandwidth spectrum of 
Pilot B and the dual spectrum input of Pilot A. 

The power spectrums change amplitude with dynamic configuration. 
Figure 34 and 35 show power spectrum for Pilots A and G for Configuration 14, 
rated Level 1 for approach but Level 2 for flare and landing. These specta 
show not only that the low frequency portion of the spectrum is eliminated 
because steady stick inputs are not required to fly the configuration, but it 
also shows that a pilot flys a poorly rated configuration less aggressively 
than one rated Level 1. The dual mode tendency of Pilot A seems to be pre- 
served, whereas Pilot G, who was the smoothest flying pilot for Configuration 
1 still flies less aggressively, choosing a low frequency input spectrum. By 
commanding inputs both near and above the short period natural frequency, it 
can be speculated that Pilot A is evaluating or anticipating both the e and y 
vehicle response, while Pilot G is flying almost totally with respect to 
flight path. 

Although no judgement can be made with respect to preferred or 
superior flying technique, it appears that Pilot A may be most versatile in 
terms of adapting to vehicle dynamic configurations because of his dual-mode 
stick activity spectrum tendency. The task for the flight control system 
designer is to try to design the system that will accomodate the widest 
possible range of piloting technique because it is clear that all these pilots 
changed their piloting technique only minimally as a function of configu- 
ration. They each seemed to evaluate the configuration with respect to their 
individual inherent technique. Configuration 1, a short and long term angle 
of attack configuration and Configuration 8, a higher 1/t 02 short term e 
command configuration but long term a command configuration seemed to best 
accomodate the wide range of piloting techniques used in this program. 
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Figure 31. PSD, CONFIGURATION 1, PILOT A 
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(PSD * FREQ) PILOT STICK ACTIVITY (FORCE) 
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Figure 33. PSD, CONFIGURATION 1, PILOT G 


5-23 



FREQ) PILOT STICK* ACTIVITY (FORCE) 




(PSD * FREQ) PILOT STICK ACTIVITY (FORCE) 


K 

O 

-I 

a 

or 

u 

> 

o 


CL 

CL 

a 


N. 

CD 

0) 

K 

J 

u. 


FES Touc*4*c***j DOTE: 418. RUN: 12. BUBJ: 



Figure 35. PSD, CONFIGURATION 14, PILOT G 
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5.4 


CONCLUSIONS AND RECOMMENDATIONS 


Conclusions 

1. The result of the experiment strongly suggest that the pilot wants the 
pitch attitude to be useable in determining the flight path. The dif- 
ference between attitude and flight path angle, or angle of attack, 
appears to be the most important variable associated with the flying 
qualities of an airplane. Because the pilot generally does not sense 
angle of attack, the angle of attack response should be "well behaved". 

2. The flying qualities specification, MIL-F-8785(C) specifies short period 
frequency requirements of the angle of attack response of the vehicle as 
a function of n/a and defines what is meant by a "well behaved" angle of 
attack response in the short term. 

3. Results suggest that pilots prefer a pitch rate command system only if 
the angle of attack response is well behaved as defined by the Wp, vs n/a 
requirements of MIL-F-8785(C). Because the angle of attack response of 
a pitch rate command system is dominated by I/t^ in the short term, 
pilot preference for a pitch rate command system is a function of 1/t 02 . 

4. In the long term, a primary requirement appears to be that pitch angle 
and flight path angle changes be in harmony, i.e.» 

A relative high frequency phugoid appears not to be objection- 
able if a(t) i 0 in the long term, i.e., if As = by. 

Pilot preference is for pitch rate command/attitude "hold" only 
if precise attitude "hold" also results in precise flight path 
angle "hold". 

5. Results of the flight test program suggest that the Category C precision 
requirements, as they appear in MIL-F-8785(C), are adequate for approach. 
Higher precision, perhaps as defined by the Category A requirements, 
appear to be required for the flare and landing segment of the task. 
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6. The results suggest that a wide range of piloting technique can be accom- 
modated by* 

1. An angle of atack command configuration showing large 
phugoid mode residues in e and y. 

2. A pitch rate command configuration showing very little 
phugoid mode residues in a and y. 

Recommendations 

1. Validation of the lower u>n vs n/o CAP boundary is critical to flying 
qualities and flight control system design. The vehicle for this experi- 
ment should be configured both as an angle of attack and pitch rate 
command system with no phugoid mode residue in either angle of attack or 
pitch rate in order to separate phugoid and short period modal effects. 

2. The effect of phugoid mode residue in a(t) for both angle of attack and 
pitch rate command systems is an important element in the establishment 
of design of criteria requirements. Indications are that a phugoid mode 
effect will be substantial. A flight experiment to define allowable 
phugoid mode residues is highly recommended. 

3. The values of 1 Aq 2 as a function of w S p to establish the boundary 
between angle of attack and rate command preference should be verified. 

4. Additional command configurations directed toward the optimum use of 
multiple controllers is critical to future flight control system design. 

5. A strong flight control system design trend is toward the higher order 
rate command, attitude hold system. Proper design criteria for this type 
of system should be established. 

a) Direct lower order (4th) and higher order (5th) comparisons 
should be made. 

b) Phugoid mode residue limits and 1A© 2 requirements should 
be established. 
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6. A similar study of the y(y) command system configuration should be under- 

taken. Because the y system seems to be gaining in popularity, the 
requirements for excellence should be established. 

7. A flight control system design trend is toward decoupled velocity or 

“auto throttle” behavior. Because of this trend, an investigation of 
decoupled y/tN systems should be undertaken. 

8. Associated with recommendation (7) is the requirement for flying quali- 
ties or dynamic behavior that would best allow a pilot to fly a wind 
shear profile. What FCS configurations satisfy the shear/flying quali- 
ties requirements? 

9. An angle of attack command system is an a(t) regulator and could be 

expected to reduce the vehicle angle of attack response sensitivity in 
turbulence. A rate command system would reduce attitude excursions in 
turbulence. An important flying qualities area of investigation involves 
the design of the command/response configuration in the presence of 
turbulence and wind shears. 
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Section 6 

PREDICTIVE CRITERIA RESULTS AND ANALYSES 


6.1 TIME DOMAIN RESULTS AND ANALYSES 

6.1.1 Introduction 

In the past, time delays have been considered, for the most part, 
autonomous contributors to flying qualities. Attempts were made to specify 
various levels of time delay as being acceptable or not acceptable for a given 
flight task. This has resulted in some confusion in that examples keep 
croping up that seem to be exceptions to the rule,* i.e., large aircraft in 
general and, more recently, the X-29. In both of these cases, Level 1 perfor- 
mance was achieved with time delays that were excessive according to most of 
the literature. 

Experience with the Calspan Learjet has indicated that the effect of 
a given amount of time delay is dependent on a number of parameters? i.e., 
open-loop frequency and damping, command gain, and feel system dynamics. A 
hypothesis was developed that the common denominator was the sensitivity of 
the aircraft to pilot inputs. This program provided an opportunity to examine 
these time delay effects in a formal manner. 

The results show that as pitch sensitivity is increased, tolerance to 
time delay decreases. In fact, by a proper choice of lower pitch sensitivity, 
Level 1 performance could be maintained in the flared landing task with time 
delays from 150 ms to delays in excess of 300 ms. With higher sensitivity, 
configuration with Level 1 performance at 150 ms degraded to Level 2 at 200 
ms. 


The time delay data from this program has allowed development of 
parameters that provide the designer with a method of predicting time delay 
effects using simple measurements of computer generated time histories. These 
time delay and sensitivity metrics have been applied to the previously devel- 
oped time domain predictive criteria and the result is a more mature criteria 
which, when applied to seven flying qualities programs, successfully predicted 
104 of 129 configurations by flying qualities Level (81%). The flying quali- 
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ties of 77 (60%) of these configurations were predicted to within +1 unit of 
the average Cooper-Harper rating. 

6. 1.1.1 Data Base Description 

Data Base Composition 

The original intent of this program was to refine the time domain 
predictive criteria of the previous pitch rate program (Reference 5) by 
expanding the data base of that program and to include time delay effects. 
The first part of this effort was to be a search of existing data bases and 
develop a flight program test matrix that would fill the gaps of existing data 
rather than duplicating previous data. 

The search revealed that much of the existing data base was not 
applicable, and further that very little usable data was practically available 
outside programs previously conducted by Calspan. in order to be usable for 
the time domain predictive criteria refinement i 

• Pilot ratings must be available for the flared landing task . 

• The pilot ratings must apply to the longitudinal dynamics. 

• Time histories must be available, or practical to reconstruct, 
for angle of attack, pitch rate, pilot station normal accel- 
eration, pitch attitude, flight path, controller input, etc. 

• The number of variables in the flight configuration must be 
limited so that cause and effect may be determined. 

In many cases, the pilot assigned an overall rating for approach and 
landing and the flared landing task rating was not available. In other cases 
although the program was directed at longitudinal dynamics, the task was up- 
and-away, not landing. Until recently, programs did not utilize angle of 
attack time histories in the analysis and they were not practically available. 
Other previous programs did not limit the number of variables (one tries to 
get as much out of each flight test configuration as possible) or allowed the 
pilot freedom in selecting sensitivity. Consequently, although most previous 
data can be useful in testing a predictive criteria, very little previous data 
was found to be useful in developing one. 
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During the data search it became apparent that time delay data must 
be generated by this program. It was also determined that creating data for 
criteria refinement alone was not totally efficient. The time delay test 
matrix could be accomplished with a minimum of nine (9) and a maximum of 
eighteen (18) configurations. There were other research efforts that could 
achieve data from this program and the configurations from these other 
programs could also be valuable for criteria development outside the time 
delay area. Consequently, a flight control program was conducted in parallel 
with the criteria development effort and the results of both programs were 
used for criteria development and refinement. 

The test matrix was consequently composed of nineteen (19) flight 
control configurations and thirteen (13) time delay configurations for a total 
of thirty two (32) configurations. 

Data Base Discussion 

The details of the time delay and flight control test matrix are 
discussed in Section 2. The coverage of this matrix was excellent in that 
only three (3) of the thirty-two (32) configurations did not receive multiple 
evaluations. The average coverage was three evaluations, conducted by three 
different pilots, per configuration. This allowed a much more complete 
normalization of the data than is possible in most flying qualities programs 
and enhanced trends which led to more valid analysis. 

Without the increased coverage provided by additional sponsors during 
the flight phase of the program, the analysis would have been only partly 
successful. 

6.1.2 Time Domain Results 

The results of this program, as in most flying qualities research 
programs, was pilot rating data and comments. On-board digital recording 
served primarily to assure that the model was correct and model following was 
achieved. The pilot comments were recorded on video and audio tape as well 
as being manually recorded by the flight test engineer. Comments were also 
monitored by the safety pilots and any inconsistencies with comments versus 
ratings were discussed on the spot. This provided faster learning curves for 
evaluation pilots who were not highly experienced with the rating scale. 
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6. 1.2.1 Pilot Ratings 


Individual pilot ratings are shown in Figure 36 and 37. Figure 36 
depicts the ratings for the flight control configurations and Figure 37 
depicts ratings for the time delay/sensitivity configurations. (In Figure 37 
the configurations are presented out of numerical order to better present 
trends as sensitivity is increased at each level of time delay.) 

A brief study of these figures will demonstrate quite vividly the 
importance of multiple evaluations by multiple pilots in flying qualities 
research. This is shown particularly well in Figure 37, where one knows 
intuitively (and from previous flight test) that performance should decrease 
with increases in time delay and sensitivity from the norm. Pilot B follows 
the expected trend quite closely. Pilot A does for the most part but in the 
case of Configuration 28, he rates it an HQR of 3 in one evaluation and an HQR 
of 8 in another (both evaluations under similar atmospheric conditions). In 
the case of Pilot C, configurations off nominal sensitivity pay a large 
penalty and in some cases overshadow time delay effects. Pilot F in 
Configurations 25, 22, and 28 shows a revese trend by preferring higher 
sensitivities as time delay is increased. 

When one is faced with analyzing data that resulted from single eval- 
uations this sort of pilot scatter can lead to non-results. One is then faced 
with an agonizing study of pilot comments in an attempt to identify causes of 
pilot variance, (there are many which are quite valid) and even if the causes 
are found the result is still a rather poor base for quantitative analysis. 
On the other hand average values from a larger pilot population and repeat 
evaluations can result in very meaningful data. One can manually fair in 
averages on Figure 37 and see definite trends. Figure 38 is a plot of the 
average pilot ratings of Figure 37 and clearly shows valid trends and addi- 
tionally quantifies the results. 

An additional value of repeat evaluations and multiple pilots is the 
ability to recognize anomalous data. An example of this can be seen in Figure 
36. Of the 14 evaluations conducted by Pilot G, seven were at wide variance 
with other evaluations. There are a number of explanations for the variance, 
however, the usefulness of this type data to criteria development is 


6-4 




(NMOCJHOnOi 3UV1d) UDH 


6-5 


Figure 3« INDIVIDUAL PILOT RATINGS, FLARE AND TOUCHDOWN (FLIGHT CONTROL CONFIGURATIONS) 
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Figure 37 INDIVIDUAL PILOT RATINGS, FLARE AND TOUCHDOWN (TIME DELAY CONFIGURATIONS 
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questionable. In a data base of this size one has the luxury of recognizing 
the data as anomalous and discard it for the purpose of criteria quantifi- 
cation. Another case can be seen with pilot A and Configuration 28. In one 
instance an HQR of 3 was assigned and in another instance an HQR of 8. The 
explanation in this case is rather straight forward. If the pilot does not 
make inputs or makes only open-loop type inputs time delays will not affect 
performance. In some cases pilots can get by with this technique. On the 
other hand, when the same pilot is forced into the loop for some reason; i.e., 
turbulence, slightly higher level of aggressiveness, etc., time delays can 
greatly affect performance. It would be unconservative to develop time delay 
criteria based on data resulting from open-loop pilot techniques. 

In any event, all data is shown for the purposes of other types of 
analysis or to allow the reader to take issue with the writer’s selection of 
anomalous data points. 

6. 1.2.2 Discussion of Results 

The primary results for the purposes of this analysis were those 
obtained from the time delay matrix. These results were required to quantify 
time delay and sensitivity effects. The results of the flight control matrix 
were needed to test the refined predictive criteria. 

The pilot ratings of the time delay matrix, Figure 37, were averaged 
and plotted versus time delay at various levels of pitch sensitivity in Figure 
38. Pilot rating scatter is also depicted. As expected, the data showed 
strong correlation between pitch sensitivity, q*/lb (deg/sec 2 /lb), and time 
delay (max slope intercept of pitch rate response to column force input). A 
study of Figure 38 shows that at the minimum time delay (150 ms) slight 
penalties were paid for sensitivities lower and higher than nominal, however, 
all three sensitivities resulted in Level 1, or borderline Level 1, perfor- 
mance. At the high sensitivity any increase in time delay resulted in a rapid 
decrease in flying qualities performance (higher HQR), and Level 1 performance 
could be maintained only out to approximately 180 ms. At the normal sensi- 
tivity, Level 1 performance was maintained out to approximately 220 ms. At 
the low sensitivity Level 1 performance could be maintained out to 270 ms and 
three of the six pilots who evaluated the low sensitivity/high time delay 
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configuration achieved Level 1 and borderline Level 1/Level 2 performance at 
350 ms. 

It appears there are trade-offs available to the designer between 
time delay and pitch sensitivity for the flared landing task. 

6.1.3 Time Domain Analysis 

The first goal' of the analysis was to develop time delay and pitch 
sensitivity metrics that would quantify the effects and have applicability to 
the time domain predictive criteria. The second goal was to refine the 
criteria using these metrics. The final goal of the analysis was to apply the 
refined criteria to a data base acquired from a number of previous programs in 
order to test its viability. 

6. 1.3.1 Time Delay and Sensitivity Metric Development 

Study of Figure 38 provided some insight into the parameters that 
should be addressed in the desired metrics « 

• Slope of sensitivity lines of Figure 38. 

The slope of the nominal sensitivity line is 0.025 HQR/ms. 

This was compared to data from References 1, 3, 4, 5, and 
6. The flying qualities degradation with time delay was 
remarkably consistent and matched favorably with the value 
of 0.025 HQR/ms from Figure 38. A final value of 0.02 
HQR/ms was chosen as a median value for use in the metric. 

a Maximum acceptable time delay; i.e., time delay threshold. 

Experience with various in-flight simulators and using 
these simulators to train test pilot students to identify 
and quantify time delays and effects led the writer to 
initially choose a value of around 100 ms as an acceptable 
threshold for the flared landing task, consequently, 100 
ms was chosen as a metric value. 
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• Sensitivity effect in time delay. 

Figure 38 shows a degradation, at a given time delay, in 

flying qualities with sensitivities greater than 0.42 q/lb 

(deg/sec 2 /lb) and delayed degradation at sensitivities 

less than 0.42 q/lb. Consequently, it was determined that 

• . 

a metric "modifier" of q/lb should be used (0.45 being a 
- , 0.45 

nominal value). 

As a result of the above, the time delay metric developed was« 

q/lb 

TO* = 0.02 (TDq - 100) ^745 

Where 1 TO* = time delay metric (A HQR). 

TDq = time from wheel controller pitch force step 
input to maximum slope intercept of the resul 
tant pitch rate response (ms). 

q/lb = pitch sensitivity measured by the maximum slope 
of the pitch rate responses to a step input 
divided by the pitch force (deg/sec 2 /lb) . 

and when TD q is less than 100 ms, let TO' = 0 



Consideration was given to other time delay measurements; i.e., max 
slope intercept of the angle of attack response, gamma response, N Z p response, 
etc. The use of angle of attack time delay proved promising, however, the 
results when applied to the data base were no better than those of pitch rate 
time delay, and as pitch rate time delay is an established parameter in the 
literature it was chosen as the metric. 


A sensitivity metric was required due to the fact that off optimum 
sensitivity alone can degrade flying qualities. Figure 38 showed that an 
optimum sensitivity for wheel controllers was close to 0.42 deg/sec 2 /lb. Also 
from Figure 38, near Level 1 performance was attainable when sensitivity was 
limited to ±0.2 deg/sec 2 /lb. Consequently, the sensitivity metric developed 
was> i i 


q* 


q/lb - 0.45 
5 ^ 


where* 

q* = sensitivity metric (A HQR). 

q/lb = pitch sensitivity (deg/sec 2 /lb as 
previously defined). 

0.45 = nominal optimum sensitivity 


0.2 = an allowable sensitivity range 


Analysis of the data base of References 1, 3, 4, 5, and 6 disclosed 
that time delay effects could be predicted reasonably well for wheel control- 
lers. This is because there is a larger data base for wheel controllers and 
this particular program used wheel controllers to specifically address time 
delay and sensitivity. Programs using other controllers; i.e., Reference 1, 
LAHOS, using a center stick and Reference 3, NLR, using a side stick, varied 
many parameters including sensitivity. In order to properly identify sen- 
sitivity effects specifically, data must be used from a program structured to 
identify sensitivity effects, specifically. 

An attempt was made to identify sensitivity levels for a center stick 
controller (Reference 1, LAHOS). In the LAHOS program the pilots were allowed 
to choose sensitivity while other powerful variables were also being intro- 
duced. Consequently, it was very difficult to separate sensitivity effects 
alone. 
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A first approximation for center stick controllers resultedi 

q/lb 

TD * = 0.02 (TDq - 100) 

q/lb - 2.8 

q s — n? — 

In the first -approximation it appeared that there is a 6 to 1 
increase in sensitivity required for center stick controllers versus wheel 
controllers. 

Subsequently in the analysis, application of these center stick 
controller metrics with the predictive criteria to the LAHOS data will show 
fair predictions of pilot ratings, however, the results were the worst of the 
data sets tested. This is felt to be due to the lack of adequate sensitivity 
data for center stick controllers, and to the large number of single eval- 
uations of the configurations. 

6. 1.3.2 Predictive Criteria Refinement 

A review of the original time domain predictive criteria (Reference 
5) would showi 

PHQR = 1.7 a* - 1.44 N£p + 0.55 T£ + 3.9 

where i 

PHQR = Predicted Handling Qualities Rating 


- a ss /o I 
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and 


|ot*| + 0.05 
|N* I + 0.05 

I z pl 

a 

when: T * > 6 let T* = 6, and T is defined by 

ct , ot a 

/ 

and when: < 1 let TJ = 0 

(The 0.05 terms prevent unrealistic contributions of T a when a 
or N Z p are very small or zero. The (T a - 1) term indicates that 
T a = 1 sec is a near optimum value for flared landing task.) 

3.9 = A bias term. 

In the pitch rate program of Reference 5 pitch sensitivity was held 
essentially constant at 0.45 deg/sec 2 /lb and time delay was 170 ms except for 
a few configurations. Consequently, the criteria was not valid for configu- 
rations of other sensitivities and time delays. The 3.9 bias term included 
the effects of 170 ms time delay and near optimum sensitivity. 

The first step of the refinement was to simply add the new metrics 
and modify the bias term. A “criteria perfect" flight control system should 
yield an HQR of 2 (similar to Configuration 2-1 of the LAHOS program, 
Reference 1). In the “criteria perfect" airplane the angle of attack response 
is flat at the steady state; i.e., g*= 0, the pilot is located near the center 
of rotation; i.e., N z *p = 0, T a is one second or less; i.e., T a * = 0, the time 
delay is 100 ms or less; i.e., TD* = 0, and if the pitch sensitivity is nomi- 
nal, q* = 0. Consequently: 

PHQR = 0- 0 + 0 + 0 + 0 + Bias = 2 
Therefore, Bias = 2 

The refined time domain criteria then becomes: 

PHQR *= 1.7 a* - 1.44 N Zp + 0.55 T<J + TD* + q* + 2.0 

(With the terms as previously described and TD* and q* modified 

as required to accomodate center stick or wheel controllers) 
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The refined criteria was then applied to a number of data bases from 
this and previous programs (results presented in a subsequent section). The 
intent was to refine the parameters of the criteria as necessary. As it 

turned out the criteria as presented provided predictions well below pilot 

scatter levels and no further modifications were made. 

6. 1.3.3 Application and Limitations of the Time Domain Predictive Criteria 

The criteria was originally designed to apply to flight control 

designs that fall within or close to Level 1 boundaries of the MIL-F-8785 
t*Yi S p and N z /a requirements (the logic was that a flight control designer can 
normally locate the poles and zeros where he pleases, consequently, this 

criteria was not meant to handle low damping ratios or excessively high or low 
frequencies. The problem was that the MIL Spec was originally designed to 
handle conventional pole/zero locations, i.e., conventional short period and 
phugoid roots, and newer designs (pitch rate command, etc.) were meeting MIL 
Spec, requirements as far as short period pitch rate pole locations are 

concerned, but having problems with flared landing performance. The first 
criteria, that resulting from Reference 5., identified the problems associated 
with non-conventional pole/zero locations. It also could identify effects of 

pilot location. An unplanned asset was its ability to identify low short 

period frequency problems by virtue of the T a » term. It would not accomodate 
low damping ratios, time delays, pitch sensitivity variations, or excessively 
high short period frequencies. 

The revised criteria can accomodate time delay and pitch sensitivity 
effects. A by-product of the sensitivity parameter is that excessively high 
short period frequencies will pay a penalty by virtue of resultant higher 

sensitivities, however, this has not been adequately tested. 

There are four areas of known limitations to the criteria two of 
which have been previously discussed i 

1. Lightly damped configurations are not applicable by design. 

2. Divergent configurations are not applicable by design. 
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3. “Decoupled" configurations are not applicable. During this 
program Configuration 5 through 10 were decoupled in nature. 
That is they were angle of attack command and pitch rate com- 
mand (by virtue of the flat steady state response of these 
parameters). This was accomplished by use of an additional 
controller, in, this case direct lift flaps, which provide 
increased lift with out changing the angle of attack. In this 
program, lift could be changed at constant speed and constant 
angle of attack by the direct lift flaps. Use of the criteria 
for decoupled configurations will require further modification. 
Decoupled configurations were shown on the curves but were not 
used for testing the criteria in subsequent sections of the 
analysis. 

4. Pitch controllers other than wheel and center stick are not yet 
applicable due to the lack of sufficient sensitivity data. In 
fact, the sensitivity parameters for center stick are not yet 
sufficiently refined. The NLR data base, Reference 3, could 
not be used in testing the criteria due to lack of sensitivity 
data for the controller used in that experiment. 

It should be re-emphasized that the criteria is applicable to 
the flared landing task. It was not designed to accomodate up-and-away 
tasks such as cruise or air-to-ground tracking, etc. It has only 
limited application to the approach task and has never been tested for 
the non-flared landing task. 
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6.1. A Time Domain Predictive Criteria Testing 


The criteria was tested using flight data from seven programs 
(including the current landing program). Configurations from these programs 
that were not appropriate to the criteria were not tested; i.e., divergent or 
lightly damped configurations, decoupled configurations, controllers other 
than wheel or center stick, and task other than flared landing. 

6.1.A.1 Landing Program Results 


The criteria was applied to the data of the present landing program 
in the form* 

PHQR '= 1.7 ct' - l.AA N4 + 0.55 T4 + TD' + q' + 2.0 

where a*, N£p and Tq[ are as previously defined, and TO* and q* are adjusted 
for wheel controllers, i.e., 

q/lb 

TD * = 0.02 (TDq - 100) ^3 


and 


q/lb - 0.A5 


0.2 


The calculations for the predicted pilot ratings and the actual pilot 
ratings are shown in Table 1A. The results are plotted on Figures 39 and AO, 
Figure 39 reflecting only the time delay matrix data and Figure 40 showing the 
total program results. The curves are a plot of average HQR from the flight 
data (AHQR) versus predicted HQR (PHQR')* 


Figure 39, the time delay matrix (13 configurations), shows the 
application of the criteria to Level 1 flight control systems, both angle of 
attack command and pitch rate command. The results are summarized by» 


CONFIGURATIONS 
PREDICTED TO 
WITHIN 1 HQR 

CONFIGURATIONS 
PREDICTED TO 
WITHIN 2 HQR 

CONFIGURATIONS 
PREDICTED BY FLYING 
QUALITIES LEVEL OR 
WITHIN 1 HQR 

11 of 13 

13 Of 13 

13 of 13 

or 

or 

or 

85% 

100% 

100% 
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Figure 39 


REFINED TIME DOMAIN CRITERIA APPLIED TO THE LANDING PROGRAM 
(TIME DELAY CONFIGURATIONS) 




AHQR (AVERAGE HQR) 



PHQR' (PREDICTED HQR) 


• = "DECOUPLED" CONFIGURATIONS (NOT APPROPRIATE FOR THE CRITERIA) 


Figure 40 REFINED TIME DOMAIN CRITERIA APPLIED TO THE LANDING PROGRAM 
(ALL CONFIGURATIONS) 
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Table 14 

TIME DOMAIN CRITERIA APPLICATION 
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Figure 40, all configurations applicable to the criteria (26 config- 
urations), shows the application of the criteria to many different flight 
control designs and effective time delays. The results are summarized by* 


CONFIGURATIONS 
PREDICTED TO 
WITHIN 1 HQR 

CONFIGURATIONS 
PREDICTED TO 
WITHIN 2 HQR 

CONFIGURATIONS 
PREDICTIONS THAT 
MISSED BY MORE 
THAN 2 HQR 

CONFIGURATIONS 
PREDICTED BY FLYING 
QUALITIES LEVEL OR 
WITHIN 1 HQR 

20 Of 26 

25 Of 26 

1 Of 26 

26 of 26 

or 

or 

or 

or 

77% 

96 % 

4% 

100% 


6. 1.4.2 Pitch Rate Program Results 

The refined criteria was applied to all configurations of the Pitch 
Rate Program (Reference 5), again using wheel controller sensitivity factors. 


Calculations are shown in Table 15 ana Figure 41 shows the application 
of the refined criteria to the original data base. The results are summarized 
by: 


CONFIGURATIONS 
PREDICTED TO 
WITHIN 1 HQR 

CONFIGURATIONS 
PREDICTED TO 
WITHIN 2 HQR 

CONFIGURATIONS 
PREDICTIONS THAT 
MISSED BY MORE 
THAN 2 HQR 

CONFIGURATIONS 
PREDICTED BY FLYING 
QUALITIES LEVEL OR 
WITHIN 1 HQR 

17 of 27 

24 of 27 

3 of 27 

21 of 27 

or 

or 

or 

or 

63% 

89% 

11% 

78% 

(Previous 
Criteria 14 
of 27 or 52%) 

(Previous 
Criteria 23 
of 27 or 85%) 

(Previous 
Criteria 4 
of 27 or 15%) 

(Previous 
Criteria 19 
of 27 or 70%) 


Although most of the pitch rate configurations were flown at the same 
sensitivity and time delay there were a few at greater time delays and off 
nominal sensitivity. Consequently, the refined criteria improved the predic- 
tion results. 
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Figure 41 REFINED TIME DOMAIN CRITERIA APPLIED TO THE PITCH RATE 
PROGRAM 
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Table 15 

LAHOS PROGRAM - TIME DOMAIN CRITERIA APPLICATION 
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6. 1.4.3 LAHOS Program Results 

The LAHOS data' (Reference 1) was the most difficult to utilize due 
primarily to the lack of sensitivity data for center stick controllers. In 
LAHOS the pilots were given freedom to choose pitch sensitivity. At the same 
time there were variations in other major parameters. Optimum sensitivity was 
estimated by fairing straight lines through data that contained many other 
variables. The nominal sensitivity chosen was 2.8 q/lb (deg/sec 2 /lb) as 
opposed to 0.45 q/lb for wheel controllers, or a ratio of approximately 6 to 
1. Sensitivity variation was set at 1.2 versus 0.2 for wheel controllers. 

Consequently, the criteria applied to the LAHOS data was* 

PHQR ** 1.7 a* - 1.44 N£ p + 0.55 T<$ + TO* + q» + 2.0 

wherei a*, N£ p and T a # are as previously defined 

q/lb 

and i TO* = 0.02 (TDq - 100) ^ g 

1 q/lb - 2.8 I 


Pilot trends were very difficult to analyze from the LAHOS data. The 
author feels that this is due in large part to the size of the test matrix. 
Twenty one (21), of the thirty two (32) configurations that were tested for the 
criteria were single evaluations} i.e., one pilot, one evaluation. If the 
reader will review Figures 36 and 37 it will be apparent why single evaluations 
are quite risky for research purposes. If one were limited to a single random 
choice of the pilot ratings from each evaluation of Figures 36 and 37, it 
becomes obvious why some flying qualities programs reach indefinite conclu- 
sions. One could have obtained by random choice, results that show that as 
time delay was increased for a given configuration flying qualities get better! 
In fact, this trend was found in some of the LAHOS data. 
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In any event, all the LAHOS data that was appropriate to the criteria 
was used and no effort was made to attempt to identify anomalous data. The 
results were plotted on Figure 42 with the supporting data in Table 16. The 
results are summarized by: 


CONFIGURATIONS 
PREDICTED TO 
WITHIN 1 HQR 

CONFIGURATIONS 
PREDICTED TO 
WI THIN 2 HQR 

CONFIGURATIONS 
PREDICTIONS THAT 
MISSED BY MORE 
THAN 2 HQR 

CONFIGURATIONS 
PREDICTED BY FLYING 
QUALITIES LEVEL OR 
WITHIN 1 HQR 

13 of 32 

23 Of 32 

9 of 32 

24 of 32 

or 

or 

or 

or 

41* 

12 % 

28* 

75 % 
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Figure 42 REFINED TIME DOMAIN CRITERIA APPLIED TO THE LAHOS PROGRAM 
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Table 16 

LAHOS PROGRAM - TIME DOMAIN CRITERIA APPLICATION 
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6. 1.4.4 Large Aircraft Program Results 

The Large Aircraft program (Reference 4) was used as a source of data 
for the original time domain criteria of the Pitch Rate program, in the 
original application these configurations with time delays in excess of 170 ms 
were excluded as they exceeded the scope of the original criteria. Also, when 
used to test the original criteria, no criteria parameter addressed sensi- 
tivity. 


In this case the refined criteria (modified for wheel controller) was 
applied to thirty four (34) configurations, and included time delays from 170 
ms to 410 ms. Those configurations excluded were those that were divergent 
and those where the pilot ratings were based on the approach and not the 
flared landing task. 


The results are plotted on Figure 43, with supporting data in Table 
17. The results are summarized by» 


CONFIGURATIONS 
PREDICTED TO 
WITHIN 1 HQR 

CONFIGURATIONS 
PREDICTED TO 
WITHIN 2 HQR 

CONFIGURATIONS 
PREDICTIONS THAT 
MISSED BY MORE 
THAN 2 HQR 

CONFIGURATIONS 
PREDICTED BY FLYING 
QUALITIES LEVEL OR 
WITHIN 1 HQR 

18 of 34 

31 of 34 

3 of 34 

24 of 34 

or 

or 

or 

or 

5396 

9196 

996 

7196 

(Previous 
Criteria 5 
of 16 or 3196) 

(Previous 
Criteria 13 
of 16 or 8196) 

(Previous 
Criteria 3 
of 27 or 1996) 

(Previous 
Criteria 10 
of 27 or 6396) 


These results are considered to be a significant improvement over the 
original criteria, particularly considering the much more broad scope of 
application. 
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Figure 43 REFINED TIME DOMAIN CRITERIA APPLIED TO THE LARGE 

AIRCRAFT PROGRAM 
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Table 17 

AIRCRAFT PROGRAM - TIME DOWIN CRITERIA APPLICATION 
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Table 17 (Cont’d) 

LARGE AIRCRAFT PROGRAM - TIME DOMAIN CRITERIA APPLICATION 
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6. 1.4.5 NT-33A/Ames Study Results 

The NT-33A/Ames Study is a yet unpublished program sponsored by NASA 
Ames to investigate differences between ground and in-flight simulation pilot 
ratings. The test matrix consisted of three configurations from the LAHOS 
program (Reference 1), 2-1, 2-7, and 6-1. In addition, two levels of time 
delay (108 ms and 144 ms) were added to the inherent time delay of 100 ms for 
Configuration 2-1 for a total of five (5) configurations. 

This data provided an excellent check on previous testing of the 
criteria (Figure 42) on LAHOS configurations. However, in this case there 
were numerous evaluations by different pilots for each configuration. (A 
comparison of predicted HQR (PHQR*) values from Figure 42 will show slight 
variation due to small differences in the time histories from the Ames program 
and the LAHOS work.) 

As in the previous LAHOS application the criteria used TO* and q* 
parameters adjusted for center stick controller sensitivities. The results 
are plotted on Figure 44 with the supporting data in Table 18. The results 
are summarized byi 


CONFIGURATIONS 
PREDICTED TO 
WITHIN 1 HQR 

CONFIGURATIONS 
PREDICTED TO 
WITHIN 2 HQR 

CONFIGURATIONS 
PREDICTED BY FLYING 
QUALITIES LEVEL OR 
WITHIN 1 HQR 

4 of 5 

5 of 5 

4 of 5 

or 

or 

or 

80* 

100* 

80* 


These results are a significant improvement over the previous testing 
against LAHOS configuration. The author feels that this is due in large part 
to the relatively high number of repeats which was made possible by a smaller 
test matrix. 
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AHQR (AVERAGE HQR) 



10 987654321 

PHQR' (PREDICTED HQR) 

Figure 44 REFINED TIME DOMAIN CRITERIA APPLIED TO THE NT-33/ AMES STUDY 
(NOT YET PUBLISHED) 


I 
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Table 18 

T-33/AMES STUDY - TIME DOMAIN CRITERIA APPLICATION 


M 
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6. 1.4.6 SST Program and X-29 Data Results 

The majority of SST data (Reference 12) resulted from divergent 
configurations that are inappropriate for the criteria, however, three 
configurations were stable and were applied to the criteria. (These three 
configurations were also tested on the original time domain criteria.) The 
refined criteria applied in this case used wheel controller sensitivity 
parameters. 

Some unpublished X-29 data was available for comparison and was 
applied to the criteria using center stick controller sensitivity parameters. 

The results are plotted on Figure 45 with the supporting data in 
Table 19. The results are summarized by* 


CONFIGURATIONS 
PREDICTED TO 
WITHIN 1 HQR 

CONFIGURATIONS 
PREDICTED TO 
WITHIN 2 HQR 

CONFIGURATIONS 
PREDICTED BY FLYING 
QUALITIES LEVEL OR 
WITHIN 1 HQR 

5 of 5 

5 of 5 

5 Of 5 

or 

or 

or 

100* 

100* 

100* 

(Previous Criteria 
Results 2 of 3 
or 67*) 

(Previous Criteria 
Results 3 of 3 
or 100*) 

(Previous Criteria 
Results 3 of 3 
or 100*) 


6. 1.4.7 Criteria Testing Summary 

The overall results of the time domain criteria testing are presented 
in Table 20. A total of 129 configurations were tested, of these, the pilot 
ratings of 77 (60*) were predicted to within one HQR, 113 (88X) were predicted 
to within two HQR and 16 (12%) had predictions that missed by more than two 
HQR. Predictions by flying qualities level, or within one HQR, were accurate 
on 104 of the 129 (or 81*). 
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Figure 45 REFINED TIME DOMAIN CRITERIA APPLIED TO THE SST PROGRAM 
AND X-29 UN-PUBLISHED DATA 


Table 19 

SST PROGRAM AND X-29 DATA - TIME DOMAIN CRITERIA APPLICATION 




Table 20 

REFINED TIME DOMAIN CRITERIA TESTING RESULTS 
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129 77/60% 113/88% 16/12% 104/81% 



6.2 


FREQUENCY DOMAIN RESULTS AND ANALYSES 


In order to provide additional insight into the pilot rating trends 
obtained in this program, a number of frequency domain flying qualities 
criteria were applied to the configurations of this study. Linear models of 
these configurations were used to generate pitch rate, pitch attitude, and 
attitude rate closed-loop frequency response data. These data were then used 
to assess the utility of a particular criterion with respect to predicting the 
flight results. It should be noted, however, that the following analyses do 
not provide an in-depth assessment of each configuration’s flying qualities. 
Instead, they serve as an attempt to describe the overall pilot rating trends 
obtained with this data base. 


Four frequency domain predictive techniques were considered; lower 
order eqivalent systems, pitch attitude bandwidth, altitude rate bandwidth, 
and Neal-Smith pitch attitude pilot lead compensation. 


6.2.1 Low-Order Equivalent Systems 


The configurations of this study were evaluated using low-order 
equivalent systems analysis (Reference 19). Closed-loop frequency responses 
of each configuration were fit with one of the following low-order modelst 


-TS 


es 


k q (s + 1 A )e 
q _2 

s 2 + 2c us + to* 
sp sp sp 


V s + 1/T o )e 


-TS 


es 


s 2 + 2r w s + o>* 
sp sp sp 


(Frequency range = 0.1 - 10.0 rad/sec) 


Selection of a low-order model for a given configuration was based upon its 
high-order, closed-loop frequency response shape. Configurations without 
phugoid resonant peaks in their pitch rate frequency response were matched to 
the pitch rate low-order transfer function with a fixed value of I/t^. Those 
configurations exhibiting phugoid resonant peaks in their pitch rate frequency 
response were matched to the low-order alpha transfer function. This stemmed 
from the fact that the equivalent system models did not include the phugoid 
response in their transfer functions and any phugoid dynamics of significant 
magnitude would result in unacceptably high levels of mismatch. 
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The low-order equivalent system results for the configurations of 
this study are presented in Table 21. Although satisfactory low-order fits 
were made to most of the configurations, the results did not adequately 
discriminate between the good and bad configurations. All equivalent short 
period frequencies were Level 1 (Reference 7). The equivalent short period 
damping ratios were Level 1, with the exception of the overdamped configu- 
rations (2, 4, 6, 8, and 10) which were Level 2. The equivalent time delays 
were primarily in the Level 2 range, except for the time delay configurations 
with additional delays which were Level 3 and below (Level 4). 

The results in Table 21 were to be expected, since the configurations 
of this study were originally designed to be low-order models. In addition, 
the low-order models did not include the phugoid response. The only addi- 
tional insight is provided by the equivalent time delay results, which show 
degradation of all configurations to the Level 2 range in the best cases. 
This does not correlate well with the documented pilot ratings comments. 
Overall, only 38X of the configurations were accurately predicted by level 
using the equivalent system approach. 

6.2.2 Pitch Attitude Bandwidth Criterion 

The pitch attitude bandwidth criterion (Reference 20) was used to 
provide additional insight into the pitch characteristics of the configuration 
in this study. The criterion was applied to the closed- loop frequency 
response of each configuration with the bandwidth defined as the crossover 
frequency at which the phase margin is 45° or the gain margin is 6 dB, 
whichever frequency is lower. In addition, an estimate of the time delay of 
each configuration is calculated to bet 

x = — ((> + 180°)/(57.3u>) 

where u is defined to be two times the frequency at 180° and {> is the phase at 
id. This criterion was formulated using the configurations from Reference 1. 

The pitch attitude bandwidth results for the configurations in this 
study are presented in Table 22. The level ratings in Table 22 were based 
upon the level boundaries established in Reference 20, which were determined 
using the pilot ratings for the fighter aircraft configurations from Reference 
1. As a result, only 44% of the configurations were accurately predicted by 
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Table 21 

LOW-ORDER EQUIVALENT SYSTEM RESULTS 


ZETA I LEVEL 





OMEGA 

(R/S) LEVEL' 


TIME 

DELAY** 

(SEC) 



0.148 

0.248 

0.348 

0.148 

0.248 

0.348 

0.160 

0.360 

0.360 

0.360 


LEVEL 


2 

2 

2 

2 

2 

2 

2 

2 

2 

2 


OVERALL PILOT 
PREDICTED RATING 
LEVEL LEVEL 



% 1 A ®2 terms fixed, l/i a terms free. 

** Referenced to stick force. Feel system and actuator included. 

*** Configurations B, 21-28 fit in the frequency range of 0.3-10.0 rad/sec. 
**♦* Level 4 is worse than Level 3 
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Table 22 

PITCH ATTITUDE BANDWIDTH CRITERION RESULTS 


CONFIG 

NO 

BANDWIDTH 
(Phase Margin) 
(rad/sec) 

BANDWIDTH 
(Gain Margin) 
(rad/sec) 

TIME DELAY 
(sec) 

PREDICTED 

LEVEL* 

1 

2.36 

2.86 

0.114 

2 

2 

3.01 

3.65 

0.109 

1 

3 

2.38 

2.87 

0.114 

2 

4 

2.96 

3.65 

0.109 

2 

5 

2.07 

2.59 

0.114 

2 

6 

2.24 

3.13 

0.109 

2 

7 

2.13 

2.65 

0.114 

2 

8 

2.12 

3.08 

0.109 

2 

9 

2.38 

2.87 

0.114 

2 

10 

3.01 

3.65 

0.109 

1 

11 

1.16 

1.44 

0.226 

3 

12 

1.23 

1.11 

0.308 

3 

13 

1.53 

2.08 

0.119 

2 

14 

1.36 

2.38 

0.108 

2 

16 

3.06 

3.24 

0.117 

2 

17 

2.44 

2.77 

0.126 

2 

18 

2.48 

2.81 

0.126 

2 

19 

2.52 

2.84 

0.126 

2 

20 

2.60 

2.90 

0.126 

2 

B 

2.20 

2.71 

0.114 

2 

21 

2.20 

2.71 

0.214 

3 

22 

2.20 

2.71 

0.314 

3 

23 

2.20 

2.71 

0.114 

2 

24 

2.20 

2.71 

0.214 

3 

25 

2.20 

2.71 

0.314 

3 

26 

2.20 

2.71 

0.114 

2 

27 

2.20 

2.71 

0.214 

3 

28 

2.20 

2.71 

0.314 

3 

17+L/t 

3.13 

3.21 

0.127 

2 

22A 

3.13 

3.21 

0.317 

3 

25A 

3.13 

3.21 

0.317 

3 

28A 

3.13 

3.21 

0.317 

3 


PILOT 

RATING 

LEVEL 



Level determined by boundaries established in Reference 18 
for category flight phases. 
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level. Considering the significant differences between fighter and transport 

aircraft, the pitch attitude bandwidth results in Table 22 were combined with 

similar results from the configurations of Reference 5 to produce suggested 

boundaries for transport aircraft configurations. These results are presented 

in Figure 46 and include 58 generic transport configurations between the two * 

programs. The suggested level boundaries in Figure 46 result in 70 % of the 

configurations being correctly predicted by level, a marked improvement from 

the boundaries used in Reference 20. The location of the level boundaries in 

Figure 46 relative to those used for fighter aircraft in Reference 20 bring 

out two important points. First, pilots flying transport aircraft can accept 

a lower bandwidth in pitch and still achieve desired performance. Second, 

transport aircraft can have higher values of allowable time delay relative to 

fighter aircraft, however, this is dependent on the types of controllers used 

and sensitivities selected. 

6.2.3 Altitude Rate Bandwidth Criterion 

Control of altitude rate in the flared landing task was analyzed for 
the configurations of this study through a single-loop closure technique used 
in Reference 5. A single-loop closure was performed on altitude rate around 
the following pilot model* 

Y p ^= K p e“ TS (t = 0.25 sec) 

The closed-loop performance was determined through the use of a Nichols chart 
over the frequency range of 0.1 to 10.0 rad/sec. The altitude rate bandwidth 
was the frequency at which 3 dB of closed-loop resonance corresponded with 
-90° of closed- loop phase. This procedure was performed on a Nichols chart 
through the manipulation of the pilot’s bandwidth and gain. This technique, 
although obtaining limited success in Reference 5, showed promise when applied 
to the configurations of this study. 

The altitude rate bandwidth and corresponding pilot gain results are 
presented in Table 23. A correlation was made between the amount of pilot 
gain required to achieve the stated closed-loop performance and the altitude 
bandwidth. These results are presented in Figure 47. As expected, configu- 
rations that allow the pilot to increase his gain in the control of altitude 
rate (for a given bandwidth) were generally considered to be Level 1. Figure 
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Table 23 

ALTITUDE RATE BANDWIDTH RESULTS 


CONFIGURATION 

NO. 

PILOT GAIN 

BANDWIDTH 

(rad/sec) 

1 

2.32 

1.65 

2 

4.05 

1.43 

3 

2.40 

1.68 

4 

3.57 

1.35 

5 

1.49 

1.35 

6 

1.75 

1.13 

7 

1.58 

1.43 

8 

1.49 

1.00 

9 

2.48 

1.65 

10 

4.05 

1.38 

11 

0.61 

0.83 

12 

0.78 

1.05 

13 

0.68 

1.18 

14 

0.78 

1.03 

16 

3.55 

1.93 

17 

1.56 

0.78 

18 

1.71 

0.83 

19 

1.90 

0.95 

20 

2.17 

0.95 

B 

1.97 

0.90 

21 

1.74 

0.83 

22 

1.54 

0.75 

23 

3.30 

0.90 

24 

2.90 

0.83 

25 

2.58 

0.75 

26 

1.32 

0.90 

27 

1.17 

0.83 

28 

1.02 

0.75 

17+L/t 

2.10 

1.20 

22A 

1.60 

0.98 

25A 

1.30 

0.90 

28A 

1.71 

1.00 
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47 also shows proposed level boundaries that result in 63% of the configu- 
rations being accurately predicted by level, however, only 53% of the time 
delay/sensitivity configurations were predicted by level using this technique. 

6*2.4 Neal- Smith Analysis 


The Neal-Smith criterion (Reference 21) was applied to the data base 
and correlated with the results from Reference 5 to provide insight into the 
pitch characteristics of generic transport aircraft. The criterion is based 
upon a single-loop closure performed on pitch attitude using the following 
pilot model i 


Y 




(T lead s + *^ T iag s + D 


The pilot model operates on a pitch-attitude error signal that is the dif- 
ference between the commanded attitude and the aircraft attitude. The pilot, 
through the parameters he is observing, tries to achieve a certain "standard 
of performance" which is defined by a specified closed-loop bandwidth. The 
bandwidth is defined by the 90° closed-loop phase requirement. At frequencies 
below the bandwidth, the pilot attempts to minimize steady-state pitch atti- 
tude tracking errors as defined by a minimum low-frequency droop (typically 
no more than -3 dB). The pilot also attempts to minimize the closed-loop 
resonant peak, which minimizes oscillatory tendencies in pitch attitude. The 
pilot model is adjusted so that the -3 dB droop and the -90° of the closed- 
loop phase conditions are met for a given bandwidth while the closed-loop 
resonance is minimized. These parameters then provide a measure of the 
compensation with which the pilot closes the loop in pitch attitude. After 
the closed-loop conditions are met, closed-loop resonance and pilot compen- 
sation are plotted on a Neal-Smith parameter plane and correlated with the 
pilot ratings for the flared landing task. 


The results of the Neal-Smith criterion are presented for each 
configuration in Appendix F. A bandwidth frequency of 2.0 rad/sec appeared to 
best represent the configurations of this study (in terms of predicting the 
flying qualities levels). This correlates with the pitch attitude bandwidth 
results in Reference 5. Analysis of the data showed that variations in pilot 
rating appeared to be dependent upon the amount of pilot lead or lag required 
to achieve the closed-loop requirements. The closed-loop resonance at 2.0 
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rad/sec was Level 1 for all configurations and was not a factor in the pilot 
Tatings. 


A distinct correlation was made between pitch attitude pilot compen- 
sation and pilot ratings, as shown in Figure 48. The results of Figure 48 
show that pilot lead compensation less than or equal to approximately 25° is 
required for Level 1 flying qualities of transport aircraft in the flared 
landing task. In addition, the results show a linear degradation in pilot 
ratings for increasing pilot compensation. Configurations from Reference 5 
and selected configurations for Reference 4 support this trend. However, of 
the time delay/sensitivity configurations only 54% were predicted by level and 
54% predicted to within +1 HQR. This trend and correlation was not evident in 
the Neal-Smith analyses performed in References 4 and 5. It is also 
interesting to note that the 25° "elbow** for Level 1 flying qualities in 
Figure 48 is consistent with the findings in Reference 22 regarding pitch 
attitude pilot compensation. Although pitch attitude is not the only control 
variable in the flared landing task, it is an essential part of any closed- 
loop flared landing technique employed by pilots. Through the combination of 
configurations selected from these three data bases, 73% of the configurations 
were predicted by level and 65% were predicted within 1 pilot rating. 



0NI1VU lOlld 
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6.3 


COMPARISON OF TIME DOMAIN CRITERIA AND FREQUENCY DOMAIN CRITERIA 
PERFORMANCE 


The time domain criteria was tested against seven sets of data and 
129 configuration which included fighter type aircraft with stick controllers. 
The frequency domain criteria were tested against one to three sets of data 
and 58 configurations, the majority of which were contained in the present 
landing program and the previous pitch rate program and which consisted of 
medium transport aircraft with wheel controllers. All of the 58 frequency 
domain configurations were included in the 129 configurations of the time 
domain analysis. In order to obtain a direct comparison between criteria the 
specific 58 frequency domain configurations were pulled from the time domain 
results. The comparison of the frequency domain and time domain criteria 
results of these configurations are shown in Table 24. 

The low-order equivalent systems technique predicted 38% of the 
landing program configurations by level as compared to 100% prediction by 
level of the time domain criteria. The low-order equivalent systems method 
does not appear applicable to the flared landing task. 

The pitch attitude bandwidth technique was only 44% successful in 
level prediction using previous boundaries but when new boundaries were drawn 
based on 58 configurations of the landing and pitch rate programs 70% of the 
configurations fell inside the revised boundaries. This result is promising, 
however, washout prefilter effects could not be predicted and more data is 
required in the area of high bandwidth and higher time delays to more accu- 
rately describe the boundaries. By comparison the time domain criteria 
predicted 89% of these configurations by level. 

The altitude rate bandwidth was applied to the configurations of the 
landing program. New boundaries were established which allowed 63% of the 
configurations to fall within levels, however, time delay and sensitivity 
configurations were not accurately predicted and 3 of 6 Level 3 configurations 
were missed. By comparison the revised time domain criteria predicted 100% of 
the configurations by level. 
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Table 24 


COMPARISON OF FREQUENCY DOMAIN AND 
TIME DOMAIN PREDICTIVE TECHNIQUES 


PREDICTIVE 

METHOD 


Low Order 
Equivalent 
Systems 


Pitch Attitude 
Bandwidth 
Using Fighter 
Boundaries 


Pitch Attitude 
Bandwidth 
Using New 
Boundaries 


Altitude Rate 
Bandwidth 
Using New 
Boundaries 


Neal-Smith 
Pilot Lead 
Compensation 
Using 
Suggested 
Boundaries 


LANDING 

PROGRAM 

(32 Configs Total) 


38% 

Predicted by Level 
(32 configs tested) 


44% 

Predicted by Level 
(32 configs tested) 


PITCH RATE 
PROGRAM 

(27 Configs Total) 


LARGE AIRCRAFT 
PROGRAM 

(34 Configs Total) 


Predicted by Level (58 configs tested) 
Low frequency without Filters 
not accurately predicted 


Predicted by Level (32 configs. tested) 
Time Delay and Sensitivity 
not accurately predicted 




73% 

Predicted by Level 
65% within +1 HQR 
19% missed by >2 HQR 
4 Level 3’s missed 
(26 configurations 
tested) Time Delay 
and Sensitivity not 
accurately predicted 


77% 

Predicted by Level 
60% within +1 HQR 
14% missed by >2 HQR 
no Level 3*s missed 
(22 configurations 
tested) 


75% 

Predicted by Level 
50% within +1 HQR 
25% missed by >2 HQR 
no Level 3's missed 
(4 configurations 
tested) 


Refined 
Time Domain 
Criteria 


100 % 

Predicted by Level 
77% within ±1 HQR 
4% missed by >2 HQR 
no Level 3*s missed 
(26 configurations 
tested) 


78% 

Predicted by Level 
63% within ±1 HQR 
11% missed by >2 HQR 
2 Level 3’s missed 
(22 configurations 
tested) 


75% 

Predicted by Level 
50% within ±1 HQR 
0% missed by >2 HQR 
no Level 3*s missed 
(4 configurations 
tested) 
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The Neal-Smith criteria which was applied to configurations from the 
landing program ( 26 configurations), pitch rate program (22 configurations) 
and the large aircraft program (4 configurations) showed significant correla- 
tion with pitch attitude pilot compensation. Boundaries could be drawn 
(Figure 48) that allowed correct predictions of 75% of the configurations by 
level and 62% within +1 HQR. This is the most promising of the frequency 
domain techniques tested, however, time delay/sensitivity configurations of 
the landing program were not accurately predicted. By comparison the revised 
time domain criteria predicted 88% of the above configurations by level, 69% 
within +1 HQR and of the time delay/sensitivity configurations 100% were pre- 
dicted by level and 85% were predicted to within +1 HQR. 

In the above comparisons the frequency domain criteria were applied 
to primarily two basic data sets containing 58 configurations. Assumptions 
were then made? i.e., Neal-Smith pilot bandwidth, etc., and boundaries were 
drawn that best matched the criteria with the data, i.e., empirical boundary 
definition. These results were then compared to the revised time domain 
criteria which was applied to the same data. These comparisons are shown in 
Table 24. The time domain criteria was empirically developed from the pitch 
rate program, and the revised time domain criteria from the landing program. 
The criteria was then tested on 5 other data sets (Table 20) for a total of 
129 configurations that included fighter and transport aircraft. As described 
in the previous paragraph on a one to one comparison the revised time domain 
criteria was more accurate than the Neal-Smith pilot lead compensation method 
(the most promising of the frequency domain methods) by all comparison metrics 
and considerably more accurate with time delay and sensitivity configurations. 
In fact, there has been no other criteria observed that will account for 
sensitivity. 

When the refined time domain criteria was applied to the additional 
data sets the overall results of 129 configurations listed in Table 20 were« 
predicted by level - 81% 
predicted within +1 HQR - 60% 
configurations missed by more than 2 HQR - 12% 
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These total results were more accurate than the best of the frequency 
domain methods that were empirically fitted to a 58 configuration data base. 
The time domain criteria has fewer limitations, and is especially accurate for 
predicting time delay and sensitivity effects. 

6.4 WASHOUT INVESTIGATION 

Configurations 17 through 20 were evaluated to determine the effect 
of a washout prefilter on a rate command-attitude hold system. Configuration 
17, which was identical to Configuration 1-2-2 of the Reference 5 TIFS/^itch 
study, was chosen as a typical rate command system. It had recieved pilot 
ratings of 5|, 7, and 8 in that reference study. Washout pre filters with 
frequencies of .05, .1, and .2 rad/sec were added with Configurations 18, 19, 
and 20, respectively. 

It was postulated that Configuration 17 would again receive border- 
line Level 2-3 ratings as it did in the previous program, and then the poten- 
tial improvement in flying qualities with various degrees of washout in the 
command path could be evaluated. As it turned out in this program, this con- 
figuration received pilot ratings all between 2 and 4 with an average of 2.6 
from the five separate pilot evaluations. The maximum washout configuration 
(20), with an = .2 rad/sec, was evaluated by four pilots and also received 
pilot ratings between 2 and 4 with an average of 2.8. The only significant 
pilot comments (see Appendix D) were that the no washout configuration (17) 
held flight path better in the approach phase, but that the maximum washout 
configuration (20) felt more natural with the required aft forces, in the 
flare. 


The predictive criteria in this study and Reference 5 indicate that 
the baseline rate command system, Configuration 17, should indeed be rated 
near the Level 1-2 border as it was in this program. The results from this 
study also show that the addition of a washout pre-filter does not degrade the 
flying qualities but just requires the pilot to hold aft forces in the flare. 
A better investigation into washout filters would require a baseline rate 
command system that was definitely a Level 2 or 3 configuration. 
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Calibration steps for Configuration 1-2-2 of Reference 5 were re- 
examined and found to be the same as those for Configuration 17 of the present 
study. It is not known why this configuration was rated worse in the previous 
study, other than the pilot comments indicate a floating overcontrol tendency 
in the flare which did not seem to bother the pilots in the present study. 

6.5 CONCLUSIONS 

Conclusions presented here relate to the predictive criteria results. 
Conclusions from the interpretation of the results using MIL-F-8785(C) as a 
flight control design criteria are presented in Section 5.4. 

6.5.1 Time Delay and Sensitivity Considerations 

There is strong and consistent correlation between time delay effect 
on pilot performance and pitch sensitivity for the flared landing task 
(observations from other programs and experience in in-flight simulators 
indicates strongly that a similar correlation exists for the roll axis and 
other tasks as well). This correlation has been quantified to the extent that 
the designer now has metrics that will allow him to predict the effects on 
flying qualities of various time delay and sensitivity values. 

The results of this program also call attention to the importance of 
proper sensitivity selection. Changes in flight control schemes, i.e., addi- 
tion of prefilters, changes in feedback gains, or any factor that will change 
the shape of the time response, requires proper adjustments of command gain to 
maintain sensitivity within desired bounds. For wheel controllers in the 
flared landing task these bounds are now well defined and the designer has 
metrics available that allow him to quantify, in flying qualities terms, the 
effect of design changes. For center stick controllers the bounds are not yet 
as well defined. For other controllers more data is required. 




6.5.2 Time Domain Criteria Considerations 


The time domain criteria is considered to be sufficiently refined to 
be of use to the designers for flying qualities performance predictions of 
pitch flight control systems in the flared landing task. 

The accuracy of the criteria has been tested against a significant 
number of diverse programs and flight control system designs. In addition, it 
includes significant closed-loop items that previous criteria could not 
account for (sensitivity factors, etc.). 

The time domain criteria developed from the program requires 
the following items for HQR predictions « 

computer generated time histories of stick force (F s )> 
angle of attack (a), pilot station Z-axis acceleration 
(N Z p), and pitch rate (q) responses to a step-in, step-out 
elevator force command. 

a pencil 

a ruler 

a hand calculator 
five minutes of time 

6.5.3 Frequency Domain Criteria Considerations 

1. Low-order equivalent system results did not adequately discriminate 
between the good and bad configurations in this study. Only 38% of 
the 32 configurations of this program were predicted by level. This 
is largely due to the fact that most of the configurations had signi- 
ficant phugoid response characteristics, and these were not identified 
by the low-order models. 

2. the pitch attitude bandwidth level boundaries were revised to reflect 
transport aircraft configurations. Using these boundaries 70% of the 
58 configurations from this and the pitch rate programs were predicted 
by level. Configurations with low-frequency washout prefilters were 
not accurately predicted by this criterion. 


3. The altitude rate bandwidth criterion yielded satisfactory results, 
with 6356 of the 32 configurations of this program predicted by level. 
The characteristics of the time delay configurations in this study 
were not accurately predicted by this criterion. 

4. The Neal-Smith criterion showed the most promising results. 7336 of 52 
configurations (26 from the landing program, 22 from the pitch rate 
program, and 4 from the large aircraft program) were predicted by 
level using this technique. The criterion showed a strong correlation 
between pilot lead compensation and pilot rating, again with the 
exception of the time delay configurations. 

6.5.4 Overall Predictive Criteria Considerations 

The revised time domain criteria provided the most accurate results 
of all criteria tested. It was tested against more configurations with more 
diverse characteristics than the other criteria and has fewer limitations. 
Additionally, it can account for flight control variations and significant in 
these are time delay and sensitivity. This criteria has matured to the point 
that it can be used as a predictive tool for the flight control designer. 

The Neal-Smith criteria pitch attitude pilot lead compensation tech- 
nique showed the most promise of the frequency domain criteria, however, it 
has not been applied to a large enough data base and should be modified to 
accurately account for sensitivity effects. 

6.6 RECOMMENDATIONS 

Recommendations presented here relate to the predictive criteria 
results. Recommendations from the interpretation of the results using MIL-F- 
8785(C) as a flight control design criteria are presented in Section 5.4. 

6.6.1 Time Delay Sensitivity 

It has been shown that controller sensitivity has a strong effect on 
pilot interaction with time delay. Other characteristics in the command path 
such as feel system frequency, damping, deadband, friction, etc., may have 
similar effects. A flight program should be conducted to quantify these 
effects. 


6.6.2 Time Domain Criteria 


This program has provided quantitative data on wheel controller pitch 
sensitivity requirements for the flared landing task. Other flight programs 
should be conducted to investigate sensitivity requirements of wheel, stick, 
and side stick controllers in both the pitch and roll axes and in other tasks 
as well as flared landings. 

6.6.3 Frequency Dona in Criteria 

Future efforts in frequency domain analysis should include the 
pilot’s closed-loop gain as an integral part of any predictive criterion, in 
general, the criteria applied in these analyses did not adequately predict the 
pilot rating trends obtained from the time delay configurations. It is recom- 
mended that a frequency domain criterion to be developed that consistently 
predicts the results from such configurations. 
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